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ABSTRACT 


In 1987 a NASA panel recommended the creation of the Mission to Planet Earth. This 
mission was intended to apply to remote sensing experience of the space community to earth 
remote sensing to enhance the understanding of the climatalogical processes of our planet and to 
determine if, and to what extent, the hydrological cycle of Earth is being affected by human 
activity. One of the systems required for the mission was a wide scanning, high gain reflector 
antenna system for use in radiometric remote sensing from geostationary orbit. 

This work describes research conducted at Virginia Tech into techniques for beam 
scanning offset Cassegrain reflector antennas by subreflec.tor translation and rotation. 
Background material relevant to beam scanning antenna systems and offset Cassegrain reflector 
antenna system is presented. A test case is developed based on the background material. The 
test case is beam scanned using two geometrical optics methods of determining the optimum 
subreflector position for the desired scanned beam direction. Physical optics far-field results are 
given for the beam scanned systems. The test case system is found to be capable of beam 
scanning over a range of 35 half-power beamwidths while maintaining a 90% beam efficiency or 
50 half-power beamwidths while maintaining less than 1 dB of gain loss during scanning. 
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Chapter 1 

INTRODUCTION 

1.1 Overview 

In 1987 a committee of the National Aeronautics and Space Administration chaired by 
Sally Ride proposed the Mission to Planet Earth [1]. The concept of this long-term experiment 
is to draw upon the space community’s vast experience with remote sensing from unmanned 
inter-planetary probes to assess better the changing climate of Earth [1]. High resolution remote 
sensing is crucial to the success of this project [2]. It would be advantageous to place these 
remote sensing packages in geostationary orbits to allow real-time tracking of developing micro- 
scale weather systems such as convective cells and to decrease scene revisit time below what can 
be achieved with a reasonable constellation of low earth orbit platforms. Unlike infrared and 
visual light imaging, high resolution microwave radiometry has previously been performed from 
only low-earth orbits because very large antennas are required at geostationary orbit for good 
resolution. Additionally, a radiometric system must be able to repoint accurately its beam 
reasonably quickly to achieve the required tracking performance and scanning speed to minimize 
revisit time. This beam repointing speed requirement presents a further complication because 
the large antennas which are required generally cannot be slewed fast enough without disturbing 
both the antenna structure and its shared spacecraft bus. From 1987 through 1993 the Satellite 
Communications Group at Virginia Tech worked with the Antenna and Microwave Research 
Branch at NASA Langley Research Center to develop reflector antenna systems which achieve 
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high gains, and therefore narrow beam widths, while being capable of wide angle beam scanning 
without main reflector motion. 

1.2 Project Motivation 

The original science objectives of the Mission to Planet Earth required a 100-meter class 
reflector antenna capable of operating with high beam efficiencies over a frequency range of 10 
to 60 GHz to achieve the necessary resolution [2]. Revisions to the intended geophysical science 
applications of the radiometric experiment [3] and studies of the size/complexity characteristics 
of large space antennas [4] have since reduced the antenna size requirement to a 25 meter 
diameter projected aperture. This antenna size will produce a resolution of 20 km at the sub- 
satellite point at 18 GHz and still allow the spacecraft system to be launched by a single STS 
mission [4]. 

As mentioned above, beam scanning a reflector antenna of this size by main reflector 
slewing is impractical due to mechanical considerations. To avoid main reflector mechanical 
motion, beam scanning must be achieved by a combination of feed and/or sub-optic motions. 
Purely electrical beam scanning for a high gain, multi-band, wide scanning reflector antenna has 
been ruled out because of the immensely complex feed array which would be needed [5]. Also, 
scanning schemes which use feed motion are not feasible because even a single feed per band 
radiometer is quite heavy and extremely gain sensitive and so should not be moved. Because of 
these restrictions, the antenna must be capable of wide angle beam scanning by sub-optics 
motion. The simplest reflector antenna configurations which are capable of beam scanning by 
sub-optics motion are the dual reflector Cassegrain and Gregorian systems. From geometrical 
optics reflector antenna theory, comparable performing Cassegrain and Gregorian systems can be 
synthesized. Cassegrain reflector antenna systems are more compact than an equivalent 
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Gregorian system because the subrefiector is located between the main reflector and the focal 
point in a Cassegrain system but beyond the focal point in a Gregorian system. For this reason, 
this work focuses on Cassegrain reflector antenna systems. The need for an offset main reflector 
is the final geometry restriction imposed by the mission requirements. This requirement is 
generated by the high beam efficiency (90%) needed for accurate radiometric measurements [3). 
The desired characteristics of this wide scanning, high gain reflector system are shown in Table 
1 . 2 - 1 . 

Krichevsky and DiFonzo [6] developed a method for designing offset Cassegrain reflector 
antenna systems for use on multi-beam communications satellites. Unfortunately, as shown in 
Section 3.3, their configuration scans by feed motion and does not place any limits on reflector 
size. In 1990, Peter Foldes proposed the Type 6 reflector antenna as a compromise system to 
meet these requirements [7], This reflector antenna beam scans by tilting its small subreflector. 
Analysis of this system at Virginia Tech using the GRASP7 reflector antenna physical optics 
package suggested that its scan range would become extremely limited at higher operating 
frequencies as shown in Section 3.4. Virginia Tech began research into the scanning properties 
and optimum scanning of offset Cassegrain reflector antennas in cooperation with the Antenna 
and Microwave Research Branch at NASA Langley Research Center. This thesis presents the 
results of that research. 
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Table 1.2-1. Required characteristics of a wide scanning, high gain, geostationary 
radiometric reflector antenna system 


Operating frequency 

20 - 00 GHz 

Maximum half-power beam width (HPBW) 

0.04“ 

Beam efficiency (BE) 

>90% 

Scan range 

±7“ 

Scanning mechanism 

suboptic motion or 
simple array feed 

Aperture efficiency (f ap ) 

>70% 

Main reflector diameter (D M ) 

< 25 m. 

/ D 2 \ 

Areal efficiency | * ) 

\ d m 2 + DsV 

as high as possible 
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Chapter 2 

REFLECTOR ANTENNA SYSTEM GEOMETRY 


The electrical properties of reflector antenna systems depend primarily on the ratio of 
the focal length of the reflector to its diameter. This quantity, F/D, determines the level of edge 
taper due to geometrical effects, the required directivity of the feed antenna, and, for offset 
systems, the cross-polarization level. Further, of particular interest to this study of scan 
behavior, the scanning performance of reflector antenna systems degrades as F/D decreases for a 
given aperture size. F / D cannot be decreased to much less than unity because then the surface 
area of the reflector increases at a much greater rate than the projected aperture area. Since the 
feed must be placed at the focal point and supported by a boom, the constraints on the 
minimum F/D value require a deep structure which increases the mechanical complexity and 
weight. The mechanical depth of a reflector antenna system can be reduced by using a 
secondary reflector to form an image of the feed at the focal point of the main reflector 
paraboloidal surface. In addition, the use of folded optics allows the feed to be located in a 
more convenient position. 

The two classical types of dual reflector antenna systems are the Gregorian system and 
the Cassegrain system. In a Gregorian reflector antenna system, an ellipsoidal subreflector is 
mounted farther from the main reflector surface than the focal point of the main reflector 
paraboloid. Although this method allows a longer electrical focal length for a given main 
reflector paraboloid, the system is still quite deep since the subreflector must be mounted 
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beyond the focal point. The focal length, and therefore F/D of the main reflector surface, of this 
type of system is still limited by the requirement of mounting the subreflector beyond the focal 
point. The Cassegrain reflector antenna system, which is the system of interest here, consists of 
a hyperboloidal subreflector mounted between the main reflector surface and the focal point of 
the main reflector paraboloid. In this system, the length of the reflector antenna system is 
reduced compared to the Gregorian system since the subreflector is mounted closer to the main 
reflector surface. This difference allows the main reflector surface of a Cassegrain reflector 

antenna system to have a higher F/D value but still allows the antenna to be built with a 
relatively short support structure for the subreflector. 

This chapter explains how the Cassegrain reflector antenna system evolves from the 
prime-focus paraboloidal reflector antenna system. The concept of the equivalent prime-focus 
paraboloid is introduced for the Cassegrain reflector antenna system. A method for determining 
the geometry of an offset Cassegrain reflector antenna system which has an axi-symmetric 
equivalent paraboloid is given. Finally, the test case which will be studied is specified. 


2.1 The Prime-Focus Paraboloid Reflector Antenna System 

The prime-focus paraboloid reflector antenna system consists of a paraboloidal main 
reflector and a feed located at the focal point of the main reflector. The equation that defines 
the main reflector surface, as shown in Figure 2.1-1, is 


_ x m 2 + yM 2 


( 2 . 1 - 1 ) 


where 


F m = focal length of the paraboloid. 


The focal point of the paraboloid is located on the z-axis and the vertex of the paraboloid is 
located at the origin. In general, the physical reflector is a section of the parent paraboloid and 
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Aperture plane (A) 




Figure 2.1-1. Prime-focus paraboloidal reflector antenna system geometry. 
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can be either axi-symmetric or offset. In the case that the physical reflector is offset, the feed is 
located at the focal point on the z-axis but is tilted to evenly illuminate the physical reflector. 
A paraboloidal main reflector is used because a spherical wave leaving a point source located at 
the focal point will be transformed into a planar wavefront exiting the aperture of the antenna. 
This requirement is equivalent, in a geometrical optics (GO) sense, to the requirement that all 
rays traced from the feed point exit the aperture of the antenna system perpendicular to the 
aperture plane after reflection from the paraboloidal surface [1], The correct main reflector 
shape is verified by substituting the incident and reflected ray unit vectors. ? 2 and f„ at the 
main reflector and the unit normal of the surface, n M , into Snell’s law. 


For the prime-focus reflector antenna system shown in Figure 2.1-1, a ray from the feed 
point, {0, 0, F m }, to any point on the reflector surface, |x M , y M , z M J, is 

*2 = { X M *> y ’ ( z M' F m) z }- (2.1-2) 

Since the length of this ray is 

|F 2l = V^ X M 2 + >M 2 + ( z M* F m) 2 • (2.1-3) 

the unit vector of a ray from the feed point to any point on the reflector surface is 

r 2 - (\/ x M + y M + ( z M- F m ) 2 ) { X M X > y\iy. ( z M- F m) z }• (2.1-4) 

The unit vector of a ray exiting the aperture in the z-direction is 


The required surface shape can be verified by substituting (2.1-4) and (2.1-5) into Snell’s 
law. Snell’s law, expressed in vector notation, is 


r, - r 2 - 2 (r 2 -n M )n M 

where 


( 2 . 1 - 6 ) 


ff M - the inward unit normal of the main reflector paraboloid at the point of reflection. 
An inward normal vector at the paraboloidal main reflector surface is 
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- .(Aij 
n M - \~dT x 

dz ™A 

’ dz y ' dz z y 

(2.1-7) 

For the paraboloidal surface given by (2.1-1) 


-d x M - X M 

ST-TF^ 


(2.1-8c) 

-^M_ - 
dz 2F m 


(2.1-8c) 

5z m _ , 

nr ~ 1 


(2.1-8c) 

so 



- _ r ‘ x m - 

nM -K J 

. -y M - A 
2^ y ’ Z | 

(2.1-9) 

Since the length of this normal vector is 


In- 

! + y M 2 + 4 F m 2 

(2.1-10) 

1 n M I ~ 


the unit vector normal to the paraboloidal surface is 


= (v x M 2 

+ yM 2 + 4F M 2 ) {- x M x --yMy* 2 F m z }- 

(2.1-11) 


Substituting the incident ray unit vector (2.1-4), the reflected ray unit vector (2.1-5), and the 
inward normal unit vector at the point of reflection (2.1-11) into Snell’s law (2.1-6) and 
expressing the result as equations for each of the x, y, z components gives: 


0 = 


0 = 


1 = 


X M 


\/ x M 2 + y\f 2 + ( z M -F m) 

\/ x m 2 + y\i 2 + ( z m _f m ) 2 
( f m' z m) 


_-2(f 2 -n M ) 


- X M 


\/ x m 2 + yM 2 + 4F 


M 


\/ x M 2 + >M 2 + ( z M -F m) J 


“ 2 ( F 2 - ) 


2 ( f 2 ■ ) 


~ I'm 


\/ x M 2 + yM 2 + 4F 1 
2F 


M 


r M 


\/ x M 2 + + 4 F m" 


where 


( F 2 ‘ ) = (\/ x M 2 + yM 2 + ( z M _F m) 2 \/ x M 2 + yM 2 + 4F M 2 ) 

[- X M 2 ~ I'M 2 + 2F m( z M _F m)]- 


(2. 1-1 2a) 
(2. 1- 12b) 
(2.1-12c) 


(2.1-13) 

Equation (2.1-12) is simplified by manipulating the radical in the denominator of the r 2 terms 
into the form of the n M terms. Expanding the radical in f 2 arK ^ substituting (2.1-1) for z M into 
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(2.1-12) leads to 


V^M 1 + Jm 1 + Um-Fm) 2 = 


N 


,,, 2 , 2 , [ X M 2 + - V M 2 p V 

M +>M \ 4F m - F mJ 


(2.1-14) 


4F m F M 2x M 2+16 f ’m 2 >’m 2 +( x m 2 +>'m 2 - 4 77 ) (2.1-15) 

= 4T^ v/ ,6F M 2x M 2 + I«F M a y M * + x M 4 + 2 x M 2 y M 2 

- 8F M 2 *M 2 - 8F M 2 yM 2 + >M 4 + 16F m “ (2.1-16) 


~ 4 F m \/ 8F M 2x M 2 + 8F M 2 y\1 2 + X M 4 + *- x M 2 > M 2 


+ y M 4 + i6 f m 4 


(2.1-17) 


= 4T^ \/( Xm2 + yM * + 4 Fm2 ) 2 (2.1-18) 

which is the form of the radical in the n terms of (2.1-12). Substituting (2.1-18) into (2.1-12) 
and (2.1-13) gives 


where 


ft _ 4F M x M 

■ - 2 ( f i • n ) 

X M 2 + yM 2 + 4F M 2 

0 _ 4F MyM 

- 2 ( F; - n ) 

x m 2 + >m 2 + 4F m 2 

4 F M ( F M ‘ Z M ) 

- 2 ( r j • n ) ■ 

X M 2 + yM 2 + 4F M 2 

0 

( r 2 ’ n M ) 4F M / „ 

- x m 2 - yM 2 


V^ x M 2 + - V M 2 + 4F m 2 


- 


\/ x m 2 + y\i 2 + 4 f m 2 

- f m 

V^ X M 2 + yM 2 + 4 f m 2 


(2.1-19a) 

(2.1-19b) 

(2.1-19c) 


( 2 . 1 - 20 ) 


Defining 

= X M 2 + ^M 2 + 4 F M 2 

and using C m (2.1-19) with (2.1-20), yields the following equations that correspond to the 
Cartesian component equations of (2.1-12) 


( 2 . 1 - 21 ) 
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o = V - + 2 F m (z„.F m )]x m 

(2. l-22a) 

0 = 4F My M + 8lM[. Xm 2 - y M * + 2 F m (z m -F m )] y M 

(2.1-22a) 

1 - 4F m( 2 M- F m) + ! F m[. Xm 2 . y M 2 + 2 F m (z m - F m )] 2 F m . 

(2. l-22a) 


Substituting (2.1-1) into the bracketed term common to (2.1-22a), (2.1-22b), and (2.1-22c) gives 
[' = [" X M 2 ‘ +- F M Z M'“ F M ] 

= i[-2x M 2 - 2y M 2 + x M 2 + y M 2 - 4 F M 2 ] 

= y [ x m 2 + >m 2 + 4Fm 2 ] = y- 


(2.1-23) 

(2.1-24) 

(2.1-25) 


The satisfaction of the x and y components of (2.1-12) by the paraboloidal main reflector shape 


is verified by substituting (2.1-25) back into (2.1-22) 

„ _ 4F M X M 4F M x M 

°~~~ < < 

4F M y M 4F MyM 

r - ' < 


1 = 


4F m (z m -F m )^ 8F m 2 


(2.1-26a) 

(2.1-26b) 

(2.1-26c) 


< < ' 

The satisfaction of the z component of (2.1-12), shown in (2.1-26c), is verified by substituting 


(2.1-1) into (2.1-26c). This substitution and simplification gives 

i = liM (Sm+ f m > = ( SrJ 1 " 1 + F «) 


(2.1-27) 


4F m ( x m 2 + y» 

~V 4 


yM + 4F m I. C 


" M 


c 


(2.1-28) 

The equalities of (2.1-26a), (2.1-26b), and (2.1-28), derived by enforcing Snell’s law at the 
reflector surface, prove that the paraboloid focuses the rays from a point source feed at the focal 


point to infinity as a prime-focus reflector antenna system. Although this method of GO 
raytracing can be used for synthesis, it is desirable to have a simpler synthesis method that does 
not require taking derivatives at the surface. 
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2.2 The Levi-Civita Theorem 


The proper reflector surface shape for a prime-focus antenna system can also be 
determined by imposing an equal path length constraint on each ray traced from a feed point at 
the focus to the reflector surface and leaving the reflector antenna system in a direction 
perpendicular to the aperture plane. The Levi-Civita theorem [1] states that, for a set of rays 
which exit a single reflector antenna system mutually parallel and strike a reference plane, the 
following two conditions are equivalent: 

(1) The path length of every ray from the source to the reference plane must be equal. 

(2) Snell’s law is satisfied at the reflector surface. 

Since the rays in a focused prime-focus reflector system exit the system mutually parallel, 
present a planar phase front in the aperture, and undergo a single reflection in the antenna 
system, this theorem can be used to determine the proper shape for its reflector surface [1], 

The reflector shape of a prime-focus paraboloidal reflector antenna system can now be 
found by determining the path lengths for rays traced from the feed point to the aperture after 
reflection from the reflector surface. The reflector is assumed to have a feed point at z = F on 
the z-axis and to have a parent surface which passes through the coordinate system origin as 
shown in Figure 2.1.1. As in Section 2.1, where the shape was verified by enforcing Snell’s law 
at the reflector surface, the desired direction of propagation is along the z-axis. Unlike in 
Section 2.1, the correct reflector shape will be found directly rather than verified. The length of 
any ray from any point on the reflector surface to the aperture plane is 

l R i|= z A* z M ( 2 . 2 - 1 ) 

where 

Z A = the location of a point in the aperture plane along the z-axis 
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{ X M’ z m} = a P°' nt on reflector surface. 

As shown in Figure 2.1-1, the distance from the point feed at the focus. {0, 0, F M }, to any point 
on the reflector surface is 

| R-a | = \/ x m 2 + ^m 2 + ( z m* f m ) 2 • (-•-•-) 

To satisfy condition (1) of the Levi-Civita theorem 

I Rj I + I R 2 1 = Z A - Z M + \J X M 2 + ^M 2 + ( z M‘ f m) 2 = L p (2.2-3) 

where 


L p = the constant path length for all rays. 
Rewriting (2.2-3) and squaring both sides 

X M 2 + yM 2 + ( z M‘ F m ) 2 = L p ' Z A + Z M- 
X M 2 + >M 2 + ( Z M " ^m ) 2 = ( L p ' Z A + z m) 2 - 


(2.2-4) 

(2.2-5) 


After expanding the squared terms, (2.2-5) results in 

X M 2 + J^M 2 + Z M 2 * 2F M Z M + ^M 2 = ^p 2 ‘ 2F p z A + *“^p z M + Z A 2 * “ Z A Z M + z M 2 { 2 - 2 ‘6) 
Evaluating (2.2-5) at the coordinate origin (0, 0, 0) gives 

L P = Z A + F M- t 2 - 2 ' 7 ) 

Substituting this expression for the constant L p into (2.2-6) gives 

X M 2 + >M 2 + Z M 2 ' 2 F M z M + F m 2 = Z A 2 + 2 F M Z A + FM 2 ‘ 2 z a 2 * 2 F m z a + 2 z A z M + 


2 F m z M + Z A - 2 z A z M + Z M 


A ' “ * M 4 A t “A 4 M 
2 


( 2 . 2 - 8 ) 


After simplification, (2.2-8) becomes 


X M 2 + I'M 2 _ 

TTS 


(2.2-9) 


This is the equation for a paraboloid as given by (2.1-1). This establishes the paraboloid as the 
proper shape for the reflector of a prime-focus reflector antenna system and shows that this 
shape can be developed from either Snell’s law or the equal path length constraint. The 
simplicity of the equal path length condition compared to the verification using Snell’s law 
presented in Section 2.1 is striking and is of great value. It suggests that the equal path length 
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should be used when examining the more complex dual reflector antenna system. 


2.3 The Canonical Cassegrain Reflector Antenna System 


The canonical Cassegrain reflector antenna system consists of a paraboloidal main 
reflector surface, a hyperboloidal subreflector surface mounted with one of its focal points 
coinciding with the focal point of the main reflector surface shape, and a feed located at the 
second focal point of the hyperboloidal subreflector. The equation which defines the paraboloid 


of revolution that forms the main reflector surface, as shown in Figure 2.3-1, is given by (2.1-1) 

The general second order equation which defines the subreflector surface, as shown in Figure 2.3- 
I, is 


Ax s 2 + Bx s y s + Cy s 2 + Dx s + Ey s + F 

where 


G Zg 2 + H Zg + I Xg Z s + J y s z s 


(2.3-1) 


A = 4(/? 2 -a 2 ) 
B = 


F = 0 2 -4ft 2 (x n 2 + y f2 2 + f2 2j 

G = 4 (a 2 - f> 2 ) 


C = 4( 7 2 -a 2 ) 

D = 4(/?0 + 2a 2 x f2 ) 

E = 4 (yO + 2o 2 y f2 ) 

= vj 7 + ( y n- y f2 7 + ( z n - Hi f 


H — 4 ^60+ 2ft 2 zj. 2 j 
I = -8 (36 
J = -8y6 


^ - ( X f 2 - X rx) j = (y fJ -y n ) * = ( Zf2 -z fI ) 

9 ~ X fl 2 + y f | 2 + z fl 2 - x f2 2 - y f2 2 - Zf2 2 . ft 2 
fl, f2 = the two focal points of the hyperboloid, 
e = the eccentricity of the hyperboloid. 

Th* details of the derivation of this equation have been omitted due to thei, length. I„ , he 
Cassegrain reflector antenna system, the hyperboloidal suhreflecto, surface is used because the 
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Figure 2.3-1. Cassegrain reflector antenna system geometry. 
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focal points of the hyperboloid are conjugate points. That is. an image of the feed placed at the 
second focal point of the hyperboloid is formed at the first focal point of the hyperboloid which 
is located at the focal point of the paraboloidal main reflector surface. The subreflector surface 
is chosen to be on the branch of the hyperboloid which is closer to the focal point of the main 
reflector. As with the prime-focus paraboloidal reflector antenna system, in general, the 
Cassegrain reflector antenna system can be either axi-symmetric or offset. In either case, the 
physical main reflector and subreflector are sections of the surface of the parent paraboloid and 
hyperboloid, respectively. 

The imaging property of the hyperboloidal subreflector can be simply shown through 
the use of an extension of the Levi-Civita theorem discussed in Section 2.2. Xianzhong [3] has 
extended the Levi-Civita theorem to include arbitrary focused dual reflector antenna systems. 
Xianzhong starts by assuming a set of rays which emanate from a point source at the focal point 
and reflect through the dual reflector system shown in Figure 2.3-2 [3], The rays exit the 

aperture mutually parallel and strike a perpendicular reference plane. Xianzhong considers three 
restrictions on the set of rays: 

(1) The path length of every ray from the source to the reference plane must be equal. 

(2) Snell s law is satisfied at the main reflector surface. 

(3) Snell s law is satisfied at the subreflector surface. 

If the rays satisfy any two of these conditions, then the remaining restriction is satisfied 
automatically [3]. The proof uses an arbitrary dual reflector system including an aperture plane, 

A, a main reflector, M, and a subreflector, S. The equations for these surfaces then are [3] 

z s = z s(x s , y s ) (2.3-2) 

Z M = z m( x M’ >m) (2.3-3) 

z M = * a (where z a >[z M ] max ). (2.3-4) 

where x, y, and z are the coordinates of the intersection of the ray with the surface indicated by 
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Main reflector (M) 


Aperture plane (A) 



Figure 2.3-2. General doubly reflected ray geometry as used by Xianzhong [3] 
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the subscript and z a is the z-coordinate of the intersection of the ray with the aperture plane. 
The z-axis is taken to be normal to the reference plane and the rays which strike the plane are 
parallel to the z-axis [3]. The unit vectors r 3 , r 2 , and f, are the unit vectors of one of the rays 
which is reflected through the system as indicated in Figure 2.3-2 [3], If Snell's law is enforced 


at the subreflector and main reflector, then 


(*s\ 

(cosa x - cos/? x ) 


\**s) 

(cos P z - COSQ z ) 

(2.3-5a) 

(to. s \ 

_ (COSOy - COS^y) 


vW 

(cos/? z - cosrt z ) 

(2.3-5b) 

. the subreflector and 


(to m\ 

cosp x 


\toj 

(1 - COS0 Z ) 

(2.3-6a) 

( 9z m\ 

COS^y 


UiJ 

~ (1 - C0Sa z) 

(2.3-6b) 

the main reflector (3j. The constant path length through the system, L , is 

J 


w 

X S + ys 2 + z s 2 + y ( x m* x s) 2 + (y M -y s ) 2 + (* M -*s) 2 + (v* M ) • 

(2.3-7) 


Xianzhong solves the system by taking the partial derivatives of (2.3-7) considering 
either x M and y M or x s and y s to be the independent variables in the system. If the main 
reflector parameters, x M and y s , are considered independent, then 

(ai) = (^)[ (COSOjt ■ co ^x)- («»/?, - cosa,)(|!|)] 

+ (fe) [ (cosft * * cos 0 y ) ‘ ( cos ^ ■ COSft *) (ij)] 

+ [ cos/? x - (1 - cos/?,) (^)] (2 .3. 8a ) 

fe) = fe) [ (COSCVx ‘ COS0 ^ ■ ( C ° S/? - - C ° Sa ^ (if)] 

+ (fe) [ (cosfty ' cos i 0 y) ‘ ( f0s ^ * C0Srt *) (if )] 
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(2.3-8b) 


(2.3-9a) 


COSiJy - (1 COS^) (^~|)]. 

If the subreflector parameters, x s and y s , are considered independent, then 

= [ (C “°" ' " OS0 ^ ' * C05 ^‘ ' C ° Sa ‘ > (ft)] 

ft^)['“ V(l '““‘ , * ) (ftl)] 

+ (%)[“*-<' --*>(&)] 

(®ys) ~ [^ cosay ’ " ( cot0, ' c ““’) (ft)] 

+ (^)K- (, -^>(^)] 

+ (^)H''°‘“*' , * ) (^)]' 

If the Snell’s law requirements, (2.3-5) and (2,3-6), are satisfied and substituted into (2.3-9), 
then 

(ft)=» ( ' 2 - 3 - io “» 


(2.3-9b) 


a- 


(2.3-10b) 

so the minimum path length constraint of Fermat’s principle is satisfied and the first restriction 
that Xianzhong considers is proven [3]. Alternately, substituting (2.3-5) and (2.3-6) into (2.3-8) 
gives 


dh 


P - 


dx 


= 0 


M , 


(2.3-1 la) 


(fe)- 


(2.3-1 lb) 

The other possible combinations of imposed and automatically satisfied restrictions are also 
considered by Xianzhong and are stated to be correct [3]. 

This extension of the Levi-Civita theorem makes the analysis or synthesis of the 
geometry of the canonical Cassegrain reflector antenna system significantly simpler. The Levi- 
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Cvu, theorem and the anai.vsir of Section 2 ,, or of ^ , ,, ^ ^ 

U» U satisfied at the -ta. sorface for „ y ra} . ^ pri _ rofus pitab<)|oidi| 

ranac. 0 , antenna s y stem. Beeans, th, main refieeto, surface i„ . canonic, Cassegrain ^ 
antenna spstem in par.hoioida, and the hpp.rboioidai snbrefiector p,„d„ces . spheric, phast 
mont centered on iu firs. fncnn when minated h y a point nonrce at i„ necond focnn the 

teani. o, the precedin, an., y sis of the prime-foen, parahoioida, refiector an.™,,, n.vntetn can he 

used to assure the satisfaction of condition f2) of th* v i 

( ’ f th X,a, >zl»ong extension of the Levi-Civita 

theorem for dual reflector antenna systems. 

Condition of thin ea.ension can he shown to he satisfied for the h yP erho,oida, 
nnbrenecto, hecanne the defi„i„ g rharacterintic of the h.perhoi. in that the difference the 
distances frotn a„ y point on the h y perho,a the foci in constant. A distinction in made 
between th. path ,„ g th of a r. y in the P ,i m . focnn s y stem, Lpp[ „ and the path ,e„ g th of a ,. y 

" ‘ he C “ S ' 8r ‘ in SWm ’ hpperboi, shown in Fig „ r e ,,, 

defining characteristic ran be written as 

d»-d n =V 

where 


d|I< dfr = the distances from f, and f> to a point 0 „ the snbrellecto,, 
fl, f2 = the focal points of the hyperboloid, 

2a = the distance between the foci of the hyperbola, 
e = the eccentricity of the hyperbola. 

From (2.2-3) 

N + ISrMp,,,. 

where 


(2.3-13) 


V = the T ';t ° r - - -on. spate, which in rcfiected to 

the feed point. 
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A hyperboloid.! subredccfo, located in the redeem, system with one focus .t the focn. of the 
main reflector paraboloid and the other foes at the desired feed point will reduce the length of 
an, ray reflected from the hyperboloid.! surface by the distance d„, the distance between the 
point on the subreflecto, and focal point fi, and increase the length of an, ray redacted from the 
hyperboloid.! surface by the distance d„, the distance between the point on the subreflecto, and 
focal point f2, as shown in Figure 2.3-1. Therefore, the length of any ray which reflects through 

the antenna system is 

_ 2a (2.3-14) 

| R| | + | Ft’2 | + ^f2 * d n = kpcass. ^ppf- ^ e 

Since the analysis of S«tion 2.2 shows that L pp ,., the path length of any ray in the prime-focus 
reflector antenna system, is constant, (2.3-14) shows that L_. is constant from the definition 
of the hyperbola. Thus, condition (1) of the Xianshong extension of the Levi-Civita theorem is 
satisfied. Condition (3) is also now known to be satisfied since all of the other conditions have 
Seen shown to be satisfied. Avoiding direct application of Snell's law at the subredecto, surface 
through the use of this theorem is advantageous it is more difficult than the calculation of the 

total path length. 

2.4 The Equivalent Paraboloid 

The equivalent paraboloid of a dual reflector antenna system is a mathematically 
constructed prime-focus paraboloidal reflector antenna which produces the same aperture field 
distribution as the dual reflector antenna system. The equivalent paraboloid has been found to 
be an accurate predictor of the electromagnetic performance of the modeled dual reflector 
antenna system but with the advantage of requiring far less computational effort [4]. Rusch et 
al. [4] developed a general definition of the equivalent paraboloid for offset Cassegrain and 
Gregorian dual reflector antenna systems with arbitrary aperture shape. The development 
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presented here will concentrate on the case of the canonical Cassegrain reflector antenna system. 
The notation of Rusch et al. [4] will be used through the end of this chapter to maintain 
consistency with the referenced article. 


Figure 2.4-1 shows the geometry and definition of variables which Rusch et al. [4] use in 
the derivation of the general equation of the equivalent paraboloid of a Cassegrain reflector 
system. Rusch et al. [4] use four Cartesian coordinate systems: { Xp , y p , Zp } and {Xj , ^ I# } 
which have their origin at the focal point of the main reflector surface and { x f , y f , Zf ) and { X/J , 
y 0 , z 0 ) which have their origin at the second focal point of the subreflector surface. Each 
coordinate system also has an associated spherical coordinate system. The main reflector surface 

is paraboloidal and the subreflector surface is hyperboloidal. The main reflector surface is given 
by 

_ 2F 

Pp 1 + cos 0 p (2.4-1) 

where 


F _ the focal length of the paraboloidal main reflector surface. 
The subreflector is defined by 


pf - p s = ¥ 

where 

2c = the interfocal distance of the hyperboloid, 
e = the eccentricity of the hyperboloid. 

The formulas 


Pi = 


_ c(e 2 - 1 ) 
e(ecos O 0 - 1 ) 


, e 0 | e - Ilf #1-1 

Un T = 7TT[ tan fJ • 

derived from (2.4-2), and the identity 


, » - 6 f $ l-i 

tan— ^—5 = tan-^ 


(2.4-2) 


(2.4-3) 

(2.4-4) 

(2.4-5) 
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PARABOLOIO 



Figure 2.4-1. Geometry definition of the general Cassegrain reflector antenna for the 
Rusch et al. development of the equivalent paraboloid [4]. 
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are useful in determining the geometrical optics aperture field due to a point source located at 
the second focus of the hyperboloidal subreflector. 


Rusch et al. [4] now find the aperture field at a point A due to a spherical wave 
emanating from the feed point. The magnitude of this field at a point in the aperture is given 

by 



(2.4-6) 

where 


<t>f) — the field pattern of the feed expressed in the spherical feed coordinate system. 
Rusch et al. [4] develop the amplitude term as a function of 0 f and <*> ( using the formulas 

costfp — COS0 S cos/? - sin 0 s co s<^ g sin/? 

(2.4-7) 

sin(? s = jr sin 00 

(2.4-8) 

Pj cos 00 - 2 c 
cos0 s = ^ 

(2.4-9) 

and the identities 


sin Op cos = sin Of cos <j> f cos a + cos Of sin a 

(2.4-10) 

cos Op = -sin Of cos <f>f sin a + cos Of cos a. 

(2.4-11) 

After substitution, they show that 


P s _ 1 f t 2c(l - ecos/?) _ . v 1 

PfPp epf + s »n^ f cos</> f sin(a + /?) - cos0 f cos(a + /?) 1. 

(2.4-12) 

Since the angle n has not been used yet, it can be chosen arbitrarily. Rusch et al. 

[4] set o so 

that the z r axis is aligned with the axis of the equivalent paraboloid. This selection 

causes (2.4- 

12) to reduce to 


P s 1 + COS0 f 

P{P p “ 2F eq 

(2.4-13) 

where 


F - F 1 e2 ' 1 1 

^ (e 2 + l) - 2ecos/? 

(2.4-14) 
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So (2.4-6) now has the form 


E a = 


E(*f. *r) 


(2.4-15) 


eq 


which is the aperture field distribution of a prime-focus paraboloidal reflector antenna system 
with a focal length of F eq . Finally, (2.4-12) can be solved to determine the required feed tilt 
angle, ft: 

(e 2 - 1 ) sin 0 (2.4-16) 


sin a = 


cos ft = 


(e 2 + l)- 2 e cos 0 
(e 2 + 1 ) cos 0 - 2 e 


(e 2 + l) - 2ecos/? 

These equations can be combined to form 
(e 2 - 1 ) sin 0 


tan ft = 


(e 2 -h 1 ) cos/? - 2e 


(2.4-17) 


(2.4-18) 


or 


tan §■ = e — } tan ^ . 
2 e - 1 1 


(2.4-19) 

Either (2.4-18) or (2.4-19) can be used to find the required feed tilt for an equivalent paraboloid 
of focal length F eq which will model the Cassegrain reflector antenna system. 


Rusch et al. [4] caution that this model is based on geometrical optics principles applied 
to a focused reflector antenna system. Although defocusing and diffraction effects are ignored, 
the equivalent paraboloid should provide a reasonably accurate prediction of the near-boresight 
aperture and radiated fields of the focused dual reflector system it models. In order to minimize 
these effects, Rusch et al. [4] recommend that the equivalent paraboloid model not be applied to 
systems which have a subrefiector smaller than approximately 10 wavelengths or an edge 
illumination taper of less than 10 dB on the subreflector. 


To demonstrate the accuracy of the equivalent paraboloid as a model of a focused dual 
reflector system which satisfies the restrictions of (2.4-14) and (2.4-16) through (2.4-19), Rusch 
et al. [4] considered the system shown in Figure 2.4-2. The far-field electromagnetic 
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Figure 2.4-2. Geometry definition of a Cassegrain reflector antenna system with a 
circular aperture as used by Rusch et al. [4], 
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characteristics of this system and its equivalent paraboloid were calculated using physical 
optics/surface integration. The characteristics of this system are shown in Table 2.4-1. Figure 
2.4-3 [4] shows the principal and 45’ plane far-field radiation patterns for the Cassegrain system 
and its equivalent paraboloid. The equivalent paraboloid has a gain about 0.5 dB greater than 
the Cassegrain system but otherwise predicts the co-polarized antenna pattern of the Cassegrain 
system fairly accurately down to as low as 30 dB below the peak. The accuracy of the cross- 
polarization level prediction by the equivalent paraboloid is not discussed by Rusch et al. [4]. 
Rusch et al. [4] attribute the differences between the co-polarized patterns of the systems and 
the cross-polarized component of the Cassegrain system to diffraction effects and spillover at the 
subreflector surface. 

Because the equivalent paraboloid of a dual reflector antenna system is derived under 
the assumption of a focused system, the equivalent paraboloid is a poor predictor of scan 
performance. Figure 2.4-4 [4] shows how the accuracy of the performance of the Cassegrain 
system predicted from the equivalent paraboloid degrades with beam scanning. 

2.5 The Minimum Cross- Polarization/Spil lover Condition for Offset Dual Reflector Antenna 
Systems 

The second reflector surface in a dual reflector antenna affords the antenna designer 
additional degrees of freedom in the synthesis of the antenna system. This can be of value 
during antenna synthesis to improve the performance of the antenna system. For instance, 
because the main reflector and subreflector can be shaped in a dual reflector system, the designer 
can develop antenna systems ranging from a spherical main reflector with a phase correcting 
subreflector for wide-scanning to gain-optimized systems where the subreflector and main 
reflector are shaped to give a desired aperture distribution. The equivalent paraboloid concept 
suggests a fundamental configuration selection possibility. If the equivalent paraboloid for a 
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Table 2.4-1. Reflector 


geometry used by Rusoh et al. [4] to verify the accuracy of the equivalent 
paraboloid far-field pattern predictions. 


Main reflector focal length (F) 

62.5 A 

Main reflector diameter (D) 

100 A 

Main reflector offset height (d c ) 

75 A 

Feed tilt angle (o) 

26.64° 

Subreflector axis tilt angle (/?) 

9° 

Subreflector eccentricity (e) 

1.996 

Subreflector interfocal distance (2c) 

20.41 A 

Equivalent paraboloid focal length (F ) 

179.13 A 


Feed type 

circularly symmetric 
cos q with 10 dB edge 
illumination 
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Figure 2 4-3. Equivalent paraboloid physical optics far-field pattern analysis results for 

the Rusch et al. system of Table 2.4-1. Principal and 45“ plane pattern cuts 
for the Cassegrain reflector antenna system and its equivalent 
paraboloid [4]. 

3C 
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Figure 2.4-4. Equivalent paraboloid physical optics analysis results for the Rusch et al. 

system of Table 2.4-1. Pattern cuts in the plane of scan for a Cassegrain 
reflector antenna system and its equivalent paraboloid scanned in the plane 
of offset [4j. 
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dual reflector antenna system with a circular aperture is axi-symmetric, the cross polarized field 
component radiated by the antenna due to geometrical offset effects and the antenna efficiency 
reduction to spillover effects will be simultaneously minimized [4]. 


Rusch et al. [4] developed a technique based on the equivalent paraboloid concept for 
simultaneously minimizing feed power spillover and the cross- polarized field component. Figure 
2.4-2 [4] shows the geometry used in this development. The most significant change from the 
geometry used in the development of the equivalent paraboloid is the introduction of D cq , the 
diameter of the equivalent paraboloid, and d c , the offset distance of the equivalent paraboloid. 

cq 

Because the development of the equivalent paraboloid, as presented above, is intended to be 
valid for any ray which reflects from the feed point through the dual reflector system, the 
equivalent paraboloid is shown to exist but its dimensions are not determined. Specializing this 
development to dual reflector systems which have circular apertures, allowed Rusch et al. [4] to 
show that D = D cq and also find d D , the offset height of the equivalent paraboloid. 


cq 


Rusch et al. [4] use (2.4-4) , in the y = 0 plane, to write 

1-1 


tan 




(2.5-1) 


By expanding the tangent terms and using (2.4-19), they showed that 

Of e 2 + 1 - 2e cos/? / $\ 2e • * c 0 , 

tan 2= j^TTf \ n 2/ ' ~i sm0 - (2 - 5 ‘ 2) 

Rusch et al. [4] used clockwise angles as negative and use 0 U and 0 . to represent the angle 

'eq 

between the upper and lower equivalent paraboloid edges and the z^axis. These conditions 
lead to: 


0 „ 


Or 


D eq = 2F eq 


tan — ^ - tan — ^ 

. l Li 


d % = F «. 


tan — 5^2 -f- tan — ^ 

. z 1 j 


(2.5-3) 

(2.5-4) 
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Substituting (2.4-14) and (2.5-2) into (2.5-3) and (2.5-4) gives 
D eq = 2F(tan^ - tan ^) = D 

and 


do =(d 0 -2F 0 \ 

"" V **+l - 2 e cos /? ) (-.5-6) 

where and fl u are the angles between the lower and upper edges of the physical paraboloid of 
the Cassegrain system and the z-axis as shown in Figure 2.4-2 [4J. Setting d 0 ^= 0 allows the 
solution of (2.5-6) for the required 0 for an axi-symmetric equivalent paraboloid: 

> a ngc- 2esin/? 

1 c + 1 - 2 e cos 0 (2.5-7) 


wher. the Angle #* as shown in Fig,„e 2.4-2 [4], is .he angle between .he center of .he aperture 
projected on.o .he main r.ftato, surface and the -z-axis. In general, radiated eross-polarized 
fields can he minimized by aligning the feed along the axis of the equivalent paraboloid |4J. The 
Rusch condition, (2.5-7), forces the axis of the equivalent paraboloid to align with the ,„g„|., 
center of the subreflector. This alignment will simultaneously minimize the radiated cross- 
polarized field component and the spillover [4) and will give the offset CWgr.in reflector 


antenna system unscanned characteristics similar 


to those of the axi-symmetric prime-focus 


paraboloidal reflector antenna system. 
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Chapter 3 

BEAM SCANNING IN CASSEGRAIN REFLECTOR ANTENNA SYSTEMS 

The radiation characteristics of a Cassegrain reflector antenna system are degraded 
when its main beam is scanned away from the focused boresight direction. The gain of the 
antenna and the symmetry of the pattern are reduced while the beamwidth and the sidelobe 
levels are increased. The factors which lead to this degradation must be minimized to achieve 
good performance during even moderate beam scanning. This chapter identifies the major 
factors which contribute to the degradation of antenna performance during beam scanning and 
presents two previous methods of beam scanning for the offset Cassegrain reflector antenna 
system. This material offers insight into the beam scanning properties of the offset Cassegrain 
reflector antenna system and facilitates the development of the simplified error functionals in 
Chapter 5. 

3*1 Principles of Beam Scanning 

Beam scanning in any antenna system requires a tilted a|>erture phase surface [1]. This 
phase taper can be achieved either mechanically by moving the antenna feed system or 
electrically by imposing a phase taper with an array feed. This work is intended to determine 

the optimal mechanical motions for beam scanning in an offset Cassegrain reflector antenna 
system. 

To create an asymmetric phase distribution across the aperture of a reflector antenna, 
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the antenna system feed phase center must be moved laterally away from the unsranned feed 
point. In the offset Cassegrain reflector antenna system this movement can be achieved either 
by physically moving the feed or by repositioning the subreflector to create a virtual feed 
movement. While any asymmetric aperture phase taper is sufficient for beam scanning, the 
formation of a clean scanned beam requires that two sources of beam degradation be minimized; 
these are discussed in the remainder of this section. 

3.1.1 Illumination Error Effects 

Aperture illumination amplitude error is a major source of antenna performance 
degradation during beam scanning. In an offset Cassegrain reflector antenna system, the 
significant source of illumination amplitude error during beam scanning is misalignment of the 
illumination pattern with the main reflector caused by the lateral feed translation. A small 
amount of illumination error is also introduced during beam scanning by edge illumination taper 
changes resulting from axial feed translation. Illumination misaligment, illustrated in Figure 
3. 1.1-1, can be significant in subreflector scanned offset Cassegrain reflector antenna systems 
because the required subreflector tilt repoints the central ray from the feed away from the center 
of the main reflector. 

The effects of the illumination misalignment error can be estimated for an offset 
Cassegrain reflector antenna system by inducing feed mispointing in the equivalent paraboloid 
since feed mispointing does not induce defocusing. In order to estimate the effects of 
illumination misalignment on the far-field pattern on a Cassegrain reflector antenna, uv-plane 
patterns were generated for its equivalent paraboloid using the TICRA GRASP7 numerical 
electromagnetics code using physical optics surface integration analysis. As will be discussed in 
Chapter 4, the equivalent paraboloid of the test case is axi-symmetric and has a diameter, D eq , 
of 10.63 meters and a focal length, F eq , of 42.48 meters. Viewed from the feed, the equivalent 
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paraboloid subtends a half-angle of 7.159'. For this analysis the reflector was fed by standard 
Gaussian feed patterns with reflector edge illuminations of -5, -10, -15, -20, and -25 dB. The 
edge illumination specified for the study of illumination misalignment does not include the 
spherical spreading loss edge taper, but this efTect is less than 0.04 dB. Scan induced 
illumination misalignment was simulated by pointing the center of the illumination away from 
the center of the reflector in steps of 0.2 reflector radii, R eq , up to a maximum of 2 R eq . The 
maximum illumination misalignment of 2R,*, was chosen because it causes the feed pattern to 
just miss the reflector surface. 


The actual gain of an aperture antenna, G, is given by 

4* A e 4 *a (3.1. 1-1) 

C * " *p A 2 

where A is the area of the projection of the main reflector into the aperture plane, and A is the 
wavelength of operation. The effective aperture area, A e = e^A, where e ap is the aperture 
efficiency of the antenna. Aperture efficiency can be factored as: 

(3. 1.1-2) 

£ ap ~ ~ e sp e &mp 6 <t> 

where £ is the spillover efficiency of the reflector antenna system, £ in is the illumination 
efficiency of the reflector antenna, e UBp is the illumination amplitude efficiency of the reflector 
antenna, and £ ^ is the illumination phase efficiency of the reflector antenna. Spillover efficiency 


is defined by 



where Pj(r,0) is the power density incident on the reflector aperture and P t is the total power 
from the reflector feed antenna. Illumination amplitude efficiency is defined by [2] 



( 3 . 1 . 1 - 4 ) 
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Illumination phase efficiency, e ., is given by 




II 

Aperture 


PiM) e j* dA 


1 1 


») 


J *|2 


dA 


ainp 


(3.1. 1-5) 


aperture 


where P ; (r, 6) e j * is the complex power density incident on the reflector aperture and f is the 
illumination amplitude efficiency found from (3.1. 1-4). The illumination amplitude efficiency 
term is included in illumination phase efficiency to remove the effects of aperture power 
illumination so that the phase error efficiency is 100% for a focused system. 


The far-field radiation patterns of the equivalent paraboloid with illumination 
misalignment were calculated using GRASP7 physical optics/surface integration at 20 GHz. 
Figure 3.1. 1-2 shows the gain, G, of the equivalent paraboloid as a function of normalized 
illumination misalignment distance, AR, )( for five edge illumination values. The gain loss 
caused by illumination misalignment can be separated into two loss components corresponding 
to decreased spillover efficiency and decreased illumination efficiency. As shown by Figure 3.1.1- 
3, most of the gain loss caused by illumination misalignment consists of spillover loss. 
Expressed in decibels, the decrease in spillover efficiency relative to the properly aligned case 
varies from about 13 to 36 dB at an normalized illumination misalignment distance of 2. The 
remaining gain loss, 3 to 9 dB at an normalized illumination misalignment distance of 2, is 
comprised of decreased illumination amplitude efficiency as shown in Figure 3. 1.1-4. 
Illumination phase efficiency is 100% for this case since the system is focused. The overall 
aperture efficiency, Cjkp , is shown in Figure 3. 1.1-5 as a function of normalized illumination 
misalignment distance, AR,,; this is found from the data in Figure 3.1. 1-3 using (3.1. 1-1). 

Figure .3. 1.1-6 shows the sidelobe level, SLL, as a function of normalized illumination 
misalignment distance, AR,,, for the five edge illumination values. Although the plots are not 
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Figure 3. 1.1-3. 


Test case equivalent paraboloid spillover efficiency (c ) 
normalized illumination misalignment distance (AR^). 
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Illumination amplitude efficiency (eamp), % 



Normalized illumination misalignment distance (ARn) 


Figure 3. 1.1-4. Test case equivalent paraboloid illumination amplitude efficiency (c Mlip ) as a 
function of normalized illumination misalignment distance (AR eqI1 ). 
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Aperture efficiency (tap), % 



Normalized illumination misalignment distance (ARn) 


Figure 3.1. 1-5. Test case equivalent paraboloid aperture efficiency (e ) as a function of 
normalized illumination misalignment distance (AR ). 
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Sidelobe level (SLL), dB 



Normalized illumination misalignment distance (ARn) 


Figure 3. 1.1. -6. Test case equivalent paraboloid sidelobe level (SLL) as a function of normalized 

illumination misalignment distance (AR cqil ). 
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smooth, the relative sidelobe level increases monotonically with increasing illumination 
misalignment. The roughness in the curves is caused by the consecutive blending of the 

increasing near-in sidelobes with the main lobe. Figure 3. 1.1-7 shows the peak cross-polarization 
level, 


XPOL = <UWU - (3.1. 1-6) 

as a function of normalized illumination misalignment distance, AR,,, for the five edge 

illumination values. Beam efficiency, defined as 
25 ( fl 3dB) 2 * 

J J m*) sin0 dO d<f> 


BE = 




(3. 1.1-7) 


where ^ is angle from the main beam gain peak to the the half-power point of the pattern, 
?{9, <t>) is the power pattern of the antenna, and P t is the power transmitted by the feed, is 
shown in Figure 3. 1.1-8 for several edge illumination values. In order to achieve the required 
90% beam efficiency, the edge illumination must be lower than -12 dB for a system with no 
illumination misalignment and still lower if illumination misalignment exists. 


3.1.2 Phase Error Effects 


The second source of antenna performance degradation during beam scanning is 
aperture phase error resulting from distortion of the aperture constant phase surface. The effects 
of aperture phase errors, or aberrations, on antenna performance have been extensively studied 
by Born and Wolf m connection with optical lens design [3]. Born and Wolf [1] represent 
aperture phase errors with cylindrical Zernike polynomials of p and * to define a constant phase 
surface above the aperture. For small aberrations, the first five terms of the aberration series 

are sufficient to describe the constant phase surface [1]. These terms, called the Siedel 
aberrations, are 
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Figure 3.1. 1-7. Test case equivalent paraboloid cross-polarization level (XPOL) as a function of 
normalized illumination misalignment distance (AR eqn ). 


Beam Scanning in Cm cgrain Reflector Antenna Systems 


46 



Beam efficiency (BE), % 


N X 

\ N 


\ \ 
X \ ' 
X \ 


-5 dB edge Hum. 
*1 0 dB edge Mum. 
"1 5 dB edge Mum. 
-20 dB edge Mum. 
•25 dB edge Mum. 


>V\ 

A X X 

\\\ \ 
•>\ \ X 
' \ X 
X\X 


0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Normalized illumination misalignment distance (AR„) 


Figure 3 . 1 . 1 - 8 . Test case equivalent paraboloid beam efficiency 1 BE 1 « f *• 

normalized illumination ° f 

S«»nnm« in CWmb Rnflector Aat «u Syrtnm. 



AL = apcosij) + Pp 2 + HP^cos'lO + S <p 3 cos0 + sp 4 ( > J - 1 

„d known, respectively, ns distortion, curvature of field, astigmatism, coma, and spherical 
aberration [1]. The constant aperture phase surface produced by each these aberration terms is 
shown in Figure 3. 1.2-1 [1]. The effects of these aberrations on antenna patterns are usually 
obtained by the simplification of the Zernike polynomial to a one-dimensional polynomial by 
setting 4 equal to either 0 or f [1]. With this simplification, the Siedel aberration series 


becomes 

AL(p) = ap + (/? + l)p 2 + 6 P 3 + £ P 4 


( 3 . 1 . 2 - 2 ) 


in the <p = 0 plane. 


The linear phase error term, op, doe, not degrade the antenna pattern but rather steers 
the beam away from the unscanned boresight direction [1], The amount of beam scanning 

caused by this aberration is 




. a( «A \ (*•»•«) 

= s,n (jvO 

where D M is the diameter of the antenna aperture and A is the wavelength at the frequency o 
operation [1]. This is the desired constant aperture phase surface which should be attained 
during beam scanning while minimizing the other aberration terms and illumination errors. A 


secondary pattern effect which is associated with this aberration term is a small change in the 
size of the projected aperture. For a symmetric aperture which is perpendicular to the 
unscanned beam, this change is a reduction of the projected aperture diameter which is 

proportional to cos0 o . 


The quadratic phase error term, (/? + 7 )p\ has no beam scanning effect since it 
represents a symmetric aperture phase taper. This aberration will, however, lead to an increase 
in the overall sidelobe level and pattern null-filling [1]. The exact effects of quadratic phase 
error can be determined by an analysis presented by Silver [4], This analysis gives the far-field 
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(3. 1.2-4) 


radiation pattern of the aperture from 
-1 

where u = ” ntf and f W is assumed to ** an even function whicl * re P resents the 

amplitude distribution across the aperture [1]. For small quadratic phase aberrations, the far- 

field power pattern is 

POO = 5*L{g>> + 0 ! [&,>)] 2 ) (3'-'-2-5> 

where gju) is the far-field radiation pattern with no aberration and g 0 » is the second 
derivative of the aberration free far-field radiation pattern [1]. The effects of the quadratic 
phase error term on the pattern of an aperture with a uniform amplitude are shown in Figure 
3. 1.2-2 for 0 = 0, and -y [1]. As 0 increases, the pattern continues to deteriorate with the 
main lobe eventually bifurcating but symmetry is maintained about the 0 = 0 axis [1]. 

The cubic phase error term, 6p 3 , introduces additional beam scanning since it represents 
an asymmetric aperture phase taper and also raises the peak sidelobe level [1]. Using Silver’s 
method to find the effects of the cubic phase error term on the far-field power pattern results in 
P(u) ~ ^ L {g 0 (u) + 6g 0 ’(u)} 2 (3. 1.2-6) 

where g 0 (u) is the far-field radiation pattern with no aberration and g 0 "(u) is the third 
derivative of the aberration free far-field radiation pattern [1]. The effects of the cubic phase 
error term on the pattern of an aperture with a uniform amplitude are shown in Figure 3. 1.2-3 
for 7 = 0, and y [1]. The cubic phase error term causes the main beam of the aperture to be 
steered farther away from the unscanned boresight direction, 0 = 0, as 7 increases. Also, the 
sidelobes increase on the side of the main beam away from 6 = 0 while the sidelobes on the side 
of the main beam toward 9 = 0 decrease. 

Like the quadratic phase error term, the quartic phase error term, cp 4 y has no beam 
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Figure 3. 1.2-2. Effects of quadratic phase error on an antenna pattern [1], 
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Figure 3. 1.2-3. Effects of cubic phase error on an antenna pattern [1]. 
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scanning effect since it represents a symmetric aperture phase taper. Also, like the quadratic 
phase error term, the quartic aberration will increase the overall sidelobe level and fill in the 
pattern nulls although the pattern degrades less than for a quadratic phase error of the same 
peak value [Ij. Using Silver’s analysis, as for the quadratic and cubic phase error terms, the far- 
field power pattern of the quartic phase error term is 

PW=^{ 8o ! W + £ V ,, (“)] 2 } (3. 1.2-7) 

where g c (u) is the far-field radiation pattern with no aberration and g 0 (4) (u) is the fourth 
derivative of the aberration free far-field radiation pattern [1], The effects of the quartic phase 
error term on the pattern of an aperture with a uniform amplitude are shown in Figure 3. 1.2-4 
for t - 0, f, and £ [1], As expected from the similarities between (3. 1.2-7) and (3. 1.2-5), Figure 

3.1.2- 4 is much the same as Figure 3.1.2-2 except with a smaller pattern degradation for a given 
peak aberration [1]. 

The boresight gain loss for quadratic, cubic, and quartic phase errors is shown in Figure 

3.1.2- 5 [1]. With the exception of the cubic phase error, the effect of an aberration with a given 
peak phase error decreases with increasing order. The effect of cubic phase error is greater 
because the beam peak is scanned away from boresight by the resulting asymmetric phase 
distribution. This relationship between aberration order and pattern degradation for a given 
peak phase error can be explained by noting that for a given peak phase error, the total error 
decreases with increasing aberration order as shown in Figure 3. 1.2-6. In addition to gain loss, 
as shown for each of the non-linear phase error terms, aberrations of quadratic or higher order 
also cause an increase in peak sidelobe level and can cause the filling of pattern nulls [1], Of the 

Siedei aberrations, cubic phase errors exhibit the greatest pattern degradation for a given peak 
phase error [1], 
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Figure 3.1.2-5. Boresight antenna gain loss for the Siedel aberration terms [1]. 
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Figure 3. 1.2-6. Phase error as a function of distance from aperture center for a one dimensional 

aperture. 
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3.2 Beam Scanning Effects of Equivalent Paraboloid Offset 


As discussed in Sections 2.4, a Cassegrain reflector system can be represented with 
reasonable accuracy by an equivalent paraboloid, which may or may not be offset, for offset 
Cassegrain systems. Because this equivalence is developed using geometrical optics it is only 
valid for focused systems. For a scanning offset Cassegrain reflector antenna system this means 
that unscanned performance, especially cross-polarization level, can be improved by using the 
Rusch condition discussed in Section 2.5. The additional concern for scanning offset Cassegrain 
reflector antenna systems is that scan performance might be degraded by the use of a reflector 
system geometry which has an axi-symmetric equivalent paraboloid. 

Parameters of offset Cassegrain reflector antenna systems with three equivalent 
paraboloid offsets are summarized in Table 3.2-1. These systems were studied to determine the 
effect of equivalent paraboloid offset on the scanning characteristics of offset Cassegrain reflector 
antennas. The three systems were synthesized using the main reflector of the test case which 
will be discussed in Chapter 4. The feed point was allowed to vary to create systems with either 
an axi-symmetric equivalent paraboloid or an equivalent paraboloid which was just fully offset 
above or below the axis of symmetry. The scanning performance of these systems was 
determined by using GRASP7 analysis with geometrical optics/geometrical theory of diffraction 
at the subreflector and physical optics/surface integration at the main reflector. The geometries 
of these three systems with their equivalent paraboloids are shown in Figure 3.2-1. 

3.2.1 Feed Motion Scanning 

The study of the effect of equivalent paraboloid offset on offset Cassegrain reflector 
antenna scanning characteristics was first performed using feed motion to scan the systems. 
During this trial, the feed of each system was translated in three dimensions to achieve the 
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Table 3.2.- 1. Characteristics of three offset Cassegrain reflector antenna system using the same 
main reflector but with three different equivalent paraboloids. 

Main reflector diameter (D M ) 

10.63 meters 

Main reflector offset height (H 0 ) 

7.795 meters 

Main reflector focal length (F M ) 

13.5 meters 

Subreflector interfocal distance (2c) 

6.9 meters 

Subreflector eccentricity (e) 

1.919 

Equivalent paraboloid focal length (F eq ) 

42.5 meters 

Equivalent paraboloid offset height (H^) 

0 and ± 5 meters 

Feed point (Xf, y f , z f ) 

near (-0.442,0.0,6.614) 

Feed pattern 

-15 dB subreflector 
edge illumination 
Gaussian pattern 

Frequency of analysis 

20 GHz 

i 
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Figure 3.2-1. Offset Cassegrain reflector systems and their equivalent paraboloids. 
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lowest possible aperture phase error as defined by the transmit mode raytracing error functional 
which will be described in Section 5.2. The transmit mode raytracing error functional was used 
to position the feed to allow the same error functional to be used for both feed and subreflector 
motion scanning. The systems were fed by a feed pattern which produced a -15 dB edge 
illumination for the analysis. The feed pattern was held fixed at the optimal unscanned 
pointing for each system to simplify the synthesis and analysis. 

The gain, G, of the three systems is shown in Figure 3.2. 1- 1. As for the unscanned case, 
the gain of a scanned offset Cassegrain reflector antenna is essentially unchanged by equivalent 
paraboloid offset. Spillover efficiency, c sp , shown in Figure 3.2. 1-2, is also relatively unchanged 
by equi valent paraboloid offset. The sidelobe level of the three systems could not be calculated 
because phase errors caused null filling at relatively low scan angles. The only significant effect 
of equivalent paraboloid offset on the scan performance of offset Cassegrain reflector antenna 
systems is in cross- polarization level. The cross- polarization levels, XPOL, of the three systems 
are shown in Figure 3.2. 1-3. The system with the axi-symmetric equivalent paraboloid both 
reaches the lowest cross-polarization level and generally performs the best although the systems 
with offset equivalent paraboloids each exhibit a lower cross-polarized component over a small 
portion of the scan region. This effect is caused by the change in illumination due to the 
required scanning motion. 

3.2.2 Subreflector Motion Scanning 

The study of the effect of equivalent paraboloid offset distance on the scanning 
characteristics of offset Cassegrain reflector antennas was also performed using subreflector 
translation and rotation to scan the three systems of Table 3.2-1. The subreflector position was 
determined to minimize the aperture phase error as defined by the transmit mode raytracing 
error functional which will be discussed in Section 5.2 as in the feed motion scanning study. 
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Figure 3*2. 1-2. Feed scanned spillover efficiency (e sp ) as a function of scan angle in the plane of 
offset (0 O ) for three offset Cassegrain systems with three different equivalent paraboloids. 
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Also, the feed illumination was again geometrically defined to produce a -15 dB edge 
illumination and was fixed to the optimum unscanned feed pattern pointing. 

Figure 3.2.2- 1 shows the gain, G, of the three systems as a function of scan angle in the 
plane of ofTset, 0 O . Again, the gain variation caused by equivalent paraboloid offset is small 
although the difference in scanned gain is greater than for the feed scanned case. The spillover 
efficiency, e jp , of the three systems is shown in Figure 3. 2.2-2 and, like gain, is relatively 
insensitive to equivalent paraboloid offset. Sidelobe level, SLL, could be found for the 
subreflector scanned case and is also nearly unaffected by equivalent paraboloid offset as shown 
in Figure 3.2. 2-3. As for the feed scanned case, cross-polarization level, XPOL, is the only 
pattern characteristic which was found to be strongly affected by equivalent paraboloid offset. 
Figure 3.2. 2-4 shows that the cross- polarized component is minimized in the system with an axi- 
symmetric equivalent paraboloid but can be lower in certain regions of the scan range for either 
system which has an offset equivalent paraboloid. 

3.3 Beam Scanning the Offset Cassegrain Antenna by Feed Displacement 

The conventional method of mechanical beam scanning with the offset Cassegrain 
reflector antenna system is lateral feed displacement [5]. Lateral feed displacement is in 
common use to form multiple spot beams using a single prime focus paraboloidal reflector 
antenna for satellite-to-gronnd applications [6-8]. Increasing geostationary satellite antenna 
performance requirements prompted INTELSAT to commission a study by Krichevsky and 
DiFonzo at COMSAT of beam scanning offset Cassegrain reflector antennas by optimal feed 
displacement [9]. This study developed a method to determine the optimal feed location for any 
desired main beam direction. Krichevsky and DiFonzo give results for beam scanning in the 
same offset Cassegrain reflector antenna system using lateral feed displacement so that the 
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Figure 3. 2. 2-2. Subreflector scanned spillover efficiency (f 8p ) as a function of scan angle in the 
plane of offset (0 O ) for three offset Cassegrain systems with three different equivalent 

paraboloids. 
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Figure 3.2.2-4. Subreflector scanned cross- polarization level (XPOL) as a function of scan angle 
in the plane of ofTset (0 O ) for three offset Cassegrain systems with three different equivalent 

paraboloids. 
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increased performance of their optimum feed position scanning method can be shown [5]. 

The Krichevsky and DiFonzo study of optimal beam scanning used the reflector system 
geometry shown in Figure 3.3-1 [5]. This system consists of a paraboloidal main reflector with a 
focal point at F,, a hyperboloidal subreflector with focal points at F, and F 2 , and a point source 
feed located at F 2 when the system is unscanned [5]. The coordinate system is constructed such 
that F 2 lies at the origin and the unscanned aperture plane is perpendicular to the z-axis [5], 
The subreflector rim is defined by the intersection of a cone with its vertex at F 2 and the 
hyperboloidal surface on which the subreflector lies [5]. The angle between the axis of the cone 
which defines the subreflector edge and the +z-axis is 0, and the half-angle of the cone is 0 2 [5], 
Main reflector size is determined by the area illuminated by the ray bundles which are reflected 
from the subreflector for all beam directions after the optimal feed positions are found [5], This 
method of sizing the main reflector is much less stringent because it assures that no spillover will 
occur at the main reflector. Krichevsky and DiFonzo defined the plane of lateral feed 
displacement as the plane which contains F 2 and is perpendicular to the axis of the cone which 
defines the subreflector edge [5]. 

Krichevsky and DiFonzo determined the optimum feed position for a given beam 
direction by finding the intersection of the locus of feed positions which produce the desired 
beam direction with the locus of feed positions which yield the minimum aperture phase error 
[5]. For simplicity, in the analysis presented here, that of Krichevsky and DiFonzo [5], the 
antenna system is assumed to be cylindrical and only rays which lie in the xy-plane are 
considered. The complete derivation of the optimum feed position is available in Krichevsky [9], 

The locus of feed positions which produce the desired beam direction is found by tracing 
rays transmitted from a point source feed at Q(0, y 0 , z 0 ) as shown in Figure 3.3-1 [5]. The total 
optical path length of a ray which satisfies the conditions imposed by geometrical optics at the 
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Figure 3.3-1. Antenna system geometry for the Optimum Feed Position Scanning of Krichevsky 

and DiFonzo [5]. 
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mam and subreflectors and strikes the subreflertor at Pj(0, y |t z,), the main reflector a( P 2 (0, 
y 2 , z 2 ). and the aperture plane at P 3 (0, y 3 , z 3 ) is given by 

(3.3-1) 


E L i 

i = l 


where 

= \A y i ' + ( Z i " z i-l) 2 (3.3.2) 

for i = 1, 2, or 3 [5]. The total optical path lengths and points of intersection with the aperture 
plane of the upper and lower rays in the system are given by L u and L, and (y u , zj and (y,, z,) 

as shown m Figure 3.3-2 [5]. From these parameters, the beam direction can be approximated 
by 

y u - y\ (3.3-3) 

for small scan angles and feed displacements [5j. 

Krichevsky and DiFonzo next present the series expansions of the total optical path 
length and transmitted ray/aperture plane intersection: 

l = l (0) J2 E c ,„,k y 0 k z o'"' k 


111=0 k=0 


and 


y = y (0) 


v m 


E E T m,k >0 k z 0 m ‘ k 

111=0 k=0 


(3.3-4) 


(3.3-5) 


where L (0) and y (0) are equal to L and y when the point source feed is at the unscanned feed 
point F 2 [5]. The value of v in the series expansions is the order of the desired approximation: 

a) first order approximation: 1 

H 

b) second order approximation: 1 

where R is the displacement of the point source feed from the unscanned feed point F 2 and F H is 
half the interfocal length of the parent hyperboloid of the subreflector [5]. By substituting (3.3- 
4) and (3.3-5) into (3.3-3), Krichevsky and DiFonzo find the series expansion for the locus of 
feed positions which produce a constant beam direction: 
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Figure 3.3-2. Krichevsky and DiFonzo [5] definition of beam scan angle. 
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°= E E A m.k >O k Z 0 "‘- k 

ni=0 k=0 


(3.3-6) 


where A m,k 18 a •»*“ of coefficients which are defined in terms of the geometrical parameters of 
the offset Cassegrain reflector antenna system [5]. The coefficients, A ln>k , and their derivation 
are presented in Krichevsky and DiFonzo [5] and Krichevsky [9]. The constant beam direction 
feed locus can be found from (3.3-6) by solving for y 0 in terms of z 0 , a, and the coefficients, 


Y m,k* 


Ftgure 3.3-3 shows several constant beam direction feed position loci for the offset 
Cassegrain reflector antenna system summarized by Table 3.3-1 [5]. The first-order, linear 
approximation to (3.3-6) for each constant beam direction feed locus is shown in Figure 3.3-3 by 
a dashed line which is labeled with the angular beam displacement caused [5], The second-order 
approximation to (3.3-6) is shown by a solid line which crosses the corresponding first-order 

approximation [5], Figure 3.3-4 shows the beam pointing error as a function of the beam 

direction predicted by the constant beam direction loci shown in Figure 3.3-3 (5). The beam 

directions used to determine the accuracy of the constant beam direction feed loci were 

calculated using a computer based pattern analysis code [5], The second-order approximation is, 
as expected, more accurate with a maximum beam pointing direction error of less than 0.05* 
compared to a maximum error of about l.T for the first-order approximation [5]. 

Krichevsky and DiFonzo [5] calculated the optimum feed position locus by dividing the 
subreflector and tracing n rays from the unscanned feed point to the subreflector with the same 
angular separation between the rays. This, in effect, divides the antenna for which the optimum 
feed position locus is being determined into n-1 small antennas [5). The substitutions 
*2,11 ~ 7 T 

(3.3-7) 

and 

* 

1,111,11 — IT ’ 1) + 0, - O.y 

' 1 7 7 Q\ 
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Figure 3.3-3. First and second order approximations to the Constant Beam Direction Locus of 
Krichevsky and DiFonzo [5] for the system of Table 3.3-1. 
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Table 3.3-1. Characteristics of the Krichevsky and DiFonzo [5] offset Cassegrain reflector 
antenna system used for Optimum Feed Position Scanning. 


Main reflector focal length (Fj) 
Subreflector focal length (F 2 ) 
Subreflector interfocal distance (2c) 
Feed angle subtended by subreflector 


100 A 
40A 
20A 
34.38° 
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SECOND ORDER 


t n= ion 



Figure 3.3-4. Beam pointing error given by physical optics analysis for the Constant Beam 
Direction locus of the system of Table 3.3-1 [5]. 
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where m is the ray number from 1 to allow the use of the constant beam direction solution 
method to find the beam direction for each of the subdivisions [5], The difference between the 
scan beam directions of the subdivisions can be minimized by using the constraint 

G(y °’ 2o) = >-~ £ K,,„ - «) 2 j (3.3-9) 

where a is the overall beam direction since the beam direction for each of the subdivisions is 

a n., n [5]- Since the solution must also lie on one of the constant beam direction loci, the 
complete error function is 


1(y„, 2 0 , /.) = G(y 0 , ,„) + . f v (, 0 , A, ll k )J (3.3-10) 

where „ is a Lagr.ngian multiplier [5], After trigonometric a„d .Igebreir simplification of 
(3.3-9), Krichevsky and DiFonzo found rhe following closed form solution for the optimum feed 


position: 

v m 


E E ,k V z 0 ‘"- k = o 

m= 1 k=0 


(3.3-11) 


where l m , k is a series of coefficients which are defined in terms of the geometrical parameters of 
the offset Cassegrain rellecto, antenna system [5J. The coefficients, I,„ k , and their derivation 
are presented in Krichevsky and DiFonzo [5] and Krichevsky [»]. As f„, the constant beam 
direction feed loci, v corresponds to the order of the desired approximation; either 1 fo, a first- 
order, linear approximation o, 2 for a second-orde, approximation (5J. 

The locus of optimum feed position, and the second-order approximation to the 
constant beam direction loci fo, the system of Table 3.3-1 are shown in Figure 3.3-5 (5). The 
dotted line in Figure 3.3-5 represents the possible feed positions for lateral feed displacement 
while the dashed and solid lines represent the first- and second-orde, approximations to the 
optimum feed position locus [5). The optimum feed position for each beam direction can be 
found at the intersection of the second-orde, approximation of the appropriate constant beam 
direction loci and the second-orde, approximation of the optimum feed position locus (5], The 
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required motion for optimum feed positioning beam scanning requires both a lateral and a 
longitudinal feed displacement. 

Krichevsky and DiFonzo [5] used physical optics to find the far-field patterns of the 
optimum feed position movement scanned system of Table 3.3-1. The reflector system was 
analyzed for a Potter horn feed with a radius of 2.3A at the optimum feed position for each 
beam direction [5]. The feed was repointed at each optimum feed position to align the axis of 
the feed with the angular center of the subreflector to minimize spillover [5j. Also, as mentioned 
above, the main reflector surface size was determined by allowing no spillover at the main 
reflector [5]. Figure 3.3-6 shows several scanned beams for the offset Cassegrain reflector 
antenna system with the feed for each beam located at the optimum feed position [5], Figure 
3.3-7 shows several scanned beams from the offset Cassegrain reflector antenna system fed by a 
3.1 A Potter horn which was displaced laterally [5]. In both cases, the pattern of the offset 
Cassegrain reflector antenna system tends to deteriorate less rapidly during downward beam 
scanning, but the optimum feed position scanned system shows less overall pattern deterioration. 
Note that the aperture is partly blocked by the subreflector for beams which are scanned more 
than about 5.8 below the unscanned boresight [5]. 


3.4 Beam Scanning the Offset Cassegrain Reflector Antenna by Subrcflcctor Tilt - The 
Foldes Type 6 System 

A reflector antenna system was proposed by Peter Foldes for use in the NASA Mission 
to Planet Earth radiometer design. This system, referred to as the Type 6 system, is 
summarized by Figure 3.4-1 and Table 3.4-1. The Type 6 system was designed to scan by 
subreflector tilt. For simplicity, the scanned systems were synthesized by tilting the subreflector 
to an angle and then finding the resulting beam scan from the physical optics analysis result. 
This also assured that the performance of the Type 6 system was maximized. Unfortunately, 
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Figure 3.3-6. Scanned beams of the system of Table 3.3-1 produced by the second order 
Optimum Feed Position scanning of Krichevsky and DiFonzo [5j. 


Beam Scanning in C a— r grain Reflector Antenna Systems 


80 






X-axis, m 





Table 3.4-1. Characteristics of the Type 6 reflector antenna system of Foldes [JO]. 


Main reflector focal length (F M ) 
Main reflector diameter (D M ) 

Main reflector offset height (H c ) 
Subreflector eccentricity (e) 
Subreflector interfocal distance (2c) 
Subreflector axis tilt angle (/?) 
Subreflector diameter (D s ) 

/ n 2 

Areal efficiency, ( M 

W + V 

Feed tilt angle (a) 

Intended scan range 
Frequency of analysis 
Feed pattern used for analysis 


35 meters 
25 meters 

17.5 meters 
2.81 

14.5 meters 
0 ° 

3.56 meters 

98% 

ir 

± 0.5* 

10, 20, and 40 GHz 

“15 dB Gaussian edge 
taper 


ScJ “ min * “ Reflector Antenno S yrt 
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this method for determining beam scanning parameters is not applicable to more sophisticated 
system because the multiple degrees of freedom would require excessive time to analyze as 
briefly discussed in Chapter 5. 

The feed pattern was not specified by Peter Foldes and so was chosen to produce a -15 
dB edge taper based on the beam efficiency results shown in Section 3.1. The scan performance 
results shown here were calculated using the T1CRA GRASP7 reflector antenna code. The 
system was analyzed using geometrical optics/geometrical theory of diffraction at the 
subreflector and physical optics surface integration at the main reflector. For purposes of later 
comparison, the Foldes Type 6 system was analyzed at 10, 20, and 40 GHz. 

Figure 3.4-2 shows the gain, G, of the Foldes Type 6 reflector antenna system as a 
function of frequency and scan angle in the plane of offset, 0 Q . The system displays the expected 
6 dB increase in gain for each octave increase in operating frequency. The significant increase in 
scan induced gain loss with increasing frequency indicates that the main source of error in the 
Foldes Type 6 reflector antenna system is aperture phase error. This conclusion is supported by 
the relative flatness of the curves for spillover efficiency, e sp , shown in Figure 3.4-3. As expected 
the spillover efficiency is essentially frequency independent with the slight variation being due to 
the increase electrical size of the system causing less diffraction effects. The illumination 
amplitude efficiency, £ ainp , shown in Figure 3.4-4 as a function of scan angle ill the plane of 
offset, 0 O , was estimated for the Foldes Type 6 system by using the results presented in Section 
3.1.1 because the presence of both amplitude and phase taper in the aperture field of a scanned 
reflector antenna prevent determining the amplitude efficiency factor independently. The 
illumination amplitude efficiency is essentially frequency independent and its small range 
supports the conclusion that the main source of scan loss in the Type 6 system is phase error. 
The illumination phase efficiency, of the Foldes Type 6 reflector antenna system is shown in 
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Figure 3.4-3. Spillover efficiency (e sp ) as a function of scan angle in the plane of offset (0J for 
the Type 6 reflector antenna system of Foldes. 
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Illumination amplitude 



Scan angle in the plane of offset (9 0 ), deg. 


Figure 3.4-4. Illumination amplitude efficiency (e„, p ) as a function of scan angle in the pli 
of offset ( 6 0 ) for the Type 6 reflector antenna system of Foldes. 


Beam Scanning in Cassegrain Reflector Antenna Systems 




Figure 3.4-5 as a function of frequency and scan angle in the plane of offset, 0 O . Although the 
illumination phase efficiency of the Type 6 system is fairly high at 10 GHz, at 20 and 40 GHz 
phase errors are clearly the dominant scan loss mechanism. Like the illumination amplitude 
efficiency, the illumination phase efficiency should be considered a diagnostic tool rather than an 
exact measurement because the illumination amplitude efficiency estimate was used in the 
calculation of the illumination phase efficiency. The overall aperture efficiency, £ ap , of the Type 
6 system is shown in Figure 3.4-6. 

Figure 3.4-7 shows the sidelobe level, SLL, for the Foldes Type 0 reflector antenna 
system as a function of frequency and scan angle in the plane of offset, 0 o . The sidelobe level 
should increase monotonically but blending of the increasing sidelobes with the main lobe and 
null filling cause the observed roughness in the curve. Figure 3.4-8 shows the cross-polarization 
level, XPOL, for the Type 6 system. The generally low cross-polarization level is caused by the 
long effective focal length, 101.5 meters, of the system. The offset of the equivalent paraboloid 
can be observed in the monotonic decrease in relative cross-polarization level with increasing 
negative beam scan. The beam efficiency, BE, of the Foldes Type 6 reflector antenna system is 
shown in Figure 3.4-9. As with the illumination phase efficiency, the beam efficiency is a strong 
function of frequency. Since the limit of the scan range is 90% beam efficiency, the scan range 
of the Type 6 system is about 1.5 1 ’ at 10 GHz, 0.8 at 20 GHz, and 0.6 at 40 GHz. 
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Humiliation phase efficiency (cj as a function of scan angle i 
( 0 ) for the Type 6 reflector antenna system of Foldes. 
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Figure 3.4-7. SideU.be hud (SLL) U. . fun.Uon of ecu,, „ lg ,e i„ , h e „ o(li=l (<J ^ 

Aype 6 reflector antenna system of Foldes. ' 
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Figure 3.4-8. Cross-polarization level (XPOL) as a function of scan angle in the plane of ofTset 
(0 o ) for the Type 6 reflector antenna system of Foldes. 
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Figure 3.4-9. Beam efficiency (BE) as function often,, angle j„ „l,„e of offset « ) f„, 

I ype 6 reflector antenna system of Foldes. ° 
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Chapter 4 

TEST CASE GEOMETRY SELECTION 


The selection of a test case geometry for this study is motivated by the antenna 
performance characteristics required by the Mission to Planet Earth. A NASA microwaves 
radiometric earth observation science steering panel met twice in 1990 to determine the mission 
parameters which would allow the project to make a significant contribution to the knowledge of 


climatic and meteorological phenomena [1], The requirements recommended by the panel are 
summarized in Table 1.2-1. The frequency and beamwidth restrictions require an aperture 
efficiency of 70% since a Virginia Tech study of the issues presented by large space antenna 
structures found that an antenna of up to 25 meters diameter can be launched by a single 
Shuttle Transportation System (STS) mission [1, 2]. 


As mentioned in Chapter 1, the necessity of sharing a geostationary platform with other 
experiments requires that the antenna be capable of beam scanning by sub-optics motion to 
•void disturbing other experiments. The aperture efficiency required by the desired beamwidth 
and frequency of operation effectively eliminate, conventional spherical rellecto, antenna systems 
from consideration despite their degrad.tion-free scan performance. The Cas^gr.in reflector 
antenna was selected since it is the simplest and most compact reflector configuration that 
allows beam scanning by subreflector motion. The high beam efficiency requirement indicates 

that an offset antenna configuration would be desirable to eliminate aperture blockage by ,h, 
subreflector/feed assembly. 
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The test ease main redeem, has a diameter of D M =10.6J mere,., a fed length of 
F -,3.5 meters, and an offset distanee of H,, = 7.795 meters. These dimensions were ehosen to 
allow possible verification of the synthesis procedure using a NASA Langley Antenna and 
Microwave Research Branch test article [3]. The configuration and dimensions of the test art.cle 
are shown in Figure 4-1. This test article is a 42.52% scale model of the proposed 25 meter 
radiometric earth observation reflector antenna. 

Since the results given in Section 3.2 show that the scan characteristics of Cassegrain 
reflector antenna systems are insensitive to equivalent paraboloid offset, the test case geometry 
was chosen to have an axi-symmetric equivalent paraboloid. This geometry selection should 
minimize the overall cross-polarization experienced by the system across the scan range. The 
Rusch condition for this system can be found using (2.5-7) with a center angle, 9 e , of 32.207: 

32.207 2e sin 0 (4-1) 

tan 2 e 2 + 1 - 2 e cos 0 

The eccentricity of the subreflector hyperboloid was chosen to be 1.919 to approximate the 
relative subreflector size of the 25 meter diameter Foldes Type 6 reflector antenna system which 
was originally proposed. This results in a subreflector diameter of about 1.4 meters for the test 
case. The subreflector area used in the synthesis and the GRASP7 physical optics analysis is the 
illuminated portion of the parent hyperboloid during unscanned conditions as determined by 


geometrical optics raytracing. 

The tilt angle, 0 , of the line between the hyperboloid foci can now be found to be 0 = 
3.676" from (2.6-1). The interfocal length of the hyperboloid was chosen to be 6.900 meters to 
place the feed point near the the subreflector/feed boom. These selections lead to a feed point 
located at {0.442, 0.000, 6.614} meters. The feed tilt angle, «, which will align the feed 
boresight with the angular center of the subreflector and the axis of the equivalent paraboloid 

can now be found using (2.4-19): 
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Figure 4-1. NASA Langley AMRB test article configuration and dimensions [3]. 
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(4-2) 


un? = u»(f). 

Solving (4-2) with these selections of e and fi gives « = 11.640". The test case is fed with an x- 
polarized feed. The feed pattern is defined independent of frequency to produce a 15 dB edge 
taper on the subreflector. This edge taper was chosen based on the beam efficiency results of 
Section 2.1.1. The geometry of the test case is given by Figure 4-2 and Table 4-1. 

The unscanned far-field pattern of the test case was evaluated using the TICRA 
GRASP7 reflector analysis code. The evaluation was performed by geometrical theory of 
diffraction analysis at the subreflector and physical optics/surface integration at the main 
reflector. The gain pattern, G<u,v), of the unscanned test case system is shown in Figure 4.3 
over a rectangular uv-grid of ±0.57". The cross-polarization pattern, XPOL ( u,v), of the 
unscanned test case system is shown in Figure 4.4. The far-field pattern characteristics of the 
unscanned test case are summarized in Table 4.2. 
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Figure 4-2. Plane of offset view of the final test case configuration. 
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Table 4-1. Characteristics of the final test case geometry. 


Main reflector diameter (D M ) 

10.63 meters 

Main reflector focal length (F^j) 

13.5 meters 

Main reflector offset height (H c ) 

7.795 meters 

Subreflector axis tilt angle (/?) 

3.676° 

Subreflector eccentricity (e) 

1.919 

Subreflector interfocal length (2c) 

6.9 meters 

Subreflector diameter (Dg) 

1.4 meters 

Area! ,tr,c,e„cy 

98.3% 

Subreflector scanning motions 

x, y, and z translations 
a and ft rotations 

Feed point (f) 

{-0.442, 0.000, 6.614} 
meters from the main 


reflector vertex 

Feed tilt angle (a) 

11.64° 

Feed pattern 

Gaussian pattern with 
15 dB edge taper at 
7.157° 

Frequency of operation 

20, 40, and 80 GHz 
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Table 4-2. Far-field pattern characteristics of the unscanned 


test case system. 


Far-field pattern characteristic 

Frequency, GHz 

Units 

20 

40 

80 

Gain, G 

65.67 

71.72 

77.78 

dBi 

Sidelobe level, SLL 

-29.00 

-29.32 

-29.45 

dB 

Cross- polarization level, XPOL 

-67.33 

-73.40 

-74.70 

dB 

Beam efficiency, BE 

93.28 

93.44 

93.63 

% 

Aperture efficiency, e hp 

74.44 

74.95 

75.63 

% 

Spillover effiency, e 

sp 

94.74 

94.85 

94.98 

% 

Illumination amplitude effiency, e Amp 

78.57 

79.02 

79.63 

% 

Half-power beamwidth, HPBW 

0.092 

0.046 

0.023 

0 
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Chapter 5 

SIMPLIFIED ERROR FUNCTIONALS FOR GEOMETRICAL 
OPTICS SYNTHESIS 


As discuss*! previously, the res! ruse developed |„ Chapter 4 be, m scuus by subrellector 
While the primary goal of this effort is the maximization of beam efficiency during 
ben scanning, analysis „f the pattern of the antenna determine beam efficiency |, 
impractical for use as an error functional in determining the optimal position of the subrellector 
during beam scanning. The physical optics analysis required to determine the far-field pattern 
of a Cassegrain reflector antenna system of this size requires approximately 6 minutes using a 50 
MHz 80480 compute,. An error functional of this complexity would result in an optimization 
which would require over one month of compute, time fo, a single scan direction. I, is highly 

desirable then to create a simple, computationally efficient estimate of the antenna’, 
performance. 


The desire that the error functional be easy to evaluate requires that the error functional 
be defined as a geometrical optics process with a minimal number of mathematical operations. 
Two error functionals, the correcting subrcllec.o, surface fating error function., and the 
transmit mode r.y.racing error functional, were investigated. Both of these error functionals 
involve a minimum of mathematical evaluation, each require, „„|, . ,i„ g | e summation aud i. 
tolerant of a much coarser analysis grid than is physical optics analysis. Each error functional 
was tested by using it to determine the optimum subreflector position f„, several scanned beam 
directions. Each case was then analyzed using the geometrical theory of diffraction at the 
Simplified Error Function^, for Geometric^ Optic Synthmi. 
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subreflector and physical optics/surface integration a< the main reflector using the T1CRA 
GRASP7 reflector antenna analysis package. 

S.1 The Correcting Subrefloetor Surface Fitting Error Functional 

The correcting snbreflector surface fitting error functional del ermines the optimum 

position for the unscauned subreilector by fitting the repositioned unscanned snbreflecto, to a 

correcting subreilector for .hr desired scan direction. Bec.usc the optimization is hinted on 

biting the available subject., to a full, illuminated correcting subteflector. thte error 

functional should be expired to emphasize the minimization of illumination amplitude errors at 

the expense of illumination phase errors. A simple, version of this error functional was pro, rosed 

by Kitsuregawa [1] for beam scanning by snbreflector motion in telleetor antenna systems where 

the main reflector and feed cannot be moved. The fitting erro, used b, Kitsuregawa is given by 

N . (5.1-1) 

<k= £[( S 'i-Si)-nsi] • 

where S'; is the i — point on the repositioned unscanned snbreflector, S; is the i point 
correcting snbreflector for the desired direction of scan, and n Si is the unit normal of the 
correcting snbreflector at the i * point [1]. While Kitsuregawa [1) uses this error functional in a 
system of least squares equations, it can also be used as an error estimate for optimization. 

5.1.1 Analytic Development of the Correcting Subreflector 

The test case reflector antenna system discussed in Chapter 4 is described by grids of 
points which define the parent surface of the main and subrefieelors. These grid, are found by 
using geometrical optics r.ytracil.g in the lest case geometry. As discussed in Section 2.3, 
Snell's law is satisfied at the main tefieclot and the total path length of the rays from the 
aperture plane through the system to the feed point is held constant to find the subrefieeto, 
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points. This procedure, developed by YVernlz [2]. is also used in the correcting subreflector 

surface fitting error function to find the phase error correcting subreflector for the scanned 
system. 


Figure 5.1. 1-1 shows the geometry used to determine the subreflector points. The rays 
used to define the system form a 25 x25 grid in the aperture plane which circumscribes the 

projected aperture of the antenna system. A ray incident on the main reflector from the 
aperture plane has a unit vector 

Fj = {-sin0 o cos0 o x, -sin0 o sin<6 o y, cos0 o ?}, (5. 1.1-1) 

where 6 0 is the scan angle from the z-axis and <t> 0 is the scan angle from the x-axis toward the y- 
axis as shown in Figure 5. 1.1-2. 


The length of a ray from the aperture plane to the main reflector is 
|R, | = cos0 o (z A + tan(-0 o )cos^ o x M + tan(-0 o )sin* o y M + z M ), (5.1. 1-2) 

where z A Is the z-coordinate of the intersection of the aperture plane and the z-axis and 
{ X M’ ^M’ z m} * s a P°i n t on the main reflector. 

The unit vector of a ray reflected from the main reflector is found by satisfying Snell’s 


law at the main reflector. This gives 
r 2 = r, - 2(r, -n M )n M , 

where n is the unit normal at the main reflector point given by (2.1-11) 
"M = (\/ x M 2 + >M 2 + 4F 2 ) 1 {-x M x, -y M y, 2Fx}, 


(5.1. 1-3) 


(5.1. 1-4) 


where F is the focal length of the main reflector paraboloid. Five intermediate results are next 


calculated to simplify the equations for the subreflector point. These are 

Ij = r i» ~ 2 ( r i ‘ n M ) n Mx 
T \z ' 2 (^1 

! 2 = F i y - 2 ( ? i -n M )n My , 


(5. 1.1-5) 


(5. 1.1-6) 
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Figure 5. 1.1-1. Correcting subreflector synthesis geometry. 
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Figure 5. 1.1-2. Definition of scan angles. 
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(5. 1.1-7) 


i - i /1±1 l 
13 ■ l2 v i-i 2 2 ’ 

I 4 = L p - 1 Rj | , and 

i 5 = 2 [i 4 /TTT7 ^ + M x m* *f) + yr)] 

where L p is the desired path length through the system and 


(5.1. 1-8) 


| x f , y r , z f J is the feed point. The 


subreflector point, |x s , y s , z s }, can now be found by 

* 4 2 + Mm - (*f - x m ) 2 ~ (yf ~ y\i ) 2 + z m 2 - z f 2 
Z s - I 5 + ~ Z M 2z f 

x s = ( Z S * z m) + X M> 
y$ = l 3 (z s - * M ) + y M- 


(5.1.1-9a) 
(5.1. l-9b) 
(5.1.1 -9c) 


The length of the ray from the subreflector point to the feed point is 

|R 3 |= v/(x S - x f ) 2 + (ys - yf ) 2 + ( Z S - z f ) 2 (5.1.1-10) 


so the unit vector of the ray from the subreflector point to the feed point is 


r o — 


= {t 


iSt >f >s „ 'r‘s ; 

IN IN IN 


} 


The subreflector unit normal can now be found by taking the difference between 


(5.1.1-11) 
?2 and f 3 . 


This results in 

~ - ( f 3x ~ f 2x) x + ( z 3y ~ ^2y)y +( r 3» ~ r 2») Z 
\/( *" 3x ‘ ^2x ) 2 + ( ^3y * ^2y ) + ( ^3z * ^iz ) 


( 5 . 1 . 1 - 12 ) 


The set of subreflector points and normals generated by this procedure is either saved to a data 


file in the case of test case synthesis or used in the error functional as the correcting subreflector 
during optimization. 


5.1.2 Correcting Subreflector Surface Fitting Optimization Implementation - CSSFT1 

The correcting subreflector surface fitting error functional is a slight, extension of the 
error functional of Kitsuregawa [1] given by (5.1-1). For the optimization, the error function is 
expanded to a double summation over the grid of points of the subreflector but is otherwise 
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unchanged. This results in 

= ,5 [^ S iJ ' S ij)‘"Sij] 2 . (5. 1.2-1) 

where S' t - is the ij ^ point on the repositioned test case subreflertor, S ;j is the i j Hi point on 
the correcting subreflector found using the method on Section 5.1.1, and n SiJ is the normal of 
the correcting subreflector. This error functional is then used to estimate the pattern error 
caused by the current position of the test case subreflector. To achieve this, the error functional 
is used in the two-stage optimization process which has the function blocks shown in Figure 
5.1.2-1. The FORTRAN source code which implements this process is listed in Appendix 1. 

The inner loop is a five-dimensional Powell’s method optimization [3] which positions 
the test case subreflector to best fit the correcting subreflector. The test case subreflector can be 
translated in the x-, y-, and z-directions and rotated about its center grid point in o- and 0-tilts 

as shown in Figure 5. 1.2-2. The position of a translated and rotated test case subreflector point 
is given by 

X S' - x t + x §r + COSOr r( x s - x s r ) + Slnft r sin /*r(ys * ^ + S " m r cos 0 r ( z s ’ z s r ). (5. 1.2- 2a) 
y S' = y ‘ + y Sr + C0S ^r(ys - y 5r ) - Sin0 r (z- - Zg r ), and (5.1.2-2b) 

z$< - z t + z s r + sin« r (x£ - Xg r ) + cosa r sin0 r (y- - y-^) + cas« r cos0 r (z^ - Zg r ) (5.1.2-2c) 

where { x s „ y $ „ z s ,} is the repositioned test case subreflector point, { x g , y § , z-j is the test case 
subreflector point, { Xg r , y^, Zg r j is the point about which the rotations occur, and { x t , y t , z t J 
vector of translations. The outer optimization loop is a Golden Section optimization 
process [4] which varies the path length through the antenna system to allow the calculation of a 
family of correcting subreflectors. This step is necessary because the curvature of the correcting 
subreflector is greater for higher total path lengths. 
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5.1.3 Predicted Subreflector Motions 

The optimum subredeeto, position, were dele, mined using the meeting subredeeto, 

surf.ce dtting error function.! for three se.n paths: . sc.n from 9„ = *° *° " 10 f " 

K = O'. . sen from = O' to *. = 180' for »„ = 0.5', .nd . se.n from = O’ to *„ = 180' for 
« o = 1.0*. These sc.n path. represent 30 individu.l sen directions will, . .»er.ge 
computational time of nearly 45 minutes per direction or a total time of jus, over 29 hours. An 
earlier version of the optimisation which u«d dv, linear optimisation routines instead of the 
Powell’s method inner optimisation loop was noted to be substantially slower. 

The constant 4 scan path from 9„ = -1.0* to «„ = 10* was approximated by 21 discrete 

scan directions at 0.1* intervals. The subredeeto, motion for this scan path consists of only X- 

and s- translations and o-ro.a.ion since the main beam is being scanned in the plane of offset. 

The X-, y-, and s-transl.tions required fo, this scan path me shown in Figure 5.1.31 as a 

function of scan angle in the plane of offset. The maximum translation fo, this scan path is 

about 0.25 meter at either limit of scan. The o- and flotations required fo, this sc.n path are 

shown in Figure 5.1.3-2 as a function of scan angle in the plan, of offset. The maximum 

rotation required for this scan path is less than 1 . 

/i a( d - n V and 0 =1.0* were each approximated by 11 

The constant 6 scan paths at 0 O - U.J> ana v Q - i.v 

r ± in A i — IKff at 18* intervals. The x- t y-» and z-translations 

discrete scan directions from = 0 to <t> 0 - l«U at 10 interval 

required for these scan paths are shown in Figure 5.1.3-3. The required total translations for 
these scan paths are nearly constant at approximately 0.11 meter for 0 o = 0.5" and 0.25 meter 
for 0 o = 1.0". The a- and ^-rotations required for these scan paths are shown in Figure 5. 1.3-4. 
As for the case of scanning in the plane of offset, the maximum required rotation is slightly less 

than T. 
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Figure 5. 1.3-2. Optimum subreflector «- and ^-rotations (« r , 0 r 
determined using the correcting subreflector error functional as a 

plane of offset (0 O ). 
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5.1.4 Computational Results using TICRA CRASP7 Physical Optics Analysis 

The performance of the correcting subreflector surface fitting error functional was 
evaluated by performing physical optics analysis to find the far-field radiation pattern 
characteristics of the scanned test case geometries at 20, 40, and 80 GHz. These frequencies 
were chosen to allow comparison with the results for the Foldes Type 6 reflector antenna system 
discussed in Section 3.4. Also, the two octave frequency range assists in the determination of 
the reasons for pattern degradation during beam scanning. 

The selection of geometrical theory of diffraction analysis at the subreflector and 
physical optics/surface integration analysis at the main reflector was mandated by the excessive 
time required for physical optics/physical optics analysis. At 20 GHz a dual physical optics 
analysis requires approximately 2 hours to complete on a 50 MHz 80486 computer. 
Furthermore, the time required for a PO/PO analysis increases roughly as the cube of the 
frequency. A further reduction in computing time was realized by limiting the output far-field 
pattern to an area on the uv-plane which contained a constant number of sidelobes with respect 
to frequency. This reduction in the output pattern area allowed the number of integration 
points in the main reflector surface integration to remain constant rather than quadrupling with 
each octave increase in frequency. 

Figure 5. 1.4-1 shows the gain, G, of the test case reflector antenna system as a function 
of scan angle in the plane of offset, 6 Q . As expected, the unscanned gain increases 6 dB for each 
octave increase in frequency indicating no loss of efficiency with increasing frequency. However, 
the scan loss varies greatly with frequency indicating phase error across the aperture. The 
sidelobe level, SLL, of the test rase is shown in Figure 5. 1.4-2 as a function of scan angle in the 
plane of offset, 6 q . The sudden changes in sidelobe level at the higher frequencies are caused by 
sidelobes blending with the main lobe and indicate that significant phase errors are occurring. 
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Figure 5. 1.4-1. Gain (G) of the test case system as a function of scan angle in the plane of offset 
(0 O ) for optimized scanning using the correcting subreflector surface fitting error functional. 
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jidelobe level (SLL) of the test case system as a function of scan angle in the 
) or optimized scanning using the correcting subreflector surface fitting error 
functional. 
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• • i i y POL of the test case, shown in Figure 5. 1.4-3, is below -50 dB at 

The cross-polarization level, XPOL, ot me iesi, 

tbm. frequencies. The efficiency, BE, of the ,ee. case is shown in Figure 5.1.4-4. Since 

the 90% beam efficiency point is taken to be the limit of scan, the test case has a scan r g 
0 . 86 * at 20 GHz, 0.42* at 40 GHz, and 0.22* at 80 GHz. 

The aperture efficiency, c„, of the test case system is shown in Figure 5.1.4-5. Figure 
5 ,,4-fi shows the spillover efficiency, r.„, of the scanned test case system. The feed mi.poin.ing 
raulU Section 3.1.1 we,, used with the aperture efficiency and spillover efficiencies to 

approximate the illumination amplitude and phase efficiencies, and V of the test case 

. ,47 _ n( i 5 i 4 - 8 , suggest that most of the scan 
system. These results, shown m Figure 5.1.4-7 and 5.1.4 a, gg 

induced pattern degradation is caused almost exclusively by phase error as expec 

Figure. 5.4-0 through 5...4-16 .umma.it. the analysis results for the constant » scan 
paths. In general, these results show that the worst pattern degradation occurs for scan in the 

= 90* region. 


5.2 The Transmit Mode Raytracing Error Functional 

The transmit mode raytracing error functional determines the optimum position for the 
the cross-product of the rays exiting the system wl'h a unit vector in 
This is equivalent to minimizing the variation in the path 
the repositioned subreflector [5]. The 


subreflector by minimizing 
the desired scanned beam direction, 
lengths of the rays transmitted through the system using 
transmit mode raytracing error functional is defined 

M N 


‘tin _ is I r, 'j X r ° 

i=l j=l 


(5.2-1) 


= LEI 

when. r„j is the uni. vector of the id * ray in the system and f D is a uni. vector in the desired 
ccmned Iain beam direction. Unlike the correcting snbreflecto, surface titling error functional 
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Figure 5. 1.4-3. Cross-polarization level (XPOL) of the test case system as a function of scan 
angle in the plane of offset (0J for optimized scanning using the correcting subreflector surface 

fitting error functional. 
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Beam efficiency (BE), % 



Figure 5. 1.4-4. 
plane of offset 


Beam efficiency (BE) of the test case system as a function of scan angle in the 
( 6 0 ) for optimized scanning using the correcting subreflector surface fitting error 

functional. 
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Illumination amplitude effic 
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Figure 5.1 4-7. Illumination amplitude efficiency (e ) of the test case system as a function of 
scan angle in the plane of olTset (0J for optimized scanning using the correcting subrefiector 

surface fitting error functional. 
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Figure 5. 1.4-9. Gain (G) of the test case system as a function of <j> scan angle (<j> 0 ) for 6 — 0.5' 
and 0 o = 1.0* for optimized scanning using the correcting subreflector surface fitting error 

functional. 
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Sidelobe level (SLL), dB 



Figure 5.1.4-10. Sidelobe level (SLL) of the test case system as a function of <t> scan angle (<t> ? ) 
for 6 = 0.5’ and 0 o = 1.0* for optimized scanning using the correcting subreflector surface fitting 

error functional. 
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uTfIrV - n 3 V AP 7Tl f"-T y ° f the test Case system 88 a functi °n of 4> scan angle 
(<t>j o - °-5 and 0 O _ 1.0 for optimized scanning using the correcting subreflector surface 

fitting error functional. 
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Figure 5.1.4-14. Spillover efficiency (e ) of the test case system as a function of * scan angle 
(<t> ) for 0 o = 0.5* and 0 O = 1.0* for optimized scanning using the correcting subreflector surface 
° ° fitting error functional. 
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which emphasizes illumination amplitude error, this error functional evaluates only illumination 
phase error to provide an estimate of scanned pattern degradation. 


5.2.1 Transmit Mode Raytracing Optimization Implementation - TMRT1 


Optimization of the subreflector motions for scanned beams is much the same for the 
transmit mode raytracing error functional as for the correcting subreflector surface fitting error 
functional. The single major difference between the optimization procedures is the elimination 
of the outer path length optimization loop required by the correcting subreflector surface fitting 
error functional. The general scheme for optimization using the transmit mode raytracing error 
functional is shown in Figure 5.2. 1-1. First, the subreflector is positioned using (5. 1.2- 2a) 
through (5.1.2-2c) and a geometrical optics ray is traced through the system from the feed to the 
main reflector. Next, a ray is traced from the feed point to each of the grid points on the 
repositioned subreflector. The unit vector of this ray is 


r 3 

where 


f ^T_ x s' - yf ’ y S ' ~ 

i I r jI *’ W’ 1 



(5.2. 1-1) 


I ^3 1 — \/( x S ' * x f) 2 + (>s' " yf) 2 + ( z s' ~ z f) 2 * (5.2. 1-2) 

Snell’s law is then used with this unit vector and the unit normal of the subreflector to find the 


unit vector of the ray from the subreflector to the main reflector, ? r 


Unlike in the optimization using the Correction Subreflector Surface Fitting error 
functional, the main reflector is analytically defined for the transmit mode raytracing error 
functional. This allows an exact solution for the intersection of the ray from the subreflector 
with the parent paraboloid of the main reflector. The length of the ray from the subreflector to 
the main reflector, | R 2 (, can be found by using the quadratic formula 


Ip | _ \/b 2 - 4ac - b 
l R >l 2a 


(5.2. 1-3) 
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Figure 5.2. 1-1. Transmit mode ray tracing error functional optimization algorithm structure. 
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with 


& = 


(5.2.1-4a) 

(5.2.1-4b) 

(5.2.1-4c) 


a = r 2x + r 2y 

b = 2f 2x x s i -f 2r 2y y s , - 4f 2z F 
c = x g , 2 + y s , 2 - 4 z s # F. 

The intersection of the ray and the parent paraboloid of the main reflector is then 

^ = { X S' + r 2x| ^2 |* + r 2y | ^“2 I 1 Z S' r 2z|^2| |* (5.2. 1-5) 

The unit normal of the parent paraboloid of the main reflector, can now be found from 
(2.1-11). Using this unit normal with the unit vector ? 2 in SnelFs law gives the unit vector of 
the ray which is reflected from the main reflector, f r This unit vector is calculated for each ray 
which is traced through the system. 


The other necessary component of the transmit mode raytracing error functional is the 
unit vector of a ray exiting the system in the desired scan direction, ? D . This unit vector is 
given by the negative of (5. 1.1-1), the unit vector of a ray entering the system from the desired 
scan direction. After fj is found for each of the rays in the system, the error is found from 
(5.2-1). This error is used in t lie PowelPs method optimization which repositions the 
subreflector. 


5.2.2 Predicted Subreflector Motions 

The optimum subreflector positions were determined for the test case with the transmit 
mode raytracing error functional as with the correcting subreflector surface fitting error 
functional for three scan paths: a scan from 0 A = -1.0* to 0 = 1.0° for 6 = 0*, a scan from 6 - 
0 to <f> 0 = 180 for 6 q = 0.5, and a scan from = 0 to <j> Q = 180* for 0 o = 1.0". A significant 
decrease in the computational effort required was noted with a time of slightly less than 7 
minutes per scan direction and a total time of 4 hours 29 minutes for all 39 scan directions. 
This is over a six-fold reduction in the optimization time required. 
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The required x-, y-, and ^translations for the scan path in the plane of offset are shown 
in Figure 5.2.2- 1 as a function of scan angle in the plane of offset. The a- and ^-rotations 
required for the scan path are shown in Figure 5.2.2-2 as a function of scan angle in the plane of 
offset. The scan motions produced by the transmit mode raytracing error functional are larger 
because the subreflector is not constrained to be located along the central ray reflected from the 
main reflector as in the correcting subreflector surface fitting optimization. The x-. y-. and z- 
translations and o- and ^-rotations required for the constant 6 scan paths at 0 O = 0.5' and 
$ o = 1.0’ are shown in Figures 5. 2.2-3 and 5. 2. 2-4. 

5.2.3 Computational Results using TICRA GRASP7 Physical Optics Analysis 

The scanned test case systems created using the transmit mode raytracing error 
functional were analyzed as discussed in Section 5.1.4. Figure 5.2.3- 1 shows the gain, G, of the 
scanned test case system as a function of scan angle in the plane of offset, 0 o . The gain increases 
approximately 6 dB for each octave frequency increase as with the correcting subreflector surface 
fitting error functional, but, the loss during scan at each frequency is essentially independent of 
frequency for the transmit mode raytracing error functional. This difference indicates that the 
scan loss is probably not being caused by illumination phase error. The low overall sidelobe 
levels, SLL, shown in Figure 5.2.3-2 confirm this indication. Figure 5.2.3-3 shows the cross- 
polarization level, XPOL, as a function of scan angle in the plane of offset, 0 o . Figure 5.2.3-4 
shows the beam efficiency, BE, of the test case system when optimized for beam scanning using 
the transmit mode raytracing error functional. The beam efficiency produced with this 
optimization approach is comparable with that using the correcting subreflector surface fitting 
error functional of Section 5.1 at 20 GHz and much higher at the higher frequencies across the 
scan path. The scan range as defined by a minimum beam efficiency of 90% is 0.72° at 20 GHz, 
0.82* at 40 GHz, and 0.80“ at 80 GHz. This is a scan range of 35 half-power beam widths at 80 
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Figure 5.2.2- 1. Optimum subreflector x-, y-, and z-translations (x„ y„ z.) for the test case 
system as determ.ned using the transmit mode raytracing error functional t a function ofla 

angle in the plane of offset (0 o ). 
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Subreflector rotations (ar, pr), deg. 



Scan angle In the plane of offset (0o) f deg. 


Figure 5.2.2-2. Optimum subreflector a- and 0-rotations (or r , 0 r ) for the test case system as 
determined using the transmit mode raytracing error functional as a function of scan angle in 

the plane of offset (0 O ). 
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Figure 5 2.2-4. Optimum subreflector o- and 0-rotations (« r , 0 r ) for the test case system as 
determined using the transmit mode raytracing error functional as a function of * scan angle 

(<t> ) for 0 O = 0.5* and 0 O = 1.0*. 


Simp lified Error Functionals for Geometrical Optics Synthesis 


144 



Gain (G), dBi 



hgurp 5.2. 3-1. C.ain ((1) of the test case system as a function of scan angle in the plane of offset 
(0 O ) for optimized scanning using the transmit mode raytracing error functional. 
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Sidelobe level (SLL), dB 



Scan angle in the plane of offset (0o)> deg. 


Future 5.2. 3-2. Sidelobe level (SLL) of the test case system as a function of scan angle in the 
plane of offset ( 6 a ) for optimized scanning using the transmit mode raytrac.ng error functional. 


<t;— pl.r>«l Error Functional* for Geometrical Optica Syntbeai* 


146 




■I.o -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
Scan angle in the plane of offset (0 o ) f deg. 


Figure 5. 2 .3- 3 . Cross- polarization level (XPOL) of the test case system as a function of scan 
gle in the plane of offset (0J for optimized scanning using the transmit mode raytracing error 

functional. 
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Figure 5. 2. 3-4. Beam efficiency (BE) of the test case system as a function of scan angle s in , the 
plane of offset ( 0 O ) for optimized scanning using the transmit mode raytracing error functional. 
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Scan angle in plane of offset (Go), deg. 


Figure 5. 2.3-5. Aperture efficiency (f ap ) of the test case system as a function of scan angle in 
the plane of offset (0 O ) for optimized scanning using the transmit mode raytracing error 

functional. 
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GHz. 


The aperture efficiency, e ap , of the scanned test case system is shown in Figure 5.2.3-5 
as a function of scan angle in the plane of offset, 9 0 . The spillover effiency, r ap , given by 
GRASP7, is shown in Figure 5.2.3-6. As before, the results of Section 3.1.1 were used to 
generate approximate illumination amplitude and phase efficiencies, e Mnp and £*, for the 
system. These efficiencies are shown in Figures 5.2.3-T and 5.2.3-S. As expected, the scan loss 
in the test case systems synthesized using the transmit mode raytracing error functional is 
dominated by spillover and illumination amplitude losses. 

Figures 5.2.3-9 through 5.2.3-16 summarize the analysis results for the constant 9 scan 
paths. These results show that the worst pattern degradation occurs for scan to positive 9 0 
angles in the <t> 0 = O' region for the transmit mode raytracing error functional. This behavior is 
caused by the extreme motions selected by the error functional in this region as shown in 

Figures 5.2.2- 1 through 5.2. 2*3. 
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Spillover efficiency (esp), 



- 1.0 - 0.8 - 0.6 - 0.4 - 0.2 0.0 0.2 0.4 0.6 0.8 1.0 

Scan angle in the plane of offset (Go), deg. 


Figure 5.2.3-C. Spillover efficiency (c sp ) of the test case system as a function of scan angle in 
the plane of offset ( 0 O ) for optimized scanning using the transmit mode raytracing error 

functional. 
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Illumination amplitude efficiency (eamp), % 



Scan angle in the plane of offset (0o), deg. 


Figure 5.2.3-7. Illumination amplitude efficiency of the test case system as a function of 

scan angle in the plane of offset (0 O ) for optimized scanning using the transmit mode raytracing 

error functional. 
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3> scan angle W, deg. 


G ain (G) of the test case system as a function of <t> scan angle (*„) I for 0 O - 0.5* 
‘"and i 2T.C* to a£ild scanning using the transmit mode ray tracing error funct.onal. 
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Sidelobe level (SLL), dB 



F ig U re 5/2.3-10. Sidelobe level (SLL) of ll.e test case system as a function of * scan angle (* ) 
or 0 o - 0.5 and 6 o _ 1.0 for optimized scanning using the transmit mode ray tracing error° 

functional. 
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igure Beam efficiency (BE) of the test case system as a function of <t> scan angle (</> ) 

for 0 O - 0.5 and B Q = 1.0' for optimized scanning using the transmit mode ray tracing error” 

functional. 
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Figure 5.2.3-13. Aperture efficiency (c, p ) of the test case system as a function of * scan angle 
16 1 for 6 =0 5' and 6 =1.0* for optimized scanning using the transmit mode ray tracing error 
K 0> ° ' ° functional. 
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Figure 5.2.3-14. Spillover efficiency (<r sp ) of the test case system as a function of 0 scan angle 
(<U f° r K = 0-5" and 0 o = 1.0° for optimized scanning using the transmit mode ray tracing error 

functional. 


Simplified Error Functionals for Geometrical Optics Synthesis 


159 



Illumination amplitude efficiency (eamp), % 




Illumination pt 


20 GHz, 60 « 0.5* 
40 GHz, 6o - 0.5* 
80 GHz. 0o - 0.5* 
20 GHz, 0o - 1 .0* 
40 GHz, 0o - 1 .0* 
80 GHz, 0o - 1 .0* 
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Figure 5.2.3-16. Illumination phase efficiency (f J of the test case system as a function of <f> 
scan angle (4>J for 0 O = 0.5° and 0 O = 1.0° for optimized scanning using the transmit mode 

ray tracing error functional. 

Simplified Error Functional* for Geometrical Optics Synthesis 



References 


[1] T. Kitsuregawa, "Section 2.6 Steerable Bean. Antennas." TceA«o/o ff y » 

Satellite Communications Antennas, pp. 177-178, Artec!. House, Boston, MA, 1990. 


[ 2 ] 

[3] 

[4] 

[5] 


Personal conversation with Paul Werntz. 


W.H. Press, et al„ "10.5 Direction Set (Powell’s) Methods in Multidtoen.io».r 
Numerical Recipes, pp. 294-301, Cambridge University Press, New \ ork, M, 1990. 

W.H. Press, et al., "10.1 Golden Section Search in One Dimension," Numerical Recipes, 
pp. 277-282, Cambridge University Press, New York, NY, 1990. 


B. Shen, "2.3. Applications of Fermat’s Principle,” Multiple 
Ph.D. Dissertation, pp. 13-17, Virginia Polytechnic Institute 


Reflector Scanning Antennas, 

and State University, 1993. 


Sj^.pl.r.^1 Error Function*!* for Geometric*! Optics Synthesis 


1C2 



Chapter 6 

SUMMARY 


€.1 Conclusions 

The goal of this work was to develop a high gain, wide scanning reflector antenna 
system for use in a geostationary, microwave radiometric system for the Mission to Planet 
Earth. The design constraints imposed by the radiometric requirements of this mission are 
summarized in Table 1.1. Beam efficiency, given by (3. 1.1-7), is the far-field pattern 
characteristic that is of greatest concern during beam scanning. The beam efficiency must be as 
high as possible for this mission because radiometry measures thermal noise. This means that 
the ratio of desirable noise, the scene which is being observed, to undesirable noise, the rest of 
the universe, must be as large as possible. Also, beam efficiency should be constant with 

frequency because the radiometric measurements must be performed on several frequency bands 
simultaneously. 

The canonical offset Cassegrain reflector antenna configuration was chosen to 
complement the more complex structures being designed by Werntz [1] with the tri-reflector 
nfigurations, by Shen [2] using the spherical main reflector configurations, and by 
Takamizawa [3] using shaped-reflector configurations. The test case, summarized by Figure 4.2 
and Table 4.1, was chosen to fit a 42.52% scale test article developed at NASA Langley 
Research Outer to allow possible experimental verification of the resulting beam scanning 
system. This is the simplest system which is capable of scanning by subreflector motion only 
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, nd has an areal efficiency of greater .Iran MS. The « «- 

effiency of nearly 75% "zing » «»*“ M tha " “ ^ 

T »o geometrical optica WO error functionals were uud to generate the op, - 

scanning motions of the suhrefiector. The firs, o, these, the transmit mode ra, tracing error 
functional of Section 5.7. minimises the far-fteid pattern degradation by min, miring the 

minting between the rays ending the system and a -y in the desired scan direction. The 

r o upBVV for D\< = 709A, 18 HPBW 
test case system has a 90% beam efficiency scan range of 8 HPBW 

(oI Dm = 1417*. and 35 HPBW for D M = 2835* when optimized using this error funct.onal. 
optimization of the test rase system with the transmit mode raytracing error functional fo, a 
single scan direction requires just less than 7 minutes using a 50 MHz 80480 micro-compute,. It 
is believed that this is the fust application of an error function., of this type to the 
of subr.fl.ctor induced beam scanning in a dual reflector system. 

The correcting sub, dec, o, surface fitting emu, functional of Section 5.1, .fie, 
Kitsuregaw. (4), was the second error functional used to find the optimum scanning motions of 
the suhrefiector. This error function., minimize the far-field pattern degradation by fitttng the 
unscanned suhrefiector to a phase correcting suhrefiector synthesized fo, the scan direction and a 
given path length through the system. Optimization of the us, case system ustng tins error 
functional was over sin times slower at nearly 45 minutes per uan direction. This difference in 
optimization speed can he attributed to the addition., optimization loop required to vary the 
path length. The 90% beam efficiency scan range of the US, case when optimized nsmg the 
correcting suhrefiector surface fitting ecu, functional was found to he . HPBW fo, D M = 709* 
and 14 1 7 A and 10 HPBW for D M = 2835A. 

The Foldes Type 6 reflector antenna system, discussed in Section 3.4, is the other 
suhrefiector scanned, high gain reflector system which was investigated. Although, the Type 6 
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system beam scans by subreflector tilt only it was found to be capable of a 90% beam efficiency 
scan range of 19 HPBW for D M = 833A, 21 HPBW for D M = 1667A, and 31 HPBW for D M = 
3333A. This scan performance is significantly greater than the scan performance of the test case 
when optimized by either of the error functionals at the two smaller electrical aperture 
d.ameters. Because the Type G system has a higher subreflector eccentricity, one concern raised 
by its larger scan range is that the scan range of an offset Cassegrain reflector antenna may be 
highly affected by subreflector eccentricity. The transmit mode raytracing error functional does 
have a scan range which is comparable with that of the Type 6 system at. the largest electrical 
aperture size. Also, the lower phase error observed for the transmit mode raytracing error 
functional promises better scan performance with increasing aperture size than for the Type 6 

system. The scan characteristics of these three wide scanning reflector antenna system are 
summarized in Table 6.1-1. 

The optimization approach presented in Chapter 5 should be applicable to any dual 
fi ctor antenna system. In general, the transmit mode raytracing error functional should be 
used to optimize systems with electrically large apertures. Also, the transmit mode raytracing 
error functional is more suitable for radiometric systems because the beam efficiency of the 
resulting scanned system is more constant as a function of frequency. The correcting 
subreflector surface fitting error functional is better suited to smaller aperture antennas. The 
transition between the applicable ranges of these two error functionals as implemented here 
occurs for systems of D M from 800A to 1000 A. The correcting subreflector surface fitting error 
functional is better suited for use in the optimization of systems using array feeds to compensate 

for aperture phase error because this error functional primarily maximizes illumination 
amplitude and spillover efficiency. 

Other than the original geostationary radiometric mission, two possible uses for 
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subreflector induced beam scanning are tracking of geostationary satellites by large earth station 
antennas and reduction of the actuator demands imposed by over-zenith tracking in elevation 
over azimuth pedestals. A typical large earth station antenna in the 100’ class has a half power 
beam width of about 0.06’ at C-band. These earth stations must track geostationary satellites 
because most geostationary satellites have diurnal angular motions of greater than 0.06’. An 
antenna system capable of beam scanning over the limited tracking range required for this use 
would allow the use of smaller actuators for tracking by displacing the much smaller subreflector 
instead of the entire antenna. Over-zenith tracking in elevation over azimuth pedestals currently 
requires very rapid azimuth slew when the target passes through the zenith. A system capable 
of beam scanning by subreflector motion could reduce the maximum azimuth slew rate by using 
subreflector induced beam scanning to hold the beam on target near the zenith. This would 

allow the azimuth slew to occur more gradually during the time the target is near, rather than 
at, the zenith. 

The characteristics of five scanning antenna configurations of increasing complexity are 
summarized in Table 6.1-2. The prime-focus paraboloid with scanning achieved by lateral or 
optimal feed movement is included for comparison purposes [1], The Type 1 system is the test 
case system described here and is scanned by optimal subreflector positioning as determined 
with the transmit mode raytracing error functional. The Type 2 system is the Cassegrain II tri- 
reflector system of Werntz [I], The spherical system is the two subreflector and one caustic 
mirror with feed tilt system of Shen [2]. As shown, the scan range of reflector antenna systems 
generally increases with increasing complexity. On the whole, the Type 1 system described here 
provides good overall scan range while remaining fairly simple. Also, the Type 1 system’s areal 
and aperture efficiencies are significantly higher than those of the other two systems investigated 
here a Virginia Tech. 
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Table 6.1-2. Characteristics 


of five wide scanning reflector antenna systems 


Parameter 

Prime-focus Parab. 
Feed movement 

Type 1 
(LaPean) 

Type 2 
(Werntz) 

Spherical 

(Shen) 


lateral 

Optimal 

Main reflector diameter (D M ), A 


700 

2835 

480 

1200 x 1000 

F/D p 

1.000 

1.000 

0.515 

0.519 

0.260 

Areal efficiency ( — y— —r 5 ) 

-1- ... + L>S.l / 

1.0000 

1.0000 

0.9817 

0.8818 

0.9218 

Translational degrees of freedom 

1 

2 

3 

0 

0 

Rotational degrees of freedom 

1 

1 

2 

2 

feed: 1 

mirror: 2 

Gain (G), dBi 



77.78 

62.44 

63 

Aperture efficiency (e ap ), % 



75.63 

77.13 

50 

1 dB gain loss scan range, HPBW 

8 

34x0 

50 

30 x 60 

82 
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6.2 Future Work 


Future efforts to improve the scanning capabilities of offset canonical Cassegrain 
reflector antennas should focus on the following four main areas. 

• Improving the transmit mode raytracing error functional. 

The basic form of the transmit mode ray tracing error functional could be expanded to 
include illumination error effects as well as phase error effects. While this would 
necessarily cause the optimization to be performed for a given frequency of operation, 

many antenna systems, especially for civil communications, operate over only a narrow 
band of frequencies. 

• Improving the correcting subreflector surface fitting error functional. 

The implementation of the error functional could be extended to allow the grid of points 
to be of variable size. This improvement would allow a portion of the actual 
subreflector to be fitted to the correcting subrenector. This will possibly allow a closer 
fit as the compound curved surface move relative to each other. One potential 
complication to this scheme is that the illuminated portion of the actual subreflector 
may shrink below the region where geometrical optics techniques apply. 

• Investigating the possibility of using a simple array feed system to extend the scan range. 

A simple array feed could be used to reduce the aperture phase error. If the system was 
not limited by a low spillover or illumination amplitude efficiency this could 
significantly increase the scan range. As mentioned in Section 6.1, the correcting 
subreflector surface fitting error functional seems well suited for subreflector motion 
optimization for an array fed system because this error functional primarily optimizes 
the illumination amplitude and spillover efficiencies. Bailey [5] and Smith [C] have 
investigated the possibility of using array feed systems to improve the illumination 



phase efficiency of reflector antenna systems. 

. Investigating the effects of the geometry of Cassegrain reflector antenna systems on scan 
range. The results of Section 3.2 indicate that the offset height of the equivalent 
paraboloid of an offset Cassegrain reflector antenna has little effect on its beam scanning 
behavior. However, the effects of main reflector focal length, subreflector eccentricity, 
and other parameters were not investigated due to the pre-existing geometry 
requirements. The potential importance of this line of investigation is hinted at by the 
scan range of the Foldes Type 6 system when scanned by subreflector tilt only. 
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Appendix 1 

Program listing of the Correcting Subrcflector Surface Fitting error functional 

optimization software 
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**************************************************************** ******** 
* 

* Correcting Subreflector Surface Fitting Optimization Code 

* 

* James W. LaPean, Jr. 11/18/93 

***************************************************** *************** **** 


PROGRAM CSSF 
IMPLICIT NONE 

REAL* 8 SRFMD ( 4 ) , SRFM(25 , 25 , 6) , SRFA (25, 25 , 6) , SCND(6) , SCAN (2 ) , FOC, 
C FEED ( 3 ) , X ( 4 ) , SC ( 2 ) , XP ( 4 ) , PI , SRFS (25, 25 , 6) , SFOC ( 2 , 3 ) 

INTEGER*2 ICN , I , J , T, P, TINC, PINC 
CHARACTER* 12 SUBFILE 

COMMON /REFL/ SRFMD, SRFA, FOC, FEED, SC, SRFM, SRFS 
PI=3 . 14 1582 653 5897932 384 6D0 

* Program inputs read from input file 

OPEN ( 3 , FILE- ' KSFT1 . INP ' , STATUS- ' OLD ' ) 

READ ( 3 , * ) SRFMD ( 1 ) ,SRFMD(2) ,SRFMD(3) ,SRFMD(4) 

READ ( 3 , * ) FOC , FEED ( 1 ) , FEED (2) , FEED ( 3 ) 

READ ( 3 , * ) SCND ( 1 ) ,SCND(2) ,SCND(3) ,SCND(4) , TINC, PINC 
READ(3,*)X(1),X(4) 

READ(3, *) ICN 

READ ( 3 , ' (A12) ' ) SUBFILE 

X ( 3 ) —0 . DO 

X (2 ) —0 . DO 

* Nominal subreflector data entry 

OPEN ( 4 , FILE-SUBFILE , STATUS- ' OLD ' ) 

READ (4, ' (6F12.7) ') SFOC (1,1) , SFOC (1,2) , SFOC (1,3) , 

C SFOC (2,1) , SFOC (2,2) , SFOC (2,3) 

DO 1 1-1,25 
DO 1 J-1,25 

1 READ (4 , ' (6F12 . 7) ' ) SRFA (I , J, 1) , SRFA(I , J, 2 ) , SRFA I Z *3), 

C SRFA (I , J, 4 ) , SRFA (I , J, 5) , SRFA ' , o, 

CLOSE (4) 

* Main reflector array setup 

CALL MNREF 

* Perform scan optimized fitting for scan combinations 

WRITE (*, ' (//,A1) ') ' ' 

SCAN ( 2 ) —SCND ( 2 ) 

P=1 

2 SCAN (1) -SCND (1) 

T-l 

3 WRITE (* , ' (A25 , D8 . 3 ,«7 , D8 . 3 ) ' ) ' Calculating for Theta - ',SCAN(1), 

C ' phi - ' , SCAN (2) 

SC ( 1 ) —SCAN ( 1 ) 

SC (2) -SCAN (2) 

XP (1) —X ( 1 ) 

XP ( 2 ) —X (2) 

XP ( 3 ) —X ( 3 ) 

XP ( 4 ) —X ( 4 ) 

CALL PATHOPT( ICN ,XP, SFOC) 

ICN-ICN+1 

SCAN ( 1 ) -SCAN ( 1 ) +SCND ( 3 ) 

T-T+l 

IF (T.LE.TINC) GOTO 3 
SCAN ( 2 ) -SCAN ( 2 ) +SCND ( 4 ) 

P-P+l 

IF (P.LE.PINC) GOTO 2 
STOP 
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END 


********************************************************“********““** 
* Main reflector array filling subroutine 

.***...**.*****************«****•*************** 

SUBROUTINE MNREF 

REAL*8^SRFM^( 25 , 25,6) ,SRFMD(4) , DEN , F, SRFA ( 25 , 25 , 6 ) , FEED ( 3 ) , SC ( 2 ) , 

C SRFS (25,25,6) 

COMMON R /REFL/ SRFMD , SRFA , F , FEED , SC , SRFM , SRFS 
DO 1 1-1,25 

' J 5 2 1 - ( SRFMD ( 2 ) -SRFMD ( 1 ) ) * FLOAT ( 1-1 ) / 24 . DO+SRFMD ( 1 ) 

^RFM f I ' 3) - (SRFMD (4 ) -SRFMD(3) ) * FLOAT (J-l) /24 . DO+SRFMD ( 3 ) 
SRFM(I# <?' 1) — 2^50-1 MSRFM(I ,J, 2) **2 . D0+SRFM(I ,J,3)**2.D0) / F 
DEN-SQRT ( (SRFM(I , J , 2 ) / (2 . D0*F) ) **2 . DO 
C + (SRFM ( I , J , 3 ) / (2.D0*F) ) **2 . DO+1 . DO) 

SRFM ( I , J , 4 ) — 1 • DO / DEN 

SRFM(I,J,5)— SRFM ( I, J, 2) / (2.D0*F*DEN 
1 SRFM (I , J, 6) — SRFM(I, J, 3) / (2.D0*F*DEN) 

RETURN 

STOP 

END 


***************************************** 
Path length optimization subroutine 
***************************************** 


*********************** 




* 


****** 


****** 


* 


SUBROUTINE PATHOPT ( ICN , X , SFOC) 

REAL*8 1 SRFM 1 ! 2 5 , 25,6), SRFS (25,25,6) , FEED ( 3 ) , SCAN (2 ) , PI ,™C , PATH , 
C SRFA (25,25,6) ,TRANS(6) , R, C, TOL, X (4 ) ,EPS,F(2) ,SRFMD(4) , 

C SFOC (2,3) , XI ( 5 , 5 ) , DSCAN ( 2 ) 

INTEGER* 2 ICN, I , J , N, ITER 

COMMON /REFL/ SRFMD, SRFA, FOC, FEED, SCAN, SRFM, SRFS 
PI*3 . 14 159265 3 589793 2 384 6D0 


R-0.61803399D0 

C-1.D0-R 


TOL— 1 . D-10 

Optimization initialization 
SCAN (1) -SCAN (1) *PI / 1 • 8D2 
SCAN (2) -SCAN (2) *PI/1.8D2 
Xf3)-(X(4)+X(l) ) /2.D0 
X ( 2 ) — R*X ( 3 ) +C*X ( 1 ) 

PATH— X ( 2 ) 

CALL RAYTRC (PATH) 

DO 1-1,5 

TRANS ( I) =0. DO 
DO J-l, 5 

IF(I.EQ.J) THEN 
XI (I , J) —1. DO 
ELSE 
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XI ( I , J) =0 . Du 
ENDIF 
ENDDO 
ENDDO 
N=5 

CALL POWELL (TRANS, XI , N, TOL, ITER, EPS) 

F(1)=EPS 

PATH-X ( 3 ) 

CALL RAYTRC (PATH) 

DO 1*1,5 

TRANS (I) =0. DO 
DO J— 1 , 5 

IF(I.EQ.J) THEN 
XI(I, J)=1.D0 
ELSE 

XI ( I , J) =0 , DO 
ENDIF 
ENDDO 
ENDDO 
N=5 

CALL POWELL (TRANS , XI , N , TOL, ITER , EPS ) 

F(2)=EPS 

* Optimization loop 

2 IF (ABS(X(4) -X(I) ) . GT . TOL* (ABS (X ( 2 ) ) +ABS (X ( 3 
IF ( F ( 2 ) . LT . F ( 1 ) ) THEN 11 

X(1)=X(2) 

X(2)=X(3) 

X ( 3 ) =R*X (2 ) +C*X ( 4 ) 

F ( 1 ) =F (2 ) 

PATH-X ( 3 ) 

CALL RAYTRC (PATH) 

DO 1-1,5 

TRANS ( I) =0. DO 
DO J=l, 5 

IF(I.EQ.J) THEN 
XI (I , J) -1 . DO 
ELSE 

XI(I,J)-o.DO 
ENDIF 
ENDDO 
ENDDO 
N— 5 

CALL POWELL (TRANS , XI , N, TOL, ITER, EPS) 

F (2 ) —EPS 
ELSE 

X (4 ) —X ( 3 ) 

X ( 3 ) —X ( 2 ) 

X ( 2 ) — R*X ( 3 ) +C*X ( 1 ) 

F (2 ) — F ( 1 ) 

PATH-X (2) 

CALL RAYTRC (PATH) 

DO 1=1,5 

TRANS(I)=0.D0 
DO J— 1 , 5 

IF(I.EQ.J) THEN 
XI ( I , J) —1 . DO 
ELSE 

XI (I , J) =0 . DO 
ENDIF 
ENDDO 


)))) THEN 


Appendix I 


175 



VMirj: 

CALL POWELL (TRANS , XI , N , TOL, ITER, EPS) 

F(1)=EPS 
END IF 
GOTO 2 

ENDIF • ■ translation and rotation and sys. info 

* Output of optimum subreflector transia 

DSCAN ( 1) “SCAN ( 1) *1 . 8D2/PI 

CALL^ OUT ( FEED | TRANS , DSCAN , ICN , SFOC , PATH , EPS ) 

RETURN 

STOP 

END 


*************************— ***************************** 

Correcting subreflector generation subroutine 

**************************************** 


************** 

* 

* 

******************************** 


SUBROUTINE RAYTRC(PATH) 

C SRFA (25,25,6) 

JoMMON R / REFL/ SRFMD , SRFA , FOC , FFE Rl f ® J eCTOR FS 

* FIND COMPONENTS OF RAYS INC DENT ON MAIN REFLECTOR 

Rl(l)— COS (SCAN (1) ) 

R1 (2) =— SIN ( SCAN ( 1) ) *COS ( SCAN ( 2 ) ) 

Rl(3)— SIN (SCAN (1) ) *SIN (SCAN ( 2 ) ) 

DO 3 1=1,25 

. pihd 

R2 (1)-R1 (1) -2. *RDOTN*SRFM(I, J ,4) 

R2 (2)=R1 (2) -2 . *RDOTN*SRFM(I , J, 5) 

R2 ( 3 ) “R1 ( 3 ) -2 . *RDOTN*SRFM ( I , J , 6) 

* find 

R lL-COS(SCAN(in c ^(i))* si ^ ( s C AN (2)) * S RFM ( I,J,3)-SR FM( I,J,l)) 

* FIND SUBREFLECTOR POINTS BY SETTING TOTAL PATH LENGTH 

A*(R1 (2 ) -2 . DO*RDOTN*SRFM (I , J / 5) ) / 

C (R1 ( 1) -2 . DO*RDOTN*SRFM (I , J , 4) ) 

OR1 ( 3 ) -2 . DO*RDOTN*SRFM ( I , J# 6) 

B~C*SQRT ( 1+A**2 . DO) /SQRT ( 1-C**2 . DO) 

PROD=2 DO * (K*SQRT ( 1 +A* * 2 . DO+B* *2 . DO ) +A* ( SRFM ( I , J , 2 ) -FEED ( 2 ) ) + 
o B* ( SRFM ( I , J , 3 ) -FEED ( 3 ) ) ) 

SRFS(I,J,l)=(K**2.D0+SRFM(I,J,l)*PROD 

- ( FEED ( 2 ) -SRFM (I,J,2))**2. DO , D0 

- ( FEED ( 3 ) -SRFM (I , J, 3) ) **2 . DO+SRFM ( I , J , } > * * 

-FEED ( 1) **2 . DO) / (PROD+2 • DO*SRFM(I , J , 1) -2.D0*FEED(1) ) 

SRFS ( I , J , 3 ) =B* ( SRFS ( I , J , 1 ) -SRFM ( I , J , 1 ) ) +SRFM ( I , J , 3 ) 

* FIND SUBREFLECTOR NORMAL v _ ft 

R3L=SQRT( (SRFS (I, J, l)-FEED(l) )**2.D0 
c ' + (SRFS ( I , J , 2) -FEED (2 ) ) **2 . D0+ 


C 

C 

C 
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c ( SRFS ( I , J , 3 ) -FEED ( 3 ) ) **2 . DO) 

R3(l)=(FEED(l)-SRFS(I,J,l) ) /R3L 
R3 (2 ) = ( FEED ( 2 ) -SRFS (I,J,2))/R3L 
R3 ( 3 ) = ( FEED ( 3 ) -SRFS ( I , J , 3 ) ) /R3L 
SRFS ( I , J , 4 ) =R3 ( 1 ) -R2 ( 1 ) 

SRFS ( I , J , 5 ) =R3 (2)-R2(2) 

SRFS (I, J, 6) =R3 ( 3 ) -R2 ( 3 ) 

DEN=SQRT (SRFS ( I , J, 4 ) **2 . DO+SRFS (I,J,5)**2. DO+SRFS (I, J, 6)* *2. DO) 
SRFS (I , J, 4 ) =SRFS ( I , J, 4 ) /DEN 
SRFS (I, J, 5 ) =SRFS ( I , J , 5) /DEN 
3 SRFS ( 1 , 3 , 6 ) =SRFS ( I , J , 6 ) / DEN 

RETURN 
STOP 
END 


**************************************** 

* 

* Powell's Method Optimization 

* 

**************************************** 






SUBROUTINE POWELL ( P , XI , N , FTOL , ITER , FRET) 

IMPLICIT NONE 

REAL* 8 P (5) / XI ( 5 , 5 ) ,PT(5) ,PTT(5) ,XIT(5) , FTOL, FPTT, FRET, FP, T, DEL, 

C FUNC 

INTEGER* 2 I , J, ITER, N, IBIG, ITMAX 

ITMAX=200 

FRET=FUNC(P) 

DO J=1,N 
PT ( J) =P ( J) 

END DO 
ITER=0 

1 ITER=ITER+1 
FP=FRET 
IBIG=0 
DEL=O.DO 
DO 1=1 ,N 
DO J=1 , N 

XIT (J) =XI (J, I) 

END DO 
FPTT=FRET 

CALL LINMIN(P, XIT, FRET) 

IF (ABS (FPTT-FRET) .GT.DEL) THEN 
DEL=ABS ( FPTT-FRET) 

IBIG=I 
END IF 
END DO 

IF (2 . DO*DABS (FP-FRET) .LE. FTOL* (DABS(FP) +DABS(FRET) ) ) RETURN 
IF (ITER. EQ. ITMAX) PAUSE 'Powell exceeding maximum iteration' 

DO *J= 1 , N 

PTT ( J) =2 . DO*P ( J) -PT ( J) 

XIT ( J) =P ( J) -PT ( J) 

PT ( J) =P ( J) 

END DO 

FPTT=FUNC ( PTT) 

IF (FPTT.GE. FP) GOTO 1 

T*2 . DO* (FP-2 . DO*FRET+FPTT) * (FP-FRET-DEL) **2 . DO-DEL* ( FP-FPTT) **2 . DO 
IF (T.GT.O.) GOTO 1 
CALL LINMIN(P, XIT, FRET) 
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DO J=1,N 

XI(J,IBIG)=XIT(J) 


END DO 
GOTO 1 
END 


SUBROUTINE LINMIN (P , XI , FRET) 

IMPLICIT NONE 

REAL*8^P (5)^X1 ( 5) , AX , XX, FRET,TOL, PCOH ( 50) , XICOM (50) , BX, FA , FX , FB, 
c XMIN , BRENT 

INTEGER*2 J,NCOM 
COMMON /F1COM/ PCOM, XI COM, NCOM 
TOT,*! . D-10 


NCOM=5 
DO J*l,5 

PCOM(J)=P(J) 
XICOM (J) -XI (J) 
END DO 


AX-0. 


yv c 1 

CALL MNBRAK(AX,XX,BX,FA,FX,FB,F1DIM) 
FRET-BRENT (AX, XX, BX,FlDIM,TOL, XMIN) 


DO J-1,5 

XI(J)“XMIN*XI(J) 
P(J)-P(J) +XI (J) 
END DO 
RETURN 
END 


REAL*8 FUNCTION FIDIM(X) 

REAL^PCOM^O) , XICOM ( 50) , XT ( 50) , X , FUNC 
INTEGER*2 NCOM,J 

COMMON / F1COM/ PCOM, XICOM, NCOM 
DO J-l.NCOM 

XT ( J) -PCOM ( J) +X*XICOM ( J) 

END DO 

F1DIM-FUNC (XT) 

RETURN 

END 


SUBROUTINE MNBRAK ( AX , BX , CX , FA , FB , FC , FUNC) 

SS:S I L"S?CX.PX.FB.FC.DUH.GLIMIT f GOLD.TIN».R,Q.U.OLIM,F0 1 


C FUNC 

GOLD-1 . 618034D0 
GLIMIT-IOO.DO 
TINY-l.D-20 
FA- FUNC (AX) 

FB— FUNC (BX) 
IF(FB.GT.FA) THEN 
DUM-AX 
AX-BX 
BX-DUM 
DUM-FB 
FB-FA 
FA-DUM 
END IF 

CX-BX+GOLD* (BX-AX) 


Appendix 1 


178 


§23 


FC=FUNC(CX) 

1 IF (FB.GT.FC) THEN 

R=(BX-AX)*(FB-FC) 

Q= ( BX-CX) * (FB-FA) 

* R ’ ' <2 - * S1<=N ‘ MAX ' ABS< °-«> - TI »V) .Q-R, , 

IF ( (BX-U) * (U-CX) .GT. 0. ) THEN 
FU-FUNC(U) 

IF (FU.LT.FC) THEN 
AX=BX 
FA=FB 
BX=U 
FB=FU 
RETURN 

ELSE IF (FU.GT.FB) THEN 
CX=U 
FC=FU 
RETURN 
END IF 

U=CX+GOLD* (CX-BX) 

FU=FUNC(U) 

ELSE IF ( (CX-U) * (U-ULIM) .GT. 0. ) THEN 
U=ULIM 
FU=FUNC(U) 

ELSE 

U=CX+GOLD* (CX-BX) 

FU=FUNC(U) 

END IF 
AX=BX 
BX=CX 
CX=U 
FA=FB 
FB=FC 
FC=FU 
GOTO 1 
END IF 
RETURN 
END 


implicit^none°^ BRENT ( AX ' 1 BX ' Cx ' F , TOL, XMIN) 
c r “ l ‘ 8 

INTEGER*2 ITER, ITMAX 
ITMAX=100 


CGOLD=. 3819660D0 
ZEPS=1 . 0D-10 
A=MIN (AX, CX) 
B=MAX(AX,CX) 

V=BX 


W-V 


X=V 


E=0 . 

FX=F (X) 

=FX 
=FX 

ITER=1 , ITMAX 
XM=0 . 5* (A+B) 
TOLl»TOL*ABS (X) + ZEPS 
TOL2=2 . *TOLl 
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IF (ABS(X-XM).LE.(TOL2-.5*(B-A))) GOTO 3 
IF (ABS(E) . GT.TOL1) THEN 
R«(X-W)*(FX-FV) 

Q«(X-V) * (FX-FW) 
p«(X-V)*Q-(X-W)*R 
Q*2 . * (Q-R) 

IF (Q.GT.O) P— P 
Q-ABS(Q) 

ETEMP-E 

E-D 

IF (ABS (P) .GE. ABS( . 5*Q*ETEMP) . OR.P.LE.Q* (A-X) .OR. 
& P.GE.Q* (B-X) ) GOTO 1 

D-P/Q 
U-X+D 

IF (U-A. LT.TOL2 .OR. B-U. LT.TOL2) D*SIGN (TOL1 , XM-X) 
GOTO 2 
END IF 

1 IF (X. GE. THEN 

E-A-X 
ELSE 
E-B-X 
END IF 
D-CGOLD*E 

2 IF (ABS (D) .GE.TOL1) THEN 

U-X+D 

ELSE 

U-X+SIGN(TOLl,D) 

END IF 
FU=F(U) 

IF (FU.LE.FX) THEN 
IF (U. GE. X) THEN 
A-X 
ELSE 
B-X 

END IF 
V-W 
FV-FW 
W«X 
FW-FX 
X-U 
FX-FU 
ELSE 

IF (U.LT.X) THEN 
A-U 
ELSE 
B-U 

END IF 

IF (FU.LE.FW.OR.W.EQ.X) THEN 
V-W 
FV-FW 

w-u 

FW-FU 

ELSE IF ( FU . LE . FV . OR . V . EQ . X . OR . V . EQ . W ) THEN 
V-U 
FV*FU 
END IF 
END IF 
END DO 

PAUSE 'Brent exceed maximum iterations.' 

3 XMIN-X 
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BRENT=FX 

RETURN 

END 


************************************************************************ 

* 

* Powell's Method Optimization Error Function 

* 

************************************************************************ 

REAL*8 FUNCTION FUNC (TRANS) 

IMPLICIT NONE 

REAL*8 SRFA (25,25,6) , SRFAT (25,25,6) , TRANS ( 5 ) , RMSERR , SRFS (25,25,6) , 
C SRFMD ( 4 ) , FOC , FEED ( 3 ) ,SC(2) , SRFM (25,25,6) 

INTEGER*2 I,J 

COMMON /REFL/ SRFMD , SRFA, FOC , FEED, SC , SRFM, SRFS 
CALL POSIT (SRFAT, TRANS) 

RMS ERR* 0 . DO 
DO 1 1*1,25 
DO 1 J=1 , 25 

RMSERR*RMSERR+ ( ( (SRFAT (I , J , 1) -SRFS (I , J, 1) ) *SRFS(I, J,4) ) **2.D0 
C + ( (SRFAT (I, J, 2) -SRFS (I, J, 2) ) *SRFS ( I , J, 5) ) **2.D0 

C +( ( SRFAT ( I , J, 3 ) -SRFS (I , J, 3 ) ) *SRFS(I, J, 6) ) **2 .DO) 

1 CONTINUE 

FUNC=DSQRT (RMSERR) 

RETURN 

STOP 

END 




* Type 1 Concept Subreflector Positioning Code 

* 

************************************************************************ 


SUBROUTINE POSIT (SRFAT , TRANS) 

IMPLICIT NONE 

REAL* 8 SRFA (25, 25, 6) , SRFAT (2 5 , 25 , 6) , TRANS (5) , SRFMD (4) , FOC, FEED (3) , 
C SCAN (2) ,T(5) , SRFM (25, 25, 6) , SRFS ( 25 , 25 , 6) 

INTEGER*2 I,J 

COMMON /REFL/ SRFMD , SRFA , FOC , FEED , SCAN , SRFM, SRFS 
T ( 1 ) “TRANS ( 3 ) 

T ( 2 ) =TRANS ( 4 ) 

T ( 3 ) “TRANS ( 5 ) 

T ( 4 ) “TRANS ( 1 ) 

T ( 5 ) “TRANS ( 2 ) 

DO 1 I“1 , 25 
DO 1 J»1 , 2 5 

SRFAT ( I, J, 1 ) =T ( 1 ) +SRFA (13,13,1)- 

SIN(T(4) ) *( SRFA ( I , J, 2 ) -SRFA( 13,13,2) )+ 

COS (T ( 4 ) ) *SIN (T (5) ) * ( SRFA ( I , J , 3 ) -SRFA ( 13 , 13 , 3 ) ) + 

COS (T ( 4 ) ) *COS (T(5) )*( SRFA ( I , J , 1) -SRFA (13,13,1) ) 

SRFAT (I , J, 2 ) “T (2 ) +SRFA (13,13,2)+ 

COS (T (4 ) ) * ( SRFA (I , J , 2) -SRFA (13 ,13,2))+ 

SIN (T ( 4 ) ) *SIN (T ( 5) ) * (SRFA ( I , J , 3 ) -SRFA ( 13 , 13 , 3 ) )+ 

SIN (T ( 4 ) ) *COS (T ( 5) ) *(SRFA(I,J, 1) -SRFA ( 13 , 13 , 1 ) ) 

SRFAT ( I , J , 3 ) -T ( 3 ) +SRFA (13,13,3)+ 

COS (T(5))*( SRFA ( I , J , 3 ) -SRFA (13,13,3))- 
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c 


SIN(T(5))*( SRFA( I , J , 1) -SRFA (13,13,1' ) 
SRFAT ( I , J , 4 ) *-SI J(T(4) ) *SRFA(I, J, 5) + 

C COS (T ( 4) ) *SIN(T(5) ) *SRFA ( I , J, 6) + 

C COS (T ( 4 ) ) *COS (T ( 5 ) ) *SRFA ( I , J, 4 ) 

SRFAT (I , J, 5) “COS (T(4 ) ) *SRFA ( I , J, 5) + 

C SIN (T ( 4 ) ) *SIN(T(5) ) *SRFA(I r J, 6) + 

C SIN (T (4 ) ) *COS(T(5) ) *SRFA(I, J,4) 

1 SRFAT (I , J , 6) “COS (T( 5) ) *SRFA(I , J, 6) - 

C SIN (T ( 5 ) ) *SRFA ( I , J , 4 ) 

RETURN 

STOP 

END 


************************************************************************ 

* 

* Type 1 Concept Output Code 

* 

************************ ************************ ********* *************** 


SUBROUTINE OUT ( FEED , TRANS , SCAN , ICN , SFOC, PATH , RMSERR) 

IMPLICIT NONE 

REAL* 8 FEED (3) , TRANS (5) , SCAN ( 2 ) , SRFAT ( 25 , 25 , 6) , PI , PATH , RMSERR, 

C SFOC (2 r 3) , GCOEF (10) ,U,V 

INTEGER* 2 ICN,K(3) 

CHARACTER* 12 OUTFILE 

PI-3 . 14 1592653 5897932384 6D0 

U— SIN (SCAN ( 1 ) *PI/1.8D2) *COS (SCAN (2)*PI/1.8D2) 

V— SIN ( SCAN ( 1 ) *PI/1.8D2) *SIN (SCAN (2 ) *PI/ 1 • 8D2 ) 

K(l)-48+ICN/100 

K(2)«48+(ICN-(K(1)-48)*100) /10 
K(3)«48+ICN-(K(l)-40) *100- (K(2) -48) *10 
CALL POSITF (SFOC , TRANS , GCOEF) 

CALL POSIT (SRFAT, TRANS) 

OUTFILE- 'KSFT1 XXX . MOV ' 

OUTFILE (6:6) -CHAR ( K ( 1 ) ) 

OUTFILE (7:7) -CHAR ( K ( 2 ) ) 

OUTFILE (8 : 8) -CHAR(K (3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

WRITE (4 , 9 (A10, 2F12 . 7) ' ) ' SCANNED TO ' , SCAN ( 1) , SCAN (2) 

WRITE (4 , * (5F15.10) ')TRANS(4) # TRANS(5) ,TRANS(3) ,TRANS(1) , TRANS ( 2 ) 

WRITE (4 , 9 (2F15. 10) ') PATH , RMSERR 

ENDFILE(4) 

CLOSE (4) 

* 20 GHz GRASP7 Input File 

OUTFILE- ' KSFT 1XXX . INP ' 

OUTFILE (6:6) — CHAR (K ( 1 ) ) 

OUTFILE (7:7) — CHAR (K ( 2 ) ) 

OUTFILE (8 : 8) — CHAR(K{3) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

WRITE(4 , 9 (A30) ') ****** GRASP77 INPUT FILE ****** 

WRITE (4 , 9 (A10,2F12.7) ' ) 'SCANNED TO ' ,SCAN(1) ,SCAN(2) 

WRITE (4, * ( F12 .7,412) ' ) 1 . 5D-2 , 2 , 1 , 0 , 0 

WRITE (4, * (A31) ' ) *** MAIN REFLECTOR INPUT DATA **' 

WRITE (4 , 9 (3F12.7) ')0. f 0. ,0. 

WRITE (4, * (3F12.7) ')1.,0.,0. 

WRITE (4, * (3F12.7) * )0. , 1. ,0. 

WRITE (4, ' (I2,2F12.7) ') 1,7.795,0. 

WRITE ( 4 , * (312) * ) 4 , 0 , 0 

WRITE (4, * (4F12.7) * ) 0 . , 0 . , 0 . , 13 . 5 
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WPTTF (A ' (2F12 .7) ' J ^ Jl c 5.315 

WRITE ( 4 | ' (A30) ' ) ' * * Sr REFLECTOR INPUT DATA **' 

WRITE ( 4 , ' (3F12.7) ')0. ,0. ,0. 

WRITE ( 4 , ' ( 3F12 .7) ' ) 1. ,0. ,0. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' (12, 2F12.7) ') 1, ( SRFAT (25,13,2) +SRFAT (1,13,2) )/2.D0, 

C ( SRFAT (13,25,3) +SRFAT (13,1,3) ) /2 .DO 

WRITE ( 4 , ' (212) ') 3, 0 

WRITE ( 4 , ' (5F14.7) ')GCOEF(l) / GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 
WRITE (4, ' (5F14.7) ')GCOEF(6) ,GCOEF(7) ,GCOEF(8) ,GCOEF(9) ,GCOEF(10) 
WRITE (4 , ' (5F12 .7, 15) ')0.,0.,0.,0.,0.,1 

WRITE (4, ' (2F12 .7) ' ) ( SRFAT ( 2 5 , 13 , 2 ) -SRFAT ( 1 , 13 , 2 ) ) /2 . DO , 

C ( SRFAT (13,25,3) -SRFAT (13,1,3) ) /2.D0 

WRITE (4, ' (A21) ') '** FEED INPUT DATA **' 

WRITE ( 4 , ' (3F12.7) / )FEED(2) , FEED (3) , FEED ( 1) 

WRITE (4 , ' (3F12.7) ') .964483694,0. ,-.264142395 
WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' (12;' ) 1 

WRITE (4, ' (6F12.7) 'JO.,0.,0.,0.,0.,0. 

WRITE (4, ' (F12.7, 12) ' ) 0. ,1 

WRITE (4, ' (4F12.7,3I2) ' ) 0 . ,0. ,1. ,90. ,0,0,3 
WRITE ( 4 , ' (12) ')6 

WRITE (4, ' ( 5F12 .7,13) ')-15. , -15. , 7 . 1574 , 0 . , 0. , 1 

WRITE (4, ' (A40) ' ) ' ** MAIN REFLECTOR FIELD SPECIFICATION **' 

WRITE (4 , ' (12) ' ) 1 

WRITE (4, ' (3F12.7) ' ) 7 . 7 95 , 0 . , 10 . 

WRITE ( 4 , ' (3F12.7) ') 1. ,0. ,0. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' (216, F12. 7,213) ') 60,144,0. ,0,0 
WRITE ( 4 , ' (313) 7 ) 0, -1, 1 
WRITE ( 4 , ' (F12.7) ') 10. 

WRITE ( 4 , ' (2F12.7) ' )U, V 
WRITE ( 4 , ' (213) ')0,1 
WRITE ( 4 , ' (213) ') 3, 1 

WRITE (4, ' ( 4F12 .7,216) ' ) -1 . D-2 , -1 . D-2 , 1. D-2 , 1 .D-2 , 25 , 25 
WRITE ( 4 , ' (416) ' ) 3 , 1 
OUTFILE* ' KSFT1XXX . PI 9 
OUTFILE (6:6) =CHAR (K ( 1 ) ) 

OUTFILE (7:7) =CHAR ( K ( 2 ) ) 

OUTFILE (8:8) =CHAR ( K ( 3 ) ) 

WRITE ( 4 , ' (A12) ' )OUTFILE 

WRITE (4 , '(A3 9) ' ) ' ** SUB REFLECTOR FIELD SPECIFICATION **' 

WRITE ( 4 , ' (12) ')0 

WRITE (4 , 9 (A30 ) ' ) ' *+ FEED FIELD SPECIFICATION 
WRITE ( 4 , ' (12) ')0 
ENDFILE (4 ) 

CLOSE (4) 

OUTFILE* ' KSFUVXXX . INP ' 

OUTFILE (6:6) *CHAR ( K ( 1 ) ) 

OUTFILE (7:7) *CHAR ( K ( 2 ) ) 

OUTFILE (8:8) “CHAR ( K ( 3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS* ' UNKNOWN ' ) 

WRITE ( 4 , 9 ( A2 5 ) ' ) 9 UVPROC CONTROL INPUT FILE' 

OUTFILE* 'KSFT1XXX. PI' 

OUTFILE (6:6) *CHAR ( K ( 1 ) ) 

OUTFILE (7:7) *CHAR ( K ( 2 ) ) 

OUTFILE (8:8) =CHAR ( K ( 3 ) ) 

WRITE ( 4 , ' (A12) ' JOUTFILE 
WRITE (4, ' (612) ')1, 1,1, 1,0,0 
WRITE ( 4 , ' (212) ' ) 1, 1 
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WRITER, ' <61-; ,o, i,u,0, 1,0 
WRITE ( 4 , ' ( A12 ) ' ) ' TEMP. P2 ' 

WRITE (4 , ' (12) ')3 
WRITE ( 4 , ' (214) ' ) 101 , 101 
WRITE (4, ' (712) ') 0,0, 1,1, 0,0,0 
WRITE (4 , ' (412) ')3, 0,2,0 
WRITE (4 , ' (F5 . 1) ' ) 3 . 

WRITE ( 4 , ' (F5.1) ') 10. 

WRITE(4, ' (512) ')!, 0,0,1, 1 
WRITE (4, ' (F7. 4) ') 0.1213 
ENDFILE (4 ) 

CLOSE (4) 

* 40 GHz GRASP7 Input File 

OUTFILE- ' KSFT1XXX . INP ' 

OUTFILE(6:6)=CHAR(K(l)+l) 

OUTFILE (7 : 7)— CHAR (K (2) ) 

OUTFILE (8 : 8 ) — CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS** ' UNKNOWN' ) 

WRITE(4, ' (A30) ') ****** GRASP77 INPUT FILE *****' 
WRITE(4, ' (A10,2F12. 7) ') 'SCANNED TO ' , SCAN(l) , SCAN(2) 
WRITE (4, ' (F12. 7,412) ' ) 7 . 5D-3 , 2 , 1 , 0 , 0 
WRITE (4, ' (A31) ') ' * * MAIN REFLECTOR INPUT DATA **' 


WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 

WRITE (4, 
WRITE ( 4 , 
WRITE (4 , 
WRITE (4, 
WRITE (4, 

* 

WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4, 


(3F12.7) ')0. ,0. ,0. 

(3F12.7) ') 1. ,0. ,0. 

(3F12.7) ') 0. ,1. ,0. 

(12 , 2F12 . 7 ) ') 1,7. 795,0. 

(312) ' ) 4 , 0, 0 

(4F12.7) ')0. ,0. ,0. ,13.5 

(2F12.7) ')5. 315, 5. 315 

(A30) ' ) ' ** SUB P.EFLECTOR INPUT DATA **' 

(3F12.7) ')0. ,0. ,0. 

( 3F12 .7) , )1.,0.,0. 

(3F12.7) ' ) 0 . ,1. ,0. 

(12 , 2F12 . 7 ) ' ) 1 , ( SRFAT (25,13,2) +SRFAT (1,13,2) ) /2. DO , 

( SRFAT (13,25,3) +SRFAT (13,1,3) ) /2.D0 

(212) ' ) 3, 0 

(5F14.7) ')GCOEF(l) ,GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 
(5F14.7) ' ) GCOEF (6) ,GCOEF(7) ,GCOEF(8) ,GCOEF(9) ,GCOEF(10) 
(5F12.7,I5) ' ) 0 . ,0. ,0. ,0. ,0. ,1 

(2F12.7) ') (SRFAT (25,13,2) -SRFAT(1, 13,2) )/2. DO , 

(SRFAT (13 ,25,3) -SRFAT (13 , 1, 3) ) /2 . DO 
(A21)')'** FEED INPUT DATA **' 

(3F12 . 7) ' ) FEED (2) , FEED(3) , FEED(l) 

(3F12.7) ') .964483694,0. ,-.264142395 
(3F12.7) ')0. ,1. ,0. 


WRITE (4 , ' (12) ') 1 


WRITE (4 ,'(6F12.7)')0.,0.,0.,0.,0.,0. 

WRITE (4, ' (F12 .7,12) ')0. ,1 

WRITE (4 , ' (4F12.7,3I2)')0.,0.,1.,90.,0,0,3 

WRITE(4, ' (12) ' ) 6 

WRITE (4 , ' (5F12.7,I3) ' ) -15 . ,-15. ,7.1574,0. ,0. ,1 

WRITE (4 , ' (A40) ') MAIN REFLECTOR FIELD SPECIFICATION **' 

WRITE(4, ' (12) ') 1 

WRITE(4, ' (3F12.7) ') 7. 795,0. ,10. 

WRITE(4 , ' (3F12.7) ' ) 1 . ,0. ,0. 

WRITE (4 , ' (3F12.7) ')0. ,1. ,0. 

WRITE (4 , ' (216, F12. 7, 213) ') 60, 144,0. ,0,0 
WRITE (4 , ' (313) ') 0,-1, 1 
WRITE ( 4 , ' (F12.7) ') 10. 

WRITE (4, ' (2F12.7) ')U,V 
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WRITE (4 , ' (213) ')0,1 
WRITE ( 4 , ' (213) ') 3, 1 

WRITE (4, ' (4F12.7,2I6) ' ) -5 . D-3 , -5 . D-3 , 5 . D-3 , 5 . D-3 25 25 
WRITE(4, ' (416) ' ) 3, 1 ' ' 

OUTFILE- ' KSFT1XXX . PI ' 

OUTFILE ( 6 : 6) =CHAR (K ( 1 ) +1 ) 

OUTFILE (7:7) =CHAR (K ( 2 ) ) 

OUTFILE (8:8) — CHAR (K (3 ) ) 

WRITE (4, ' (AI2 ) ') OUTFILE 

SUB REFLECT0R FIELD SPECIFICATION **' 

WRITE (4 ; ' (12 )*' ) o'** FEE ° FIEL ° SPECIFICATION **' 
ENDFILE(4) 

CLOSE (4) 

OUTFILE* ' KSFUVXXX . INP ' 

OUTFILE (6:6) — CHAR (K ( 1 ) +1 ) 

OUTFILE (7:7) -CHAR (K ( 2 ) ) 

OUTFILE (8:8) — CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE=OUTFILE , STATUS* ' UNKNOWN ' ) 

WRITE(4 , ' (A25) ' ) 'UVPR0C CONTROL INPUT FILE' 

OUTFILE-' KSFT1XXX. PI' 

OUTFILE (6:6) =CHAR (K ( 1 ) +1 ) 

OUTFILE (7 : 7 ) —CHAR (K (2 ) ) 

OUTFILE (8:8) =CHAR (K ( 3 ) ) 

WRITE (4, ' (A12 ) ') OUTFILE 
WRITE (4 , ' (612) ')1, 1,1, 1,0,0 
WRITE(4, ' (212) ' ) 1, 1 

WRITE(4, '(612) ')0, 1,0, 0,1,0 

WRITE(4, ' (A12) ') ' TEMP. P2 ' 

WRITE (4, ' (12) ')3 
WRITE (4 , ' (214) ')101,101 
WRITE (4 , ' (712) ')0, 0,1, 1,0, 0,0 
WRITE (4, '(412) ')3, 0,2,0 

LTD TPPiri / J / /tic . 


WRITE ( 4 , ' (F5 . 1) ') 3. 

WRITE ( 4 , ' (F5 . 1) ') 10. 

WRITE (4 , '(512) ')1, 0,0, 1,1 
WRITE (4 , ' (F7.4) ') 0.0606 
ENDFILE(4) 

CLOSE (4) 

80 GHz GRASP7 Input File 
OUTFILE* 'KSFT1XXX. INP' 

OUTFILE ( 6 : 6) =CHAR (K ( 1 ) +2 ) 

OUTFILE (7:7) -CHAR ( K ( 2 ) ) 

OUTFILE (8:8) —CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

WRITE (4 , ' (A30) ') '***** GRASP77 INPUT FILE ****** 

KE, ( i:-'!??;L'i;:;o“j N ,r LECTOR 1NPUT mta **' 

WRITE (4, ' (3F12.7) ') 1. ' 0. 'o. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4 , ' (I2,2F12.7) ') 1,7.795,0. 

WRITE (4, '(312) ' ) 4 , 0 , 0 

WRITE(4,' (4F12.7) ')0. ,0. ,0. ,13.5 

WRITE (4 , ' (2F12.7) ' ) 5 . 3 15 , 5 . 3 15 

' (A30) '** SUB Reflector INPUT DATA **' 
WRITE (4 , ' (3F12.7) ')0. ,0. ,0. 

WRITE (4, ' (3F12.7) ') 1. ,0. ,0. 
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WRITE (4 
WRITE (4 
C 

WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
C 

WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
WRITE (4 
OUTFILE 
OUTFILE 
OUTFILE 
OUTFILE 
WRITE (4 
WRITE (4 


(3F12 .7)')0.,1.,0. 

( 12 2F12 . 7 ) ' ) 1 , ( SRFAT (25,13 , 2) +SRFAT( 1 ,13,2))/2.D0, 

( SRFAT (13,25,3) +SRFAT (13,1,3) )/2. DO 

(212) ' ) 3 0 

5F14.7) ')GCOEF(l) ,GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 
(5F14.7)') GCOEF ( 6 ) ,GCOEF(7) ,GCOEF(8) ,GCOEF(9) ,GCOEF(10) 
(5F12.7,I5) ')0. ,0. ,0. ,0. ,0. ,1 

(2F12.7) ') ( SRFAT (25,13,2) -SRFAT (1,13,2) )/2.D0, 

( SRFAT (13,25,3) -SRFAT (13,1,3) ) /2. DO 
(A21)')'** FEED INPUT DATA **' 

(3F12.7) ' ) FEED(2) , FEED ( 3 ) , FEED ( 1) 

(3F12.7) ') .964483694,0. ,-.264142395 
(3F12 .7) ' ) 0. , 1. , 0. 

(12) ')1 

(6F12.7) ')0.,0.,0.,0.,0.,0. 

(F12.7.I2) ')0. ,1 

(4F12 .7,312) ')0.,0.,1.,90.,0,0,3 
(12 ) ' ) 6 

(5F12 .7,13) ') -15. ,-15. ,7.1574,0. ,0. ,1 
(A40) ' ) ' ** MAIN REFLECTOR FIELD SPECIFICATION **' 

(12) ')1 

(3F12.7) ')7. 795,0. ,10. 

(3F12.7) ')1. ,0. ,0. 

(3F12.7) ')0. ,1. ,0. 

(216 , F12 .7,213) ' ) 60 , 144 , 0 . , 0 , 0 
(313) ')0,-l, 1 
(F12.7) ' ) 10. 

(2F12.7) ' )U,V 

(213) ')0,1 
(213) ' ) 3 , 1 

( 4F12 . 7 , 216) ' ) -2 . 5D-3 , -2 . 5D-3 , 2 . 5D-3 , 2 . 5D-3 ,25,25 
(416) ')3,1 
KSFT1XXX. PI ' 

(6:6)«CHAR(K(l)+2) 

(7:7) “CHAR ( K ( 2 ) ) 

(8:8) -CHAR ( K ( 3 ) ) 

, ' (A12) ' )OUTFILE 

'(A39)')'** SUB REFLECTOR FIELD SPECIFICATION **' 


WRITE ( 4 , ' (12) ')0 

WRITE (4 , ' (A30) ') '** FEED FIELD SPECIFICATION **' 

WRITE ( 4 , ' (12) ')0 

ENDFILE(4) 

CLOSE (4) 

OUTFILE- • KSFUVXXX . INP • 

OUTFILE ( 6 : 6) -CHAR (K ( 1 ) +2 ) 

OUTFILE (7 :7)— CHAR(K(2) ) 

OUTFILE (8 : 8) -CHAR(K (3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

WRITE (4, ' (A25) ') ' UVPROC CONTROL INPUT FILE' 
OUTFILE- ' KSFT1XXX . PI ' 

OUTFILE (6: 6) -CHAR (K(l) +2) 

OUTFILE (7:7) —CHAR (K ( 2 ) ) 

OUTFILE (8:8) — CHAR (K ( 3 ) ) 

WRITE ( 4 (A12) ') OUTFILE 

WRITE (4, ' (612) ')1, 1,1, 1,0,0 
WRITE (4, ' (212) ' ) 1, 1 
WRITE (4, ' (612) ') 0,1, 0,0, 1,0 
WRITE (4 , ' (A7) ') 'TEMP.P2' 

WRITE (4 , ' (12) ' ) 3 
WRITE (4, ' (214) ') 101, 101 
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WRITE (4 , ' (712) ')0, 0,1, 1 , 0 , 0,0 
WRITE ( 4 , ' (412) ')3, 0,2,0 
WRITE ( 4 , ' (F5.I) ')3. 

WRITE ( 4 , ' (F5 . 1) ') 10. 

WRITE (4, ' (512) ')1, 0,0, 1,1 
WRITE ( 4 , ' (F7.4) ' ) 0 . 0303 
ENDFILE(4) 

CLOSE (4) 

RETURN 

STOP 

END 


********************* *************************************************** 
* T yP e 1 Concept Subreflector Focal Point Positioning Code 

****************************** t ** 1H „*^ t „ ltltt4114t((t0(ttt(ttt4itt ^^ 

SUBROUTINE POSITF (SFOC , TRANS , GCOEF) 

IMPLICIT NONE 

REAL*8 SRFA (25,25,6) ,SFOC(2,3) ,TRANS(5) ,SRFMD(4) ,FOC,FEED(3) , 
n ,T{5) ' CCOEF (10) , SFOCT (2,3) , ALPHA , BETA , GAMMA , 

C DELTA , C , SRFM (25,25,6) , SRFS (25 , 25, 6) 

INTEGER *2 I ' ' 

COMMON^/REFL/ SRFMD, SRFA, FOC, FEED, SCAN, SRFM, SRFS 

T ( 2 ) =TRANS ( 4 ) 

T ( 3 ) “TRANS ( 5 ) 

T ( 4 ) “TRANS ( 1 ) 

T ( 5 ) “TRANS ( 2 ) 

DO 1 1=1,2 

SFOCT (1,1) =T ( 1 ) +SRFA(13, 13, 1)- 

SIN(T(4) ) * (SFOC (1,2) -SRFA (13,13,2) ) + 

COS ( T ( 4 ) ) *SIN(T(5) ) * ( SFOC ( I , 3 ) -SRFA ( 13 , 13 , 3 ) ) + 

c. ‘ 0! \ (T < 4 > } * cos (T(5) ) *(SFOC(I, 1) -SRFA (13, 13,1) ) 

SFOCT (1,2) =T ( 2 ) +SRFA (13,13,2)+ ' 

OT „, m COS (T (4 ) ) * (SFOC (1,2) -SRFA (13,13,2)) + 

SIN(T(4) ) *SIN(T(5) ) * ( SFOC (1,3) -SRFA (13,13,3) )+ 

? IN 5 T(4,) * COS(T(5,, * (SFOC < I ' 1 )- SRFA (“. 13.D) 

SFOCT ( 1 , 3 ) “T ( 3 ) +SRFA ( 1 3 , 1 3 , 3 ) + 

COS (T ( 5) ) * ( SFOC (1,3) -SRFA (13,13,3))- 

C*DSQRT{ ( SRFA (13,13, lf-SFOC(2 *if 1 

c (SRFA (13, 13,2) -SFOC (2, 2 ) ) ** 2 .D 0 + 

C (SRFA ( 13 , 13,3) -SFOC (2 , 3) ) **2 . DO) - 

C DSQRT( (SRFA(13, 13,1) -SFOC(l, 1) ) **2.D0+ 

C (SRFA (13,13,2) -SFOC(l , 2) ) **2 . D0+ 

C (SRFA ( 13 , 13,3) -SFOC ( 1 , 3) ) **2 . DO) 

BETA=SFOCT (2,2) -SFOCT (1,2) 

GAMMA=S FOCT (2,3) -SFOCT (1,3) 

DELTA=SFOCT (2,1) -SFOCT (1,1) 

ALPHA=SFOCT ( 1 , 1 ) * * 2 . DO+SFOCT ( 1 , 2 ) * *2 . DO-t-SFOCT ( 1 , 3 ) **2 . D0- 
C SFOCT (2 , 1 ) **2 . DO-SFOCT ( 2 , 2 ) * *2 . DO-SFOCT (2 , 3 ) **2 . DO- 

C* * 2 . DO 

GCOEF (1) *=4 . DO* (BETA* *2 . D0-C**2 . DO) 

=8 . DO* BETA* GAMMA 
=4 .DO* (GAMMA* *2 . D0-C**2 . DO) 

=4 . DO* ( ALPHA* BETA+2 . D0*C**2 . DO* SFOCT (2 2 ) ) 
e 4 .DO* (ALPHA*GAMMA+2 . D0*C**2 . DO* SFOCT ( 2 , 3 ) ) 


C 

C 


c 

c 


c 

c 


GCOEF (2) 
GCOEF (3) 
GCOEF (4) 
GCOEF (5) 
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GCOEF (6) “ALPHA* *2 . DO-4 . D0*C**2 . DO* (SFOCT(2 , 1) **2 . D0+ 
C SFOCT ( 2 , 2 ) **2 . DO+SFOCT (2,3)**2. DO) 

GCOEF (7 ) -4 . DO* (C**2 . DO-DELTA**2 . DO) 

GCOEF (8) — 4.DO*(ALPHA*DELTA+2.DO*C**2.DO*SFOCT(2, 1) ) 

GCOEF ( 9 ) “-8 . DO*BETA*DELTA 

GCOEF ( 10) —8 . DO*GAMMA*DELTA 

RETURN 

STOP 

END 
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Appendix 2 

Program listing of the Transmit Mode Raytracing error functional optimization 

software 
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************************** 




* Transmit Mode Raytracing Optimization Code 

* 

* James W. LaPean, Jr. 11/18/93 

************************************************************************ 


PROGRAM TMRT1 
IMPLICIT NONE 

REAL* 8 SRFMD ( 4 ) , SRFA (25,25,6) ,SCND(6) , SCAN ( 2 ) , FOC , RMSERR , TRANS ( 5 ) , 
C FEED ( 3 ) , SC (2 ) , PI , XI (5 , 5) ,TOL,SFOC(2,3) 

INTEGER* 2 ICN , I , J , T , P , TINC , PINC , N , ITER 
CHARACTER* 12 SUBFILE 

COMMON /REFL/ SRFMD , SRFA , FOC , FEED , SC 
PI«3 . 14 15926535897932384 6D0 
TOL-l.D-lO 

* Program inputs read from input file 

OPEN ( 3 , FILE- ' TMRT1 . INP ' , STATUS- ' OLD ' ) 

READ ( 3 , * ) SRFMD ( 1 ) , SRFMD ( 2 ) , SRFMD ( 3 ) ,SRFMD(4) 

READ ( 3 , * ) FOC , FEED ( 1 ) , FEED(2) , FEED (3) 

READ ( 3 , * ) SCND ( 1 ) ,SCND(2) ,SCND(3) ,SCND(4) , TINC, PINC 

READ (3 , *) ICN 

READ ( 3 , ' ( A12 ) • ) SUBFILE 

* Nominal subreflector data entry 

OPEN ( 4 , FILE-SUBFILE , STATUS- ' OLD ' ) 

READ ( 4 , ' ( 6F12 . 7 ) ' ) SFOC (1,1), SFOC (1,2), SFOC (1,3), 

C SFOC (2,1) , SFOC (2,2) , SFOC (2,3) 


DO 1 1-1,25 
DO 1 J-1,25 

1 READ (4, ' (6F12.7) ' ) SRFA ( I , J , 1) , SRFA(I , J, 2 ) , SRFA ( I , J , 3 ) , 

C SRFA ( I , J , 4 ) , SRFA (I,J,5) , SRFA ( I , J , 6 ) 


CLOSE (4) 

* Perform scan optimized fitting for scan combinations 
WRITE (*,'(//, Al)')' ' 

SCAN ( 2 ) —SCND ( 2 ) 

P-1 

2 SCAN (1) -SCND (1) 

T— 1 

3 WRITE ( * , ' (A25 , D8 . 3 , A7 , D8 . 3 ) ' ) ' Calculating for Theta - ',SCAN(1), 

C 'Phi - ' , SCAN ( 2 ) 

SC ( 1 ) —SCAN ( 1 ) *PI/ 1 . 8D2 
SC (2) -SCAN (2) *PI/1.8D2 
DO 1-1,5 

TRANS (I) -0. dO 
DO J-1,5 

IF (I.EQ.J) THEN 
XI ( I , J) —1 . DO 
ELSE 

XI (I, JJ-O.DO 
ENDIF 
ENDDO 
ENDDO 


N— 5 

CALL POWELL (TRANS , XI , N , TOL, ITER, RMSERR) 
CALL OUT (FEED, TRANS, SCAN, ICN, SFOC) 

ICN— ICN+1 

SCAN ( 1 ) -SCAN ( 1 ) +SCND ( 3 ) 

T*T+1 

IF (T. LE.TINC) GOTO 3 
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SCAN ( 2 ) =SCAN ( 2 ) +SCND ( 4 ) 
P=P+1 

IF (P.LE.PINC) GOTO 2 

STOP 

END 


* 
* 


Powell's Method Optimization 


SUBROUTINE POWELL (P , XI , N, FTOL, ITER , FRET) 

IMPLICIT NONE 

^REAL*8 p^' XI ( 5 '5) ,PT(5) ,PTT(5) , XIT (5) , FTOL, FPTT, FRET, FP, T, DEL, 

INTEGER*2 I , J, ITER, N, IBIG, ITMAX 

ITMAX-200 

FRET=FUNC(P) 

DO J»1,N 
PT(J)=P(J) 

END DO 
ITER=0 
ITER-ITER+1 
FP=FRET 
IBIG=0 
DEL= 0 • DO 
DO 1=1, N 
DO J=1,N 

XIT(J)=XI(J,I) 

END DO 
FPTT=FRET 

CALL LINMIN(P, XIT, FRET) 

IF (ABS (FPTT-FRET) .GT.DEL) THEN 
DEL=ABS (FPTT-FRET) 

IBIG=I 
END IF 
END DO 

IF (2 . DO*DABS ( FP-FRET) . LE. FTOL* (DABS (FP) +DABS f FRFT -„ DM 

" (ITEP.EQ. ITMAX, PAUSE 'POMel! 

PTT(J)=2.D0*P(J)-PT(J) 

XIT(J) =P(J) -PT ( J) 

PT(J)=P(J) 

END DO 

FPTT=FUNC ( PTT) 

IF (FPTT . GE. FP) GOTO 1 

l? ■ 1 <n-«CT-DH.) **2 . DO-DEL* (FP-FPTT, .*2.00 

CALL LINMIN (P, XIT, FRET) 

DO J=1,N 

XI (J , IBIG) =XIT ( J ) 

END DO 
GOTO 1 
END 


SUBROUTINE LINMIN (P , XI , FRET) 
IMPLICIT NONE 
EXTERNAL F1DIM 
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REAL* 8 P(5),XI(5) , A a , XX , FRET , TOL , PCOM ( 50 ) , XICOM ( 50) , BX , FA, FX , FB . 
C XMIN, BRENT 

INTEGER*2 J,NCOM 

COMMON /F1COM/ PCOM, XICOM, NCOM 

TOL-l.D-lO 

NCOM-5 

DO J«l,5 

PCOM(J)-P(J) 

XICOM (J) —XI ( J) 

END DO 
AX-0. 

XX™1 • 

CALL MNBRAK ( AX , XX , BX , FA , FX , FB , F1DIM) 

FRET-BRENT (AX , XX , BX , F1DIM , TOL , XMIN) 

DO J-1,5 

XI ( J) -XMIN*XI (J) 

P(J)-P(J)+XI(J) 

END DO 
RETURN 
END 


REAL* 8 FUNCTION FIDIM(X) 

IMPLICIT NONE 

REAL* 8 PCOM ( 50 ), XICOM (50) , XT (50) , X, FUNC 
INTEGER*2 NCOM,J 

COMMON /F1COM/ PCOM, XICOM , NCOM 
DO J— 1 , NCOM 

XT ( J) — PCOM ( J) +X*XICOM ( J) 

END DO 

F1DIM— FUNC (XT) 

RETURN 

END 


SUBROUTINE MNBRAK (AX , BX , CX , FA , FB , FC , FUNC) 

IMPLICIT NONE 

REAL*8 AX , BX , CX , FA , FB , FC , DUM , GLIMIT , GOLD , TINY , R , Q , U , ULIM , FU , 

C FUNC 

GOLD-1. 618034D0 
GLIMIT— 100 . DO 
TINY-1 . D-20 
FA- FUNC (AX) 

FB-FUNC(BX) 

IF(FB.GT.FA) THEN 
DUM-AX 
AX-BX 
BX-DUM 
DUM-FB 
FB-FA 
FA- DUM 
END IF 

CX-BX+GOLD* (BX-AX) 

FC— FUNC (CX) 

1 IF (FB.GT.FC) THEN 

R- (BX-AX) *(FB-FC) 

Q-(BX-CX)* (FB-FA) 

U-BX-( (BX-CX)*Q- (BX-AX) *R) / (2 . *SIGN(MAX(ABS(Q-R) ,TINY) ,Q-R) ) 
ULIM-BX+GLIMIT* (CX-BX) 

IF ( (BX-U) *(U-CX) .GT.O. ) THEN 
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FU-FUUC v 

IF (FL.LT.FCj THEN 
AX=BX 
FA=FB 
BX-U 
FB-FU 
RETURN 

ELSE IF (FU.GT.FB) THEN 
CX-U 
FC-FU 
RETURN 
END IF 

U=CX+GOLD* (CX-BX) 

FU=FUNC(U) 

ELSE IF ( (CX-U) * (U-ULIM) . GT. 0 . ) THEN 
U-ULIM 
FU-FUNC(U) 

ELSE 

U=CX+GOLD* (CX-BX) 

FU-FUNC(U) 

END IF 
AX-BX 
BX-CX 
CX=U 
FA-FB 
FB-FC 
FC—FU 
GOTO 1 
END IF 
RETURN 
END 


REAL*8 FUNCTION BRENT (AX, BX, CX , F, TOL, XMIN) 

IMPLICIT NONE 

REAL*8 AX, BX, CX, F , TOL, XMIN , CGOLD, ZEPS, A, B, V, W, X, E, FX, FV, FW, XM, 
C TOLl,TOL2,R,Q,ETEMP,P,D,U,FU 

INTEGER*2 ITER, ITMAX 
ITMAX-100 
CGOLD 22 * 3819660D0 
ZEPS—1 , 0D-10 
A-MIN(AX,CX) 

B-MAX(AX,CX) 

V-BX 

w=v 

x=v 

E=0 . 

FX-F(X) 

FV-FX 

FW-FX 

DO ITER-1, ITMAX 
XM— 0 . 5* ( A+B) 

TOL1— TOL*ABS (X) +ZEPS 
TOL2-2 . *TOLl 

IF (ABS(X-XM) *LE. (TOL2-. 5* (B-A) ) ) GOTO 3 
IF (ABS(E) .GT.TOL1) THEN 
R— (X-W) * (FX-FV) 

Q* (X-V) * ( FX-FW) 

P-(X-V) *Q- (X-W) *R 
Q~2.*(Q-R) 
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IF (Q.GT.O) P--P 
Q-ABS(Q) 

ETEMP-E 

E*D _ 

IF (ABS (P) . GE. ABS ( . 5*Q*ETEMP) .OR.P.LE.Q* (A-X) .OR. 
P.GE.Q* (B-X) ) GOTO 1 

D-P/Q 

IF^U-A . LT . TOL2 . OR . B-U . LT . TOL2 ) D=SIGN (TOL1 , XM-X) 
GOTO 2 
END IF 

1 IF (X.GE.XM) THEN 

E-A-X 
ELSE 
E-B-X 
END IF 
D-CGOLD*E 

2 IF (ABS (D) .GE.TOL1) THEN 

U-X+D 

ELSE 

U— X+S IGN ( TOL1 , D ) 

END IF 
FU-F(U) 

IF (FU.LE.FX) THEN 
IF (U. GE * X) THEN 
A-X 
ELSE 
B-X 

END IF 
V-W 
FV-FW 
W=X 
FW=FX 
X-U 
FX-FU 
ELSE 

IF (U.LT.X) THEN 
A-U 
ELSE 
B-U 

END IF 

IF (FU.LE.FW.OR.W.EQ.X) THEN 
V-W 
FV=FW 

w-u 

FW-FU 

ELSE IF ( FU . LE . FV . OR . V . EQ . X . OR . V * EQ . W) THEN 
V-U 
FV-FU 
END IF 
END IF 
END DO 

PAUSE 'Brent exceed maximum iterations.' 

3 XMIN-X 
BRENT- FX 
RETURN 
END 


************************************************************************ 
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* Powell's Method Error Functional 

* 

************************************************************************ 

REALMS FUNCTION FUNC (TRANS) 

IMPLICIT NONE 

REAL* 8 SRFAT (25,25,6) , TRANS (5) , RMSERR 
CALL POSIT (SRFAT, TRANS) 

CALL XMTTRC (SRFAT, RMSERR) 

FUNC-RMS ERR 

RETURN 

STOP 

END 


************************************ *********************************** 


Type 1 Concept Subreflector Positioning Code 

************************************************* ********************** 


SUBROUTINE POSIT (SRFAT, TRANS) 

IMPLICIT NONE 

REAL* 8 SRFA (25,25,6) , SRFAT (25 , 25 , 6) ,TRANS(5) ,SRFMD(4) ,FOC,FEED(3) , 
C SCAN (2 ) , T ( 5) 

INTEGER*2 I,J 

COMMON /REFL/ SRFMD, SRFA, FOC, FEED, SCAN 
T(1)«TRANS(3) 

T ( 2 ) =TRANS ( 4 ) 

T(3)«TRANS(5) 

T(4) -TRANS (1) 

T(5) -TRANS (2 ) 

DO 1 1=1,25 
DO 1 J=1 , 2 5 

SRFAT ( I, J, 1 ) =T ( 1 ) +SRFA (13,13,1) - 

SIN(T(4))*( SRFA ( I , J, 2 ) -SRFA (13,13,2) )+ 

COS (T ( 4 ) ) *SIN (T (5) ) * ( SRFA ( I , J, 3 ) -SRFA ( 13 , 13 , 3 ) ) + 

COS (T ( 4 ) ) *COS (T ( 5) ) * (SRFA (I, J, 1) -SRFA ( 13 , 13 , 1 ) ) 

SRFAT ( I , J , 2 ) =T ( 2 ) +SRFA (13,13,2)+ 

COS (T(4))*( SRFA ( I , J, 2 ) -SRFA (13,13,2) )+ 

SIN(T(4) ) *SIN (T (5 ) ) * (SRFA ( I , J, 3 ) -SRFA( 13 , 13 , 3 ) ) + 

SIN (T ( 4 ) ) *COS (T (5) ) * ( SRFA ( I , J, 1) -SRFA ( 13 , 13 , 1 ) ) 

SRFAT ( I , J , 3 ) — T ( 3 ) +SRFA (13,13,3)+ 

COS (T ( 5) ) * (SRFA ( I, J, 3) -SRFA (13, 13, 3) )- 
SIN(T(5) ) *( SRFA ( I , J, 1 ) -SRFA (13,13,1)) 

SRFAT ( I , J, 4 ) —-SIN (T ( 4 ) ) *SRFA ( I , J , 5) + 

COS (T (4 ) ) *SIN(T(5) ) *SRFA ( I , J , 6 ) + 

COS (T (4 ) ) *COS (T ( 5) ) *SRFA ( I , J , 4 ) 

SRFAT (I , J , 5) —COS (T(4) ) *SRFA ( I , J , 5) + 

SIN(T(4))*SIN(T(5) ) *SRFA ( I , J , 6) + 

SIN (T (4 ) ) *COS (T ( 5 ) ) *SRFA(I, J,4) 

1 SRFAT (I , J, 6) -COS (T(5) ) *SRFA ( I , J, 6) - 

C SIN(T(5) ) *SRFA ( I , J , 4 ) 

RETURN 

STOP 

END 


********************************** ************************************** 
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* Type 1 Concept Transmit Mode Raytracing Error Functional 

* 

************************************************************************ 

SUBROUTINE XMTTRC (SRFAT, RMSERR) 

IMPLICIT NONE 

REAL* 8 SRFM ( 6 ) , SRFAT ( 25 , 25 , 6) , FEED(3) ,R1(6) ,R2(6) ,R3(6) ,SCAN<2) , 

C RDOTN , A, B , C, FOC, R2L, PI , RMSERR, SRFMD (4 ) , L, ANGLE ( 2 5 , 25) , 

C DEN , RA ( 3 ) , SRFA (25,25,6) 

INTEGER*2 I , J , MISS ( 25 , 25) , K 

COMMON /REFL/ SRFMD, SRFA, FOC, FEED, SCAN 

PI-3. 14 1592653 58979323846D0 

* FIND IDEAL COMPONENTS OF REFLECTED RAY FROM MAIN REF. TO APERTURE PLANE 

RA ( 1 ) —COS ( SCAN ( 1 ) ) 

RA ( 2 ) —SIN ( SCAN ( 1 ) ) *COS (SCAN(2) ) 

RA ( 3 ) —SIN ( SCAN ( 1 ) ) *SIN ( SCAN ( 2 ) ) 

DO 1 1-1,25 
DO 1 J-1,25 
MISS (I , J) -1 

* FIND RAY AND RAY COMPONENTS FROM FEED TO SUBREFLECTOR 

R3 ( 1 ) -SRFAT ( I , J , 1 ) -FEED(l) 

R3 (2 ) -SRFAT (I , J, 2 ) -FEED ( 2 ) 

R3 ( 3 ) —SRFAT ( I , J , 3 ) -FEED ( 3 ) 

L— DSQRT(R3 (1) **2 . D0+R3 ( 2 ) **2 . D0+R3 (3) **2.D0) 

R3(4)=R3(1)/L 
R3(5)-R3(2)/L 
R3 (6)-R3 (3) /L 

* FIND COMPONENTS OF RAY REFLECTED FROM SUBREFLECTOR 

RDOTN-R3 ( 4 ) * SRFAT ( I , J, 4) +R3 ( 5 ) * SRFAT ( I , J , 5 ) +R3 ( 6 ) *SRFAT ( I , J , 6 ) 
R2 (4 ) — R3 ( 4 ) -2 . *RDOTN*SRFAT (I, J,4) 

R2 (5) -R3 ( 5 ) -2 . *RDOTN*SRFAT (I, J, 5) 

R2 (6) -R3 (6) -2 . *RDOTN*SRFAT (I , J, 6) 

* FIND INTERSECTION OF RAY FROM SUBREFLECTOR WITH MAIN REFLECTOR 

A—R2 ( 5 ) **2 . D0+R2 (6) **2.D0 

B— 2 . D0*R2 ( 5 ) * SRFAT ( I , J , 2 ) +2 . D0*R2 ( 6 ) *SRFAT (I , J, 3 ) 

C -4 . D0*R2 ( 4 ) *FOC 

C— SRFAT ( I , J , 2 ) * * 2 . DO+SRFAT ( I , J , 3 ) **2 . DO-4 . DO *SRFAT (I, J, 1) *FOC 

R2L-(DSQRT(B**2 .DO-4 . D0*A*C) -B) / (2 .DO*A) 

SRFM (1) —SRFAT ( I , J , 1) +R2 (4 ) *R2L 
SRFM ( 2 ) —SRFAT ( I , J , 2 ) +R2 ( 5 ) *R2L 
SRFM ( 3 ) —SRFAT (I , J, 3 ) +R2 (6) *R2L 

* FIND NORMAL OF MAIN REFLECTOR AT RAY INTERSECTION POINT 

DEN-DSQRT ( (SRFM ( 2 ) / ( 2 . DO* FOC) ) **2 . D0+ 

C ( SRFM ( 3 ) / ( 2 . DO*FOC) ) **2.D0+1.D0) 

SRFM (4 )-l. DO/DEN 

SRFM ( 5 ) —-SRFM ( 2 ) / ( 2 . DO* FOC* DEN) 

SRFM (6) — SRFM ( 3 ) / ( 2 . DO*FOC*DEN) 

* INDICATE A MISS IF RAY HITS AN UNUSED PART OF THE PARABOLA 

IF ( SRFM ( 2 ) . GE . SRFMD ( 1 ) - ( SRFMD ( 2 ) -SRFMD (1))/1.D2) THEN 
IF ( SRFM ( 2 ) . LE . SRFMD ( 2 ) + ( SRFMD ( 2 ) -SRFMD (1) ) / 1 . D2 ) THEN 
IF ( SRFM ( 3 ) . GE . SRFMD ( 3 ) - ( SRFMD ( 4 ) -SRFMD (3) ) / 1 . D2 ) THEN 
IF ( SRFM ( 3 ) . LE . SRFMD ( 4 ) + ( SRFMD ( 4 ) -SRFMD ( 3 ) ) / 1 . D2 ) THEN 
MISS (I , J) -0 
ENDIF 
ENDIF 
ENDIF 
ENDIF 

* FIND COMPONENTS OF RAY REFLECTED FROM MAIN REFLECTOR 

RDOTN-R2 ( 4 ) * SRFM ( 4 ) +R2 ( 5 ) * SRFM ( 5 ) +R2 ( 6 ) *SRFM ( 6 ) 
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R1 (4) =R2 (4) -2 . *RDOTN*SRFM (4 ) 

R1 (5) =R2 ( 5) -2 . *RDOTN*SRFM (5) 

R1 (6) =R2 ( 6 ) -2 . *RDOTN*SRFM ( 6 ) 

* FIND CROSS PRODUCT OF ACTUAL AND IDEAL REFLECTED RAYS TO THE APERTURE 

1 ANGLE ( I , J) = (R1 (5) *RA ( 3 ) -RA ( 2 ) *R1 { 6) ) **2.D0+ 

C (R1 ( 4 ) *RA ( 3 ) -RA (1)*R1(6))**2. D0+ 

C (R1 (4) *RA(2) -RA(1) *R1 (5) ) **2.D0 

* ERROR ANALYSIS 

K“0 

RMSERR=0 . DO 
DO 2 1=1,25 
DO 2 J=1 , 2 5 

RMS ERR=RMSERR+ ANGLE ( I , J) 

2 K=K+MISS ( I , J) 

RMS ERR=DSQRT (RMS ERR) 

RETURN 

STOP 

END 


********************************** ************************************** 
* 

* Type 1 Concept Output Code 

* 

************************************************************************ 

SUBROUTINE OUT ( FEED, TRANS , SCAN , ICN,SFOC) 

IMPLICIT NONE 

REAL*8 FEED (3) , TRANS (5) , SCAN (2) , SRFAT (25 , 2 5 , 6) , PI , 

C SFOC(2,3) ,GCOEF(10) ,U,V 

INTEGER*2 ICN,K(3) 

CHARACTER* 12 OUTFILE 

PI=3 . 14 1592653 5897932 384 6D0 

U=SIN ( SCAN ( 1 ) *PI / 1 . 8D2 ) *COS(SCAN(2) *PI/1.8D2) 

V=SIN ( SCAN ( 1 ) *PI / 1 . 8D2 ) *SIN ( SCAN ( 2 ) *PI / 1 . 8D2) 

K(1)=48+ICN/100 

K (2 ) =48+ ( ICN- (K(l)-48) *100) /10 
K ( 3 ) “48+ICN- (K(l) —48) *100- (K ( 2 ) -48 ) *10 
CALL POSITF(SFOC, TRANS, GCOEF) 

CALL POSIT (SRFAT, TRANS) 

OUTFILE= ' TMRT1XXX . MOV ' 

OUTFILE (6: 6) “CHAR (K ( 1 ) ) 

OUTFILE (7 : 7 ) “CHAR (K ( 2 ) ) 

OUTFILE (8:8) “CHAR (K ( 3 ) ) 

OPEN (4 , FILE“OUTFILE, STATUS® ' UNKNOWN ' ) 

WRITE(4, ' (A10,2F12. 7) ') 'SCANNED TO ' , SCAN ( 1 ) , SCAN ( 2 , 

WRITE (4, ' (5F15. 10) ') TRANS (4) , TRANS (5) ,TRANS(3) ,TRANS(1) , TRANS (2) 
ENDFILE(4 ) 

CLOSE (4) 

* 20 GHz GRASP7 Input File 

OUTFILE® 'TMRT1XXX. INP' 

OUTFILE (6:6) “CHAR (K ( 1 ) ) 

OUTFILE (7:7) “CHAR (K ( 2 ) ) 

OUTFILE (8:8) “CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE=OUTFILE , STATUS® ' UNKNOWN ' ) 

WRITE (4, ' (A30) ') ****** GRASP77 INPUT FILE ****** 

WRITE(4, ' (A10,2F12. 7) ') 'SCANNED TO ' , SCAN ( 1 ) , SCAN ( 2 ) 

WRITE (4, ' (F12.7,4I2) ' ) 1 . 5D-2 , 2 , 1 , 0 , 0 

WRITE (4, ' (A31) ') '** MAIN REFLECTOR INPUT DATA **' 

WRITE (4, ' (3F12 .7) ' ) 0. , 0. ,0. 
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WRITE (4 , ' (3F12.7) ')1. ,0. ,0. 

WRITE ( 4 , ' (3F12.7) ')0. ,1. ,0. 

WRITE (4 , ' ( 12 , 2F12 . 7 ) ') 1,7.795,0. 

WRITE ( 4 , ' (312) '54,0,0 

WRITE (4, ' (4F12.7) ' ) 0 . , 0 . , 0 . , 13 . 5 

WRITE (4, ' (2F12.7) ' ) 5 . 3 15 , 5 . 3 15 

WRITE (4, ' (A30) ') '** SUB REFLECTOR INPUT DATA **' 

WRITE (4, ' (3F12.7) ')0. ,0. ,0. 

WRITE (4, ' (3F12.7) ')1. ,0. ,0. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' ( 12 , 2F12 . 7 ) ' ) 1, (SRFAT(25, 13 , 2) +SRFAT(1 , 13 , 2) ) /2.D0, 

C ( SRFAT (13,25,3) +SRFAT (13,1,3) ) /2.D0 

WRITE ( 4 , ' (212) ' ) 3 , 0 

WRITE (4, ' (5F14.7) ')GCOEF(l) ,GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 
WRITE (4, ' (5F14.7) ')GC0EF(6) ,GCOEF(7) ,GCOEF(8) ,GC0EF(9) ,GCOEF(10) 
WRITE (4, ' ( 5F12 .7,15) ')0.,0.,0.,0.,0.,1 

WRITE (4, ' (2F12.7) ') (SRFAT ( 25 , 13 , 2 ) -SRFAT ( 1 , 13 , 2 ) ) /2.D0, 

C (SRFAT (13, 25, 3) -SRFAT (13, 1,3) ) /2.D0 

WRITE (4 , 9 (A21) ') '** FEED INPUT DATA **' 

WRITE (4, ' (3F12.7) ' ) FEED ( 2 ) , FEED ( 3 ) , FEED ( 1) 

WRITE (4, # (3F12.7) ') . 964483694 , 0 . 264142395 
WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE ( 4 , ' (12) ') 1 

WRITE (4, ' (6F12.7) ')0.,0.,0.,0.,0.,0. 

WRITE (4, ' (F12.7,I2) ')0. ,1 

WRITE (4, ' (4F12.7, 312) ')0. ,0. , 1. ,90. ,0,0,3 
WRITE ( 4 , ' (12) ' ) 6 

WRITE (4, ' ( 5F12 .7,13) ' ) -15 . , -15 . , 7 . 1574 , 0 . ,0. ,1 

WRITE(4, ' (A40) ' ) '** MAIN REFLECTOR FIELD SPECIFICATION **' 

WRITE(4, ' (12) ') 1 

WRITE (4, ' (3F12.7) ') 7. 795,0. , 10. 

WRITE (4, ' (3F12.7) ' ) 1. ,0. ,0. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' ( 216 , F12 .7,213) ' ) 60, 144 , 0 . , 0 , 0 
WRITE (4 , 9 (313) ') 0,-1, 1 
WRITE ( 4 , 9 (F12.7) ') 10. 

WRITE ( 4 , ' (2F12.7) ')U,V 
WRITE ( 4 , ' (213) ' )0, 1 
WRITE ( 4 , ' (213) 9 ) 3, 1 

WRITE (4, ' (4F12 .7,216 ) 9 ) -1 . D-2 , -1 . D-2 , 1 . D-2 , 1 . D-2 , 2 5 , 25 
WRITE ( 4 , 9 (416) ')3,1 
OUTFILE- ' TMRT1XXX . PI 9 
OUTFILE (6:6) —CHAR ( K ( 1 ) ) 

OUTFILE (7 : 7 ) “CHAR (K ( 2 ) ) 

OUTFILE (8:8) -CHAR { K ( 3 ) ) 

WRITE (4, ' (A12) ') OUTFILE 

WRITE(4 , 9 (A39) ') SUB REFLECTOR FIELD SPECIFICATION 
WRITE (4 , 9 (12) ')0 

WRITE (4 , 9 (A30 ) ' ) 9 ** FEED FIELD SPECIFICATION **’ 

WRITE ( 4 , ' (12) ')0 
ENDFILE ( 4 ) 

CLOSE (4) 

OUTFILE* ' TMRUVXXX . INP ' 

OUTFILE (6:6) «CHAR ( K ( 1 ) ) 

OUTFILE (7:7) “CHAR ( K ( 2 ) ) 

OUTFILE (8 : 8 ) “CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE=OUTFILE , STATUS* ' UNKNOWN 9 ) 

WRITE (4 , ' (A25) 9 ) 9 UVPROC CONTROL INPUT FILE' 

OUTFILE- ' TMRT1XXX . PI ' 

OUTFILE (6:6) -CHAR (K ( 1 ) ) 
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OUTFILE (7:7) =CHAR ( K ( 2 ) ) 

OUTFILE (8:8) =CHAR (K ( J ) ) 

WRITE (4 , ' (A12) ') OUTFILE 
WRITE ( 4 , ' (612) ') 1,1, 1,1, 0,0 
WRITE ( 4 , ' (212) ')1,1 
WRITE (4, ' (612) ') 0,1, 0,0, 1,0 
WRITE (4, ' (A7) ' ) 'TEMP. P2 ' 

WRITE ( 4 , ' (12) ' ) 3 
WRITE (4, '(214)') 101, 101 
WRITE (4, ' (712) ')0, 0,1, 1,0, 0,0 
WRITE ( 4 , ' (412) ') 3, 0,2,0 
WRITE ( 4 , ' (F5. 1) ') 3. 

WRITE ( 4 , ' (F5. 1) ' ) 10. 

WRITE ( 4 , ' (512) ') 1,0, 0, 1, 1 
WRITE (4, ' (F7.4) ') 0.1213 
ENDFILE ( 4 ) 

CLOSE ( 4 ) 

* 40 GHZ GRASP7 Input File 

OUTFILE= ' TMRT1XXX . INP ' 

OUTFILE (6:6) =CHAR (K ( 1 ) +1 ) 

OUTFILE (7:7) =CHAR (K ( 2 ) ) 

OUTFILE (8:8) =CHAR (K ( 3 ) ) 

OPEN (4 , FILE=OUTFILE, STATUS= ' UNKNOWN ' ) 

WRITE (4, ' (A30) ') ****** GRASP77 INPUT FILE ****** 

WRITE(4, ' (A10,2F12. 7) ') 'SCANNED TO ' ,SCAN(1) ,SCAN(2) 

WRITE (4, ' (F12.7,4I2) ' ) 7 . 5D-3 , 2 , 1 , 0 , 0 

WRITE (4, ' (A31) ') '** MAIN REFLECTOR INPUT DATA **' 

WRITE (4, ' (3F12.7) ')0. ,0. ,0. 

WRITE (4 , ' (3F12.7) ' ) 1 . ,0.,0. 

WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE (4, ' ( 12 , 2F12 . 7 ) ') 1,7. 795,0. 

WRITE ( 4 , ' (312) ' ) 4 , 0 , 0 

WRITE ( 4 , ' (4F12.7) ')0. ,0. ,0. , 13.5 

WRITE (4, ' (2F12.7) ')5. 315, 5. 315 

WRITE(4, ' (A30) ') '** SUB REFLECTOR INPUT DATA **' 

WRITE (4, ' (3F12.7) ')0. ,0. ,0. 

WRITE (4, ' ( 3F12 .7)')1.,0.,0. 

WRITE (4, ' (3F12.7) ')0. ,1.,0. 

WRITE ( 4 , ' ( 12 , 2 F12 . 7 ) ')1, (SRFAT ( 2 5 , 13 , 2) +SRFAT ( 1 , 13 , 2 ) ) / 2 . DO , 

C ( SRFAT (13,25,3) +SRFAT (13,1,3) ) /2.D0 

WRITE(4, ' (212) ') 3,0 

WRITE (4 , ' (5F14.7) ')GCOEF(l) ,GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 
WRITE ( 4 , ' (5F14.7) ')GCOEF(6) ,GCOEF(7) ,GCOEF(8) ,GCOEF(9) ,GCOEF(10) 
WRITE (4, ' ( 5F12 . 7 , 15 ) ')0.,0.,0.,0.,0.,1 

WRITE ( 4 , ' (2F12.7) ') ( SRFAT ( 2 5 , 1 3 , 2 ) -SRFAT( 1 , 13 , 2 ) ) /2 . DO , 

C ( SRFAT (13,25,3) -SRFAT (13,1,3))/2.D0 

WRITE ( 4 , ' ( A2 1 ) ' ) ' ** FEED INPUT DATA **' 

WRITE ( 4 , ' (3F12.7) ' ) FEED ( 2 ) , FEED ( 3 ) , FEED( 1) 

WRITE (4 , ' (3F12.7) ') .96448jfa94,0. ,-.264142395 
WRITE (4, ' (3F12.7) ')0. ,1. ,0. 

WRITE ( 4 , ' (12) ') 1 

WRITE (4, ' (6F12.7) ')0.,0.,0.,0.,0.,0. 

WRITE ( 4 , ' (F12. 7, 12) ' ) 0. , 1 

WRITE (4, ' (4F12.7,3I2) ')0. ,0. ,1. ,90. ,0,0,3 
WRITE ( 4 , ' (12) ' )6 

WRITE (4, ' ( 5F12 . 7 , 1 3 ) ' ) -15 . , -1 5 . , 7 . 1574 , 0 . ,0.,1 

WRITE (4 , ' (A40) ') '** MAIN REFLECTOR FIELD SPECIFICATION **' 

WRITE ( 4 , ' (12) ' ) 1 

WRITE (4, ' (3F12.7) ') 7. 795,0. ,10. 

WRITE (4, ' (3F12.7) ' ) 1. ,0. , 0. 
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WRITE ( 4 , 
W "TE ( 4 , 
W TE ( 4 , 
WRITE (4, 
WRITE (4, 
WRITE (4 , 
WRITE (4, 
WRITE (4, 
WRITE (4, 
OUTFILE- 


' (3F12.7) ' )u. , 1 
' (216 , F12 . 1 , 2 '. 

' (313) ' ) 0,-1, 1 
' (F12.7) ' ) 10. 

' (2FI2.7) ')U,V 
' (213 ) ' ) 0 , 1 
' (2X3)')3.1 
' (4F12 . 7 , 216) ' ) 
' (416) ')3,1 
'TMRT1XXX.P1' 


. , 0 . 

) ' ) 60,144,0. ,0,0 


-5.D-3,-5.D-3,5.D-3,5.D-3,25,25 


OUTFILE (6:6) “CHAR (K ( 1) +1) 

OUTFILE (7:7) -CHAR ( K ( 2 ) ) 

OUTFILE (8:8) — CHAR (K ( 3 ) ) 

WRITE (4, ' (A12) ') OUTFILE 

WRITE (4 , ' (A39 ) ') ^ ** SUB REFLECTOR FIELD SPECIFICATION 
WRITE (4 | ' (12) ' ) 0 

WRITE (4 , ' (A3 0 ) ' ) ' ** FEED FIELD SPECIFICATION **' 

WRITE ( 4 , ' (12) ' )0 
ENDFILE(4) 

CLOSE (4) 

OUTFILE- ' TMRUVXXX . INP ' 

OUTFILE (6: 6) -CHAR (K ( I) +1) 

OUTFILE (7:7) -CHAR ( K ( 2 ) ) 

OUTFILE (8 : 8) -CHAR(K ( 3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

WRITE (4, ' (A25) ') 'UVPROC CONTROL INPUT FILE' 

OUTFILE- ' TMRT1XXX . PI ' 


OUTFILE (6:6) =CHAR(K ( 1 ) +1 ) 
OUTFILE (7 : 7) -CHAR (K ( 2 ) ) 
OUTFILE ( 8 : 8 ) -CHAR (K ( 3 ) ) 

WRITE ( 4 , ' (A12) ') OUTFILE 
WRITE ( 4 , ' (612) ') 1,1, 1,1, 0,0 
WRITE ( 4 , ' (212) ')1,1 
WRITE (4, ' (612) ') 0,1, 0,0, 1,0 
WRITE (4, ' (A7) ') 'TEMP . P2 ' 

WRITE ( 4 , ' (12) ' ) 3 
WRITE ( 4 , ' (214) ') 101, 101 
WRITE ( 4 , ' (712) ')0, 0,1, 1,0, 0,0 
WRITE ( 4 , ' (412) ')3, 0,2,0 
WRITE (4 , ' (F5.1)')3. 

WRITE (4 , ' (FS • 1) ' ) 10. 

WRITE (4,' (SI2) ') 1,0, 0,1,1 
WRITE (4, ' (F7.4) ') 0.0606 
ENDFILE(4) 

CLOSE (4) 

* 80 GHz GRASP7 Input File 


OUTFILE- 'TMRT1XXX. INP' 

OUTFILE (6:6) -CHAR ( K ( 1 ) +2 ) 

OUTFILE (7:7) — CHAR (K ( 2 ) ) 

OUTFILE (8:8) — CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE-OUTFILE , STATUS- ' UNKNOWN ' ) 

VPITE(4 , ' (A30) ') '***** GRASP77 INPUT FILE *****' 

ITE(4 , ' (A10,2F12.7) ') ' SCANNED TO ' , SCAN ( 1) , SCAN ( 2 ) 
WRITE (4 , ' (F12 .7,412) ' ) 3 . 75D-3 , 2 , 1 , 0, 0 
WRITE (4 , ' (A31) ' ) ' ** MAIN REFLECTOR INPUT DATA **' 


WRITE (4, ' (3F12.7) ')0. ,0. ,0. 

WRITE (4 , ' (3F12.7) ') 1. ,0. ,0. 

WRITE ( 4 , ' (3F12.7) ')0. ,1. ,0. 

WRITE ( 4 , ' ( 12 , 2F12 . 7 ) ') 1,7. 795,0. 
WRITE ( 4 , ' (312) ' ) 4 , 0 , 0 


**' 
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WRITE (4, 
WRITE ( 4 , 
WRITE (4, 
WRITE (4, 
WRITE (4, 
WRITE (4 , 
WRITE (4, 
C 

WRITE (4 , ' 
WRITE (4, 1 
WRITE ( 4 , ‘ 
WRITE ( 4 , ‘ 
WRITE (4 , • 
C 

WRITE (4 , ' 
WRITE ( 4 , ' 
WRITE ( 4 , ' 
WRITE (4, ' 
WRITE (4, ' 
WRITE (4, ' 
WRITE ( 4 , ' 
WRITE ( 4 , ' 
WRITE (4, ' 
WRITE ( 4 , ' 
WRITE (4, ' 
WRITE ( 4 , ' 
WRITE (4, ' 
WRITE ( 4 , ' 
WRITE ( 4 , ' 
WRITE ( 4 , ' i 
WRITE ( 4 , ' | 
WRITE (4 , ' | 
WRITE ( 4 , ' ( 
WRITE ( 4 , ' ( 
WRITE (4, ' ( 
WRITE ( 4 , ' ( 


' . ‘ - * - • ■ ) ' ) o . , 0 . , 0 . , 1 3 . 5 
' 4 2F12 . ?) ') 5. 3 1L, 5. 315 

' (A30) ' ) ' ** SUB REFLECTOR INPUT DATA **' 

' (3F12.7) ')0. ,0. ,0. 

' (3F12.7) ' ) 1 . ,0. ,0. 

' ( 3F12 . 7 ) ' ) 0 . , 1 . , 0 . 

' ( 12 , 2F12 .7) ' ) 1 , ( SRFAT (25, 13,2) +SRFAT ( 1 , 1 3 , 2 ) ) / 2 . DO , 
n (SRFAT (13,25,3) +SRFAT (13,l,3))/2. DO 

' (5F14.7) ')GCOEF(l) ,GCOEF(2) ,GCOEF(3) ,GCOEF(4) ,GCOEF(5) 

' (5F12 7 I5^) > 0 F ^0^ 'S C °0 EF( 0 ) ,GCOEF(8) » gco EF"( 9) ,GC0EF(10) 
' (2F12.7) ') (SRFAT(25, 13,2) -SRFATfl, 13,2) ) /2.D0, 

, (SRFAT (13, 2 5, 3) -SRFAT (13, 1,3) )/2.D0 

(A21) ' ) '** FEED INPUT DATA **' 

' (3F12 . 7 ) ' ) FEED ( 2 ) , FEED ( 3 ) , FEED ( 1) 

' (3F12.7) ') .964483694,0. ,-.264142395 
’ (3F12.7) ') 0. , 1. ,0. 

’ (12) ')1 

(6F12.7) ')0. ,0. ,0. ,0. ,0. ,0. 

(F12.7,I2) ' ) 0 . , 1 

(4F12.7,3I2) ')0.,0.,1. ,90. ,0,0,3 
(I2)')6 

(5F12. 7, 13) ' ) -15. , -15. ,7.1574,0. ,0. , 1 

(A40) ' ) ' ** MAIN REFLECTOR FIELD SPECIFICATION **' 

(12) ') 1 

(3F12.7) ')7. 795,0. ,10. 

(3F12.7) ')1. ,0. ,0. 

(3F12.7) ')0. ,1. ,0. 

(216, F12. 7, 213) ') 60, 144,0. ,0,0 


, 0 . , 0 . 


WRITE (4 , ' (313) ')0,-l, 1 
WRITE ( 4 , ' (F12.7) ') 10. 
WRITE ( 4 , ' (2F12.7) ')U,V 
WRITE ( 4 , ' (213) ' ) 0, 1 
WRITE ( 4 , ' (213) ') 3, 1 
WRITE ( 4 , ' ( 4 F12 . 7 , 216) ') 
WRITE ( 4 , ' (416) ' ) 3, 1 
OUTFILE= ' TMRT1XXX . PI ' 


• 5D-3,-2.5D-3,2 


•3,2. 5D-3 ,25,25 


OUTFILE (6:6) =CHAR (K ( 1 ) +2 ) 

OUTFILE (7:7) =CHAR (K ( 2 ) ) 

OUTFILE (8:8) =CHAR(K ( 3 ) ) 

WRITE ( 4 , ' (A12 ) ' ) OUTFILE 

WRITE (4 ' (12^' ) 0 J ** SUB REFLECT0R FIELD SPECIFICATION **' 

WRITE(4' ' ( I2°' 'o'** FEED FIEL ° SPECIFICATION **' 

ENDFILE(4) 

CLOSE (4) 

OUTFILE*= 'TMRUVXXX. INP' 

OUTFILE (6:6) =CHAR (K ( 1 ) +2 ) 

OUTFILE (7:7) =CHAR (K (2) ) 

OUTFILE (8:8) =CHAR (K ( 3 ) ) 

OPEN ( 4 , FILE=OUTFILE, STATUS= ' UNKNOWN ' ) 

WRITE ( 4 , ' ( A2 5 ) ' ) ' UVPROC CONTROL INPUT FILE' 

OUTFILE** ' TMRT1XXX . PI ' 

OUTFILE ( 6 : 6) =CHAR (K ( 1 ) +2 ) 

OUTFILE (7:7) =CHAR (K ( 2 ) ) 

OUTFILE (8:8) =CHAR (K ( 3 ) ) 

WRITE ( 4 , ' (A12 ) • ) OUTFILE 
WRITE(4,'(6I2)')1, 1,1, 1,0,0 
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WRITE (4 , ' (212) ')1,1 

WRITE (4, ' (612) ') 0, 1,0,0, 1,0 

WRITE ( 4 , ' (A7 ) ' ) 'TEMP . P2 ' 

WRITE ( 4, ' (12) ' ) 3 
WRITE ( 4 , ' (214) ') 101,101 
WRITE (4 , ' (712) ')0, 0,1, 1,0, 0,0 
WRITE(4,' (412) ' ) 3 , 0,2,0 
WRITE ( 4 , ' (F5.1) ' ) 3 . 

WRITE ( 4 , ' (F5.1) ' ) 10, 

WRITE (4 , ' (512) ')1, 0,0, 1,1 
WRITE ( 4 , ' (F7.4) ') 0.0303 
ENDFILE ( 4 ) 

CLOSE (4) 

RETURN 

STOP 

END 


********* 


*********************************** 


**************************** 


Type 1 Concept Transmit Mode Raytracing Subreflector Focal Point 
Positioning Code 


************************************* 


********************************* 


SUBROUTINE POSITF ( SFOC , TRANS , GCOEF) 
IMPLICIT NONE 

REAL* 8 SRFA (25,25,6) ,SFOC(2,3) ,TRANS(5) , 
C SCAN (2) , T ( 5 ) , GCOEF (10) ,SFOCT(2,3) 

C DELTA, C 

INTEGER* 2 I 

COMMON /REFL/ SRFMD , SRFA, FOC , FEED, SCAN 
T(1)«*TRANS(3) 

T(2)-TRANS(4) 

T(3)«TRANS(5) 

T ( 4 ) “TRANS ( 1 ) 

T ( 5 ) “TRANS ( 2 ) 

DO 1 1-1,2 

SFOCT ( I , 1 ) -T ( 1 ) +SRFA (13,13,1)- 
C SIN(T(4) )*( SFOC (1,2) 

C COS (T(4) )*SIN(T(5) )*( SFOC (1,3) 

C COS (T (4 ) ) *C0S (T (5) ) * (SFOC (1, 1) 

SFOCT (1,2) -T ( 2 ) +SRFA (13, 13,2)+ 

C COS (T ( 4 ) ) * (SFOC (1 , 2) 

C SIN{T(4) ) *SIN(T(5) ) * (SFOC(I, 3) 

C SIN(T(4) ) *COS (T ( 5 ) ) * (SFOC(I , 1) 

SFOCT (1,3) — T(3 ) +SRFA (13,13,3)+ 

C COS (T(5) ) * (SF0C(I , 3) 

C SIN(T(5) ) * (SFOC (1,1) 


SRFMD ( 4 ) , FOC , FEED ( 3 ) , 
.ALPHA, BETA, GAMMA, 


-SRFA(13, 13,2) )+ 
-SRFA(13, 13,3) )+ 
-SRFA (13, 13,1) ) 

-SRFA (13,13,2) )+ 
-SRFA (13,13,3))+ 
-SRFA(13, 13, 1) ) 

-SRFA(13 ,13,3))- 
-SRFA(13 ,13,1)) 


C— DSQRT ( ( SRFA (13,13,1) -SFOC (2,1))**2. D0+ 

C (SRFA (13,13,2) -SFOC (2,2)) **2 . D0+ 

C ( SRFA (13,13,3) -SFOC (2,3) )**2. DO) - 

C DSQRT ( ( SRFA (13,13,1) -SFOC ( 1 , 1 ) ) **2 . D0+ 

C (SRFA (13, 13, 2) -SFOC (1,2) ) **2.D0+ 

C ( SRFA (13,13,3) -SFOC (1,3) )**2. DO) 

BETA-SFOCT (2,2) -SFOCT (1,2) 

GAMMA-SFOCT (2,3) -SFOCT (1,3) 

DELTA-SFOCT (2,1) -SFOCT (1,1) 

ALPHA-SFOCT ( 1 , 1) **2 . DO+SFOCT ( 1 , 2 ) **2 .DO+SFOCT ( 1 , 3 ) **2 . DO- 
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SFOCT ( 2 , 1) **2 . DO-SFOCT (2,2) **2 . DO-SFOCT(2 , 3) **2 . D0- 
C**2 .DO 


GCOEF(l) 
GCOEF ( 2 ) 
GCOEF ( 3 ) 
GCOEF (4) 
GCOEF (5) 
GCOEF (6) 


=4 .DO* (BETA* *2 ,D0-C**2.D0) 

*8 . DO*BETA*GAMMA 

=4 .DO* (GAMMA* *2 .D0-C**2 .DO) 

-4 . DO* (ALPHA*BETA+2 . DO*C**2 . DO*SFOCT(2 , 2) ) 

-4 . DO* ( ALPHA*GAMMA+2 . D0*C**2 . DO*SFOCT (2,3)) 
=ALPHA**2 . DO-4 . D0*C**2 . DO* (SFOCT (2 , 1) **2 . D0+ 
SFOCT ( 2 , 2) **2 . DO+SFOCT (2,3) **2 . DO) 


GCOEF (7 ) =4 .DO* (C**2 .D0-DELTA**2 .DO) 

GCOEF (8) =-4 . DO* (ALPHA*DELTA+2 . D0*C**2 . DO*SFOCT (2,1) ) 

GCOEF (9) =-8.D0*BETA*DELTA 

GCOEF ( 10 ) =-8 . D0*GAMMA* DELTA 

RETURN 

STOP 

END 
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